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KEY WORDS: Computer analysis; Dynamic characteristics; Dynamic 


Structural analysis; Field tests; Photogrammetry; Static tests; Structural 


ynamics; Testing; Towers; Transmission lines 


ABSTRACT: A series of full- scale structural static and dynamic tests were erformed to 


characterize selected properties of a 1200-kV mechanical test line. Prior to the 1 
experiment tests, a finite-element analysis of the test tower was conducted. The 
analytical techniques used to determine the static and dynamic behavior of a 7 i 
transmission tower will provide satisfactory results. Full-scale dynamic testing using the © 
twang excitation technique proved adequate for characterizing the dynamic properties 
of a transmission line tower and conductor-bundle system. The > importance of a aes 


4 computer analysis to aid in experiment design is demonstrated. 

' REFERENCE: Kempner, Leon, Jr. (Struct. Engr., Bonneville Power Administration, 
P.O. Box 3621, Portland, Oreg. 97208), Stroud, Richard C., and Smith, Strether, 3 
“Transmission Line "lel te Brie Structural Testing,” Journal of the Structural 

Division, ASCE, Vol. 


- 16542 SHEAR STRESSES IN CIRCULARLY ree BARS 
WORDS: Finite elements; Formulas pw’ 
(mathematics); Numerical analysis; Shear stress; Stress analysis; Stress = ~~ 
stress analysis of circulary curved bars subjected to transverse forces in the plane = 4 
curvature. The formulation follows the semi-inverse method of Saint Venant, using a_ I 
stress function. Numerical results are obtained with a finite element method, which 
accepts curved bars with complicated cross-sectional geometries. 


Md. 20810), and Pilkey, Walter B., “Shear Stresses in Circularly Curved Bars,” 
of the Structural Division, ASCE, Vol. 107, No. ST10, Proc. Paper 16542, October, 


‘REFERENCE: Thasanatorn, Chirasak (Research Scientist, Hydronautics, Inc., Laurel, i 


16543 WIND LOADS ON EAVES OF LOW BUILDINGS AGS Ors 

KEY WORDS: Aerodynamics; Buildings; Canopies; Codes; Coefficients; — 
Models; Pressure; Roofs; Standards; Wind; Wind pressure 
ABSTRACT: Wind pressure loads acting on 1:12 sloped perimeter eave sections .- 
low-rise” building models have been measured in turbulent boundary layer flow, 
simulating wind over open country and suburban terrains. Both local and area load rt 
measurements have been carried out for two eave widths (5 ft and 10 ft in full scale) : 

| and three eave heights (16 ft, 24 ft, and 32 ft in full scale). Data show the + 


eaves are derived from | ‘simple ices pressure measurements on upper and lower eav £9 
Surfaces, and full correlation of the worst cases is assumed. A simplified chart — 
Pressure coefficients appropriate for design i is also suggested. 


_ REFERENCE: Stathopoulos, Theodore (Asst. Prof. of ae Centre for Building of 
., Quebec, Canada H2G ne ‘Wind Loads on Eaves of Low a 
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16592 SEQUENTIAL GEOMETRIC OPTIMIZATION 


KEY WORDS: Configuration; Geometry; Iteration; Minimum weight design; e 
Optimization; Stre Stress rati ratio; Structural analysis; Structural design; Topology; wu 
ABSTRACT: A two-phase sequence of weight minimization is —e an 
approximate method for optimizing truss geometry and topology in preliminary design. 
In the first phase, the initial geometry is held fixed while member sizes are optimized. 1 
In the second phase, the nodal coordinates are optimized. The advantages of the 4 
heuristic method lie in reducing the number of structural analyses and in the control it 
offers the designer. Through the proposed technique, a reasonable initial geometry may © 
be obtained which allows optimization of the topology of the structure. =” 
REFERENCE: Lev, Ovadia E. (Mrg., Dept. of Advanced Design and Development, 
Merritt CASES, Inc., Redlands, Calif.), “Sequential Geometric Optimization.” Journal 
of the Structural Division, ASCE, Vol. 107, No. ST10, — Paper 16592, October, 


WORDS: Buildings (codes); Composite structures; Concrete (precast); 
Concrete (reinforced); Cracks; Interfaces; Loads (forces); Reinforcement 
(stractares); Shear strength (concrete); Shear tests 


ABSTRACT: Cyclically reversing and monotonic. transfer tests for 
companion intitially cracked specimens of normal weight and lightweight reinforced 
concrete. Both composite and monolithic specimens were tested. The interfaces of all 
composite specimens were roughened as required by sec. 11.7.9 of the ACI Building 
Code, ACI 318-77. In some cases the bond at the interface was deliberately broken. In 
the case of monolithic specimens and composite specimens in which good bond is 

obtained at a roughened interface, the strength under cyclically reversing ‘shear i is about © 
80 percent of the shear transfer strength under monotonic loading. pins te Cae aor Si 
REFERENCE: Mattock, (Prof. of Civ. Engrg., Univ. of Washington, Seattle, 
Wash.) “Cyclic Shear bem and Type of Interface,” Journal of the Structural 
Division, ASCE, Vol. 107, No. Proc. Paper 16563, » October, 1981, Pp. 1945-— 


16600 CAPACITY OF CIRCULAR TUBULAR JOINTS ‘ 
KEY WORDS: Connections (joints); Equations; (mathematics); L 
Offshore structures; Specifications; Tube joints; Ultimate strength; Ultimate — able 


ABSTRACT: A total of 137 ultimate strength tests on simple T, Y, DT, and K tubular 
joints is used as a basic for development of new ultimate capacity formulas. The data ; 

are taken from a variety of sources, and only relatively large ‘geometries are 
considered. Axial tension, axial compression, in- plane bending, and out-of-plane 
bending loads are represented. The failure condition is taken as the minimum of either 
maximum load, first crack load, or load at an excessive deformation limit. Several - 
_ formulas are recommended to predict the capacity of the different joint and load types. 
| The accuracy of these formulas is then studied in a statistical manner. Predictions of 
past and present API RP 2A formulas are compared to the same data base. The new 


‘relatively simple le equations are more consistent in their level of prediction and ‘Tesult in 


REFERENCE: Yura, Joseph A. (Prof. of Civ. Engrg., The Univ. of Texas at t Austin, — : 
Austin, Tex.), Zettlemoyer, Nicholas, and Edwards, lan F., “Ultimate Capacity of 
_ Circular Tubular Joints,” Journal of the Structural Division, ASC , Vol. 107, No. 

ST10, Proc. Paper 16600, — PP. 1984 
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i 16594 CABLE REINFORCED INFLATABLE STRUCTURES a 


KEY WORDS: Cables; Computation; Deflection; Dynamic ri response; Finite 7 
+ elements; Inflatable structures; Loads (forces); Models; Snow; ‘Structural 
dynamics; Vibration; Wind — 


4 _ ABSTRACT: The influence of internal pressure and static loads | on the dynamic 
| Saas of cable-reinforced inflatable caps is studied numerically. The — 
is modelled as a three-dimensional cable-reinforced network with the action of the 
membrane neglected except for transmitting pressure load to the cables. For ea 


small static loads, including snow and dead weight, the natural frequencies are 
| Sgpeecimanety proportional to the square root of the internal pressure. The dynamic 
response of the structure under the action of static wind load differs considerably from 
that without such loading. The static wind-type load causes the natural frequencies to 
REFERENCE: Vinogradov, oO. G (Asst. Prof.,. Dept. of Mech. Engrg, Univ. of 
| Calgary, Calgary, Alberta, Canada), Malcolm, D. J., and Glockner, P. G., “Vibrations — 
4 ——— Reinforced Inflatable Structures,” Journal of the Structural Division, ASCE, vf 
107, No. ST10, Proc. Paper 16594, October, 1981, pp. 1985-1999 


OUTRIGGER-BRACED TALL BUILDINGS 
KEY - WORDS: Buildings; Cores; Deflection; Drift; Flexibility; Outriggers; 
eee Structural analysis; Structural design; Tall buildings; Wind 


_ ABSTRACT: The behavior of outrigger-braced tall building structures is studied oie - 


- jnto account the flexibility of the outriggers. Expressions are developed for the core 
moment distribution, the top drift and the optimum location of the outriggers for 
_ minimum top drift. A characteristic non-dimensional parameter is developed that takes 
_ into account the core-to-column and core-to-outrigger flexural rigidity ratios. Graphs — 
are plotted as functions of the characteristic parameter which allow the determination 
of the optimum outrigger levels for minimum drift and the corresponding core base 
-momnet for structures with up to four outriggers. The study is based on the 
assumption of a uniform structure and uniform horizontal distributed loading. 
REFERENCE: Smith, B. Stafford (Prof., Dept. of Civ. Engrg. and Applied Mechanics, | 
Mc Gill Univ., Montreal, Canada), and Salim, Irawan, “Parameter Study of Outrigger- 
Braced Tall Building Structures,” Journal of the Structural Division, ASCE, Vol. 107, 
No. ST10, Proc. Paper 16566, October, 1981, pp. 2001- -2014 


KEY WORDS: Finite "Optimization; | 


simple case of the finite element method using plane stress, constant strain triangular | 
elements. Following earlier work, the method developed neglects the constitutive — 
equations and is formulated using only node equilibrium and a “stress-type” objective 
function from plasticity theory. A form of constitutive equation subsequently appears 
through the Kuhn-Tucker conditions. As a result, the optimization algorithm used 


takes on the appearance of an analysis- redesign procedure. Two plane stress, —_ 


REFERENCE: Spillers, William (Prof., Dept. of Civ. Engrg., 
‘Polytechnic Inst., Troy, N.Y.), Optimization: Finite Element Example,” 
Journal of the Structural 107, No. ST10, ‘Proc. Paper 16 16598, 
October, 1981, pp. 2015-2025 
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Way 


| accordance with the October, 1970 action of the ASCE Board of Direction, which stated 
that all publications of the Society should list all measurements in both U.S. Customary and 
. SI (International System) units, the following list contains conversion factors to enable readers 
to compute the SI unit values of measurements. A complete guide to the SI system and its | 
= has been published by the American Society for Testing and Materials. Copies of this 
_ _ peblication (ASTM E-380) can be ased from ASCE ata price of $3.00 each; orders must 
= - All authors of Journal papers are being asked to prepare their papers in this dual-unit format. 
_ To provide preliminary assistance to authors, the following list of conversion factors and om 
_ are recommended by the ASCE Committee on Metrication. 


yards (yd) meters(m) ia 9140 
Square inches (sq square millimeters ) eu 
feet (sq ft) square meters (m 093 

Square miles (sq miles) square kilometers (km 2.59 


cubic inches (cu in. millimeters (mm ') 400 
cubic feet (cu ft) meters (m nal 0.028 
cubic yards (cu meters (m 0.765 


, 


tons (ton) mass kilograms (kg) 
pound force (Ibf) =  newtons 
kilogram force (kef) newtons (N) 
pounds per square foot (psf) pascals (Pa) 
_ pounds per square inch (psi) = kilopascals (Pa) 


gallons (gal) Tas Bie Titers pure ak 3. 
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Leon "A. M. ASCE, Richard C. Stroud,” 
During 1978, full-scale structural- dynamic ; and static tests were 
on the Bonneville Power Administration (BPA) 1,200-kV mechanical-test line. * 
_ The objectives of the test program were to: (1) Characterize the dynamic and - 
“static properties of the 1,200-kV eight-conductor bundle and one of the suspension _ 
towers; (2) improve/ verify dynamic finite element techniques used for transmis- 
sion tower modeling; and | 3) and develop experimental testing and 


@ 


In pre for 1, transmission lines, the Bonneville Power 
- Administration has constructed a mechanical- test line at Moro, Oregon. The if 
test facility, shown in Fig. 1, consists of a 1.1-mile (1.8-km) test line with _ 
six towers. It is intended to provide information on the mechanical Tesponse~ 
of various ultra- “high- voltage (UHV) insulator / conductor "designs to artificial 


_ self-supporting, lattice-steel towers, shown in Fig. 2. The tower foundations 
are rock footings that consist of stub angles grouted in bedrock which produce 
rigid boundary conditions. Special equipment such as work platforms, eo a 
ladders, and a cable elevator is installed on these towers. - nS ee 

During this series, the test line was strung with single-, twin-, and 
eight- conductor bundle configurations of 1.6-in. (41.0-mm) diam aluminum-< -con- 
ductor-steel-reinforced (ACSR) conductor. Also, two 0.75-in. (19.1-mm) in - 


‘Struct. Engr., Bonneville Power Administration, P.O. Box 3621, Portland, Oreg. 97208. 
Pres. Synergistic Tech. Inc., 20065 Stevens Creek Blvd., Cupertino, Calif. 95014. 
*Vice Pres., Synergistic Tech. Inc., 20065 Stevens Creek Blvd., Cupertino, Calif. 95014. 
Note.— Discussion open until March 1, 1982. To extend the closing date one moth, 
a written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on March 12, 1980. 
This paper is part of the Journal of the Structural Division, Proceedings of the sre 7 
_ Society of Civil Engineers, OASCE, Vol. 107, No. ST10, October, 1981. ISSN 


| 
| 
j 
- 
| 
6 
Lhe transmission structures of the test line are moditied douhbie-circuit 
| 
| | 
— 
f 
] 


overhead ground wires lightning protection and two 0. 0.5: -in. (12. 7- -mm) 
_ diam steel messenger cables supported instrumentation lead wires in the spans. — 


SPAN 1 2 3 SPAN 4 


425 


(434.34 m) 


FG. 1. Test Line Profile 


: 


er 4 


The eight- conductor bundle configuration, shown in Fig. 3 has a bundle diameter 
— of 3.5 ft (1.07 m) with a ‘16 in. (406 mm) spacing between each conductor. 
The bundle is suspended from the tower with V- “string. _ insulator assemblies. 
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are on west bundle at spacing. Spee Number 
4 is with accelerometers mounted on the and 


The on-site data- -acquisition system incorporates a Digital Equipment Corpora- 
tion (DEC) PDP-11/35 minicomputer equipped with a 256-channel multiplexer- 


“The p purpose of these tests was to characterize the dynamic properties of — 


- Tower Number 4. The experimentally determined dynamic mode shapes and 
a ff requencies were compared to the theoretical values obtained using the Structural — 
_ Preformance Analysis and Redesign (SPAR) finite- element computer program 
(Engrg. Information System, Inc., 5120 ) Campbell Ave. , Suite 240, San Jose, 
Calif. 95130). The tower coefficients were obtained from 

Oe 


‘The physical nature of the transmission- line system restricted the dynamic 
excitation alternatives for testing. g. The classical excitation procedures were 
a considered, i.e., sine dwell, sine sweep, fast sine sweep (chirp), random, impulse, — 
release from initial displacement (twang), and ambient. The twang-excitation 
method was selected because of the capability of applying large loads at several 
locations high on the tower. The twang test is performed by loading the structure _ 
toa desired level and then suddenly releasing the load. This process approximates 
a step function excitation, which has a broadband frequency spectrum. 
; _ Intheory, all dynamic modes should be stimulated by a single-point excitation. . 
However, the existence of modal points and other practical considerations == 
_ as structural nonlinearity, data-acquisition resolution, and data-analysis dynamic- 
‘range restrictions require that tests be performed with loads at several oneal 
a! _ and in several directions. The analysis of each dynamic mode can be made 
io from a test that produces the most useful response of that mode. Four tw : 


TABLE 1 .—Tower Natural Frequencies, in Hert 


A 


dongitudinal 
Longitudinal 


hed 
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STATIC STRUCTURAL TESTING = 
Two of the tests were performed by applying a load normal to and parallel 
i with the direction of the test line. These excitations were primarily intended 
to excite the transverse and longitudinal tower dynamic modes, respectively. 
aa remaining experiments were performed by onyns a load t to the tip cd 


was applied for these tests. This load produced static stresses well within the 

A pretest analysis of the tower was performed to select instrumentation and _ & 
fs locations and to give predictions of load and response levels. Three 
finite-element models" of Tower Number 4 were developed the SPAR 


program. 1 


and the tower-leg bracing. There are approx 502 members aul 177 node points — 

in each of these models. The four base nodes of the tower were assumed — 
An eigensolution was performed to estimate the resonant frequencies and — 
} associated mode shapes. The natural frequencies of the first five tower modes 


are shown i in Ta Table In general, the modified f rame and 


Structural 
damping 


, in Hertz in Hertz al 


Ist Transverse 2.37 0. 0160 


2nd Longitudinal | 264 © 37 0.020 
were higher than the frame results. The: fundamental mode found with the frame > 
analysis was not detected with the space truss model. 
_ The eigensolution used the diagonal lumped mass matrix to obtain the results — pe 
- in Table 1. In addition, the consistent mass matrix was used to obtain the 
natural frequencies of the truss model. _ These frequencies were approx 3% % higher 
than those shown i in n Table 1. me 


> "truncated to contain only the yey -of- freedom at the instrumentation positions. 
_ The mass matrix used in the orthogonality check was determined by a Guyan = 
() reduction of the full SPAR- -generated mass matrix. This analysis was used 


is approximately the unit "matrix, i.e., ones along: the diagonal and small” 
off-diagonal terms, which demonstrates the instrument distribution to be satisfac- 
tory for characterization of the analytical modes. 
The final step in the pretest analysis was to compute the dynamic response - 
_ _of the tower to each of the planned loads. This analysis was performed using — 
the Guyan-reduced mass and stiffness matrices. god 
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FIG. 5.—(a) First Longitudinal Tower Mo 
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STATIC STRUCTURAL TESTING 
and using least-squares techniques to determine a set of linear, superimposed _ 
modes which best fit the data. The Vibration n Analysis and Measurement Processor 
(STI-VAMP, Synergistic Tech. Inc. , Cupertino, Calif. ) was u was used for all data 
reduction and display. The result of the analysis was a set of mode shapes, 
modal frequencies, and damping estimates, 
Table 2 summarizes the modal characteristics of the first three tower modes 
_ determined by the theoretical and experimental analysis. The resonant frequencies 
‘determined by the theoretical analysis are approx 10% higher than the experimental - 
| determinations. Refinement of the tower model to include the effect of the 7 
_ conductor and insulator assemblies would lower the theoretical values. Other 
- - factors, which would cause discrepancies in these values, include flexibility 
and —_— of the tower joints which can not be accurately modeled. Curvefit — 


- “a experimental mode shapes. These results are for the first three tower modes. 

TABLE 3.—Tower Dynamic Response, in Gravitational Acceleration 


Degree o 


| 


Transverse 
Transverse 


Vertical 


figure represents measured degrees-of-freedom with significant wanes 
Correlation of all the and theoretical modes \ was determined 


Table 3 is a comparison of some of the measured and theoretically predi — 


acceleration responses for the tower. 


- Static load / deflection tests were performed prior to each tower dynamic test. 
The load was incremented to 10 kips (44.5 kN) in two cycles with a total of 


3 18 load steps. Deflections were recorded using Wild T-16 theodolites and two 


 close-range photogrammetric MK70 Hasselblad cameras. Fourteen theodolite 
and 43 photogrammetric targets were bas beoi 


motions to all of the modeled degrees-of-freedom by using 32 | 


— 
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_ "dimensional displacements for each photogrammetric target. With the two c camera 
_ systems used, the measurement resolution (under ideal conditions) is +0.2 in. a 
- (+5 mm). Good results were obtained with this method when displacements 
_ were on the order of 1.0 in. (25.4 mm) and greater. Ths > reduction of this 
: Fig. 6 shows typical load/deflection plots obtained from the theodolite and : 
- poogammeri data. The stiffness values for the theodolite data, 4.6 kip/in. 
ti (797 kN/mm), and the photogrammetric data, 4.9 kip/in. (860 kN/mm), show = 
 Deflections below the upper work platform (Fig. 2) were too small to detect 
with the measurement : menmen used. The photogrammetric data showed } more 


MMETRIC 


Concentric circles 


VERTICAL LOAD kips ( 
eo | 


~04 -08 12 -16 
VERTICAL DEFLECTION inches (mm) 


scatter than the theodolite results but did d verify th the effectiveness of the approach. 
It is expected that at least a twofold improvement in accuracy can be obtained 
with improved techniques, 
_ Analytical deflection results were obtained from the three finite-element models. 
gg 4 tabulates the deflections at the load point for a maximum load of 10 — 


kips (44.5 +.5 KN). The ¢ model results show good agreement with the field measured © 


evaluate its use 
| 
= 
. 
1 load and a I-kip (4.5-KN) torsional load twang test at Spacer Numbers 2 and 


= Testtype/- | Experimental | 
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y of Span Number 4. Spacer ‘Number 2 and ne: at 26% and 43% 
2 of full span respectively, measured from Tower Number4. = = i 
_ The primary instrumentation used to characterize the bulk dynamics (low- 
frequency response) of the eight- “conductor bundle was a set of accelerometers 
mounted on Spacer Numbers 2 and 3. These tra transducers were installed in a 
radial pattern with an angle of 120° between each unit, and with the bottom _ 
_ accelerometer oriented in a horizontal position. Vertical, horizontal, and rotational _ 
motions can be separated by using this transducer configuration, = 7 
- Response time histories and spectra were analyzed for a period of 60 s, b. 
_ Starting 60 s after the initial excitation. During this period, the traveling waves, 
. which dominated early conductor motion, have been damped out and ‘only 
_“‘modal’’ behavior remains. Curvefit analysis showed that the first two modes | 
were at 0.16 Hz and 0.28 Hz and the structural-damping coefficient varied 
from about 0.08, for the low-frequency modes, to 0.004 for modes above 2 " ; 
An investigation o of transmissibility between the conductor and tower response > 
for a conductor twang test was performed. . The transmissibility was found to 
be small over most of the analyzed frequencies, but, in the range of low-frequency 
tower modes, significant energy (with a gain of up to 1) can be transferred 


from m Ge conductor to the tower. 


Conouctor Static 
- "Displacements were recorded at load increments for each of the cight- -conductor 
bundle dynamic tests. In addition to these tests, a span-to-span static interaction 


TABLE 5.—Stiffness (Conductor Vertical — 


Load point Spacer Spacer #3 


Spacer #20 | (0.983 0,699 
Space #3 


was between ‘Numbers 2 and 3 (Fig. 1). A maximum load 

of 4 kips (17.8 KN) was used forthistest series, 

aa and photogrammetric data were recorded for the vertical and 

torsional tests of am Number 4. Only theodolite data was recorded for the 

a reference point on 

the conductor nrg at Spacer Numbers. 2 and 3 to measure the vertical 

7 erect Rotation targets were attached to Spacer Numbers 1, 2, and 3 


to record the bundle twist. Groups of three orthogonally oriented photogrammetric 


targets were located on Tower | Number 4 conductor-bundle yoke plate and on - 


Spacer Numbers 1, 2, and 3. Three- dimensional translations plus rotational 


displacements were from the motion of these targets. 
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STATIC STRUCTURAL TESTING 


mers 7 shows typical values of the vertical deflection measured 2 


Number 2 with the applied vertical load at this spacer. Good agreement was -_ 
obtained between the theodolite and photogrammetric data. 


load / deflection curve, for the vertical-load tests. This data shows an apparent 
in-span reciprocity between Spacer Numbers | 
Table 6 tabulates the stiffness values, i.e., slope of the torque /angular 
curve, obtained from the conductor-bundle torsion tests. The load / deflection 
_ plots of the torsion tests showed similar linear behavior to that shown in Fig. 
7 for the conductor-bundle vertical tests. 
TABLE 6.—Stiffness (Conductor Torsional eed 


Stiffness, in foot-pounds per degree = 
be complications becasue of (Joules per degree) 


Load point Spacer #1 ay Spacer 


Spacer #3 


poi int aval Center span # cate Center span 


(67.36) 


this test are presented in 


pplied at the center of each span nper the resulting vertical displacements, 
P 


span-to-span reciprocity of +10% is demonstrated by these results. —__ 04 


7 of the transmission line tower provide sufficiently accurate results. The three 


tower! showed agreement among themselves for both the Static and 


1905 
bt 
4 
| 
q 
q 
1 
| 


and modal frequency and mode- 
_ determinations agreed within 10%. The experimental determination of damping | 
showed values from 0.015-0.040 for the structural damping coefficient. sit ht 
_ 2. This test series demonstrated the twang excitation technique to be adequate 
; ‘for characterizing the dynamic properties of the tower and conductor-bundle | 
“(or low-frequency response) system despite data analysis difficulties that arise be 
_ from short-duration response. A potential improvement may be to provide a 
continuous excitation function that can be controlled in level and duration. 
_ A hydraulic shaker, driven by a sinusoidal or chirp (fast sine sweep), appears 
to be a potential solution. This excitation method would also allow for the 
= of the — as ra to a wide range of simulated ‘ “real” excitation. 
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SHEAR STRESSES IN CIRCULARLY | CurRVED Bars 


¢ By Chirasak Thasanatorn' and Walter D. Pilkey, M. ‘ASCE 


The exact stress determination of a circularly curved bar subjected to pure > 
aid bending and shear forces: in the plane of curvature requires the solution of | 
a boundary-' value problem of three- dimensional elasticity. This presents formida- _ 
ble complications because of the difficulty of satisfying the boundary conditions. 
Few elasticity problems for curved bars have been solved. Analytical solutions — 
of a curved bar with circular cross sections were developed by Golovin, as 
presented in Timoshenko (2). Also given in Ref. 2 are plane elasticity solutions — 
- for curved bars of small width. A so-called strength-of-material theory of curved 


_ bars was developed by Winkler using assumptions analogous to the Bernoulli-Euler _ 
hypothesis of straight beam theory. In this case, the equilibrium conditions | 
are satisfied in an average sense} 
This paper presents a method for the stress analysis of a circularly — 
bar subjected to a transverse end force in the plane of the bar curvature. This 
formulation follows the semi- inverse method of Saint Venant, using a stress 
function. The assumption that stresses on any cross section depend upon forces 
at that particular cross section permits the three-dimensional elasticity — 
to be reduced to determining the stresses on the cross section of the curved | 
_ bar. The s stress solutions obtained from this structural model may be eer 


“that the ¢ cross section of interest is sufficiently | far from any - points of -_ 
variation in the internal forces. Finally, a finite element method, which is a 
for a curved bar with complicated cross-sectional geometry, is used to determine | 
the stresses on the cross 
“ Torsional stresses for circularly curved bars are treated in I ef. =o 


Ww 
7 Collier i a planar curved bar of initial curvature » Ry v with one end fixed and 


s other end free, as shown in Fig. 1. The bar is subjected to a transverse 


‘Research Scientist, Hydronautics, Inc., Laurel, Md. 20810. 
Prof. of Applied Mechanics, Univ. of Virginia, " Charlottesville, Va. : 22901. eat 
_ Note.—Discussion open until March 1, 1982. To extend the closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, _ 
ASCE. Manuscript was submitted for review for possible publication on March 4, 1980. > 5 
‘4 This paper is part of the Journal of the Structural Division, Proceedings of the American -_ 
Society of Civil Engineers, ©ASCE, Vol. 107, No. ST10, October, 1981. ISSN 0044- 
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radial force passing through the e shear center ; at the free end. Assume: (1) The . 
_ material i is linearly elastic, continuous, homogeneous, and isotropic; (2) the curved 
HX bar is planar, lying in the x-z plane with the longitudinal axis coinciding with 
the geometric centroid of the cross section; (3) the curved bar can have any 
4 arbitrary cross-sectional configuration but must be prismatic; (4) the body forces © 
are small compared to the magnitudes of ‘the stresses and can be neglected; 
6) the stresses at any cross section are assumed to depend upon only the 
internal forces at that cross section; and finally (6) the normal stress distribution 
& assumed to follow a hyperbolic law in the radial direction. The curved bar 
is modeled with fixed-free ends to simplify the mathematical formulation. For 
this case the variations of internal forces along the curved bar axis | are either 


or are a function of sine or cosine functions. 


(y, 2) = LOCAL COORDINATES WITH ORIGIN: poten 
THE CENTROID OF THE BAR 
FIG. 1.—Structural Model for Direct Shear Stress Analysis 
_ From the overall “equilibrium conditions of the t bar of Fig. 1, the ‘internal 


7 forces at any cross section along the bar axis are 
W sin sin®; M, = —WRsin®; V,= Weos 

‘The twisting moment, M,, vanishes due to the assumption 1 that the force W 
passes through the shear center of the cross section. 
4 From the aforementioned assumption 5 and Eq. 1, the three-dimensional 
elasticity problem may be reduced to a two-dimensional one by expressing the — - 


o, = &,(y,z) sin 8; =5,(y, z)sin®; = 


| 
= 
< 
|g 


in which Tx, and F = functions of cross-sectional ‘coordinates 


—( y,z) only. Assumptions 5 and 6 suggest trying to satisfy the stress equilibrium 
equations of assuming a normal stress distribution in the 


R Cy + 2] ‘sin® . 

which C, = constants to > be “determined. By substituting Eq. nid 
3 into the stress equilibrium equations, the remaining stress components can 
g derived. These stress components are expressed in in terms of the two stress 4. 


UR 


ts 
[F(z) YI[C, + 2C, , F(z)] sin 


(R + @ +z) 2 2(R+ 


 2RG 

— + Cy” 


R? 


7? 


+ 2) + C,(2 ~ 


in which F(z) and G(y) are defined along the cross-sectional boundary. — a 
7 The stress expressions in Eqs. 4-9 are derived independent < of the displacements _ 


a, ¥, and w. The stress components, or more specifically the stress functions 
and YY, cannot ‘be taken as arbitrary functions of (6,y,z) but are 


Co + C, [G(y) — 2R] ¢ } cos 


= | 
| 
o 


_ to the two stress compatibility conditions to ensure the existence of maint : 
u, v, and w corresponding to these stresses. A set of consistant stress- -displacement 
relations can be derived by using the Hellinger-Reissner variational principle, 
which is a mixed principle. Then, two independent stress compatibility conditions 
can be derived from the five consistant stress-displacement equations. They 


Cy + C,y 


r 


e- 


in which r = :(R i i“ Cartesian and cylindrical coordinates a are used interchang 
ably as convenience dictates. ™ constants Co, C,, and C, in the assumed 


+ 


conditions for any section of the curved bar. 


“20, dydz= WR sin 6 


onstants C,, C,, and C, are found to be 


which J, + 2/R)| dydz; 1, = (y?/( + 2/R)) dydz; 
= + 2/R)) dydz. By substituting Eqs. 4-9 into 10 and 
the two governing Partial Gferential equations for the stress functions o and 
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(16) 
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for the s stress functions along the cross-sectional boundaries. 


‘The governing partial differential equations given in Eqs. 16 and 17 are too 


a ‘The governing equations for the stress functions, @ and W, as given in Eqs. 
16 and 17 are too complex to be suitable for closed form solutions. A variationally 
based finite element procedure will be applied to obtain the solutions. The 
problem of integrating the partial differential Eqs. 16 and 17 subject to boundary os a 
conditions of Eq. 18 may be transformed into the problem of finding the unknown a 7 
stress functions @ and WY which make stationary a functional subject to the © 4 q 


_ same bai conditions. The functional may be a in matrix form 


m= 55100)" [c1@) - 


ee 


| al surface of a curved eo: 
| 
ed to solve for the two Unknown stress Tunctions directly. ihe finite 
element formulation will be derived directly from a stationary condition mee 
= functional for which Eqs. 16 and 17 are the Euler conditions. 
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—(C, + C,2z)-v ——c - 


+R-z 
2(R+2) (R+z) (R+z) 
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- C, F(2)[F(2) 
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2G(R+z) 
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(3 +2) te 


26 (R+2) 


_ The procedure will be to idealize the cross- -sectional domain by an ni samme 
“of triangular elements. . The functional Tl may be ; approximated as the sum of ; 
functionals evaluated for — element. Mathematically, the functional may nay, be, 


ment num 


ber; and M= . the total number of elements. . 
At “node j the stress functions and their first derivatives (slopes) with respect — 


er lebom 
q 
4 


» and z are chosen as the unknown nodal parameters , and are expressed 


25) 


- 


The st stress — a iin each element m 118 be written in terms 


of their nodal parameters 


in which wi =a shape f function matrix and {®, v= a vector r of th the unknown 
iter parameters. The vector of derivatives in Eq. 20 may be — by 


differentiating | with respect toy and z, Thus 


4 in which elements ¢ of —_e By" are the second derivatives of the shape functions — 
Be respect to y and z as defined in Eq. 20. The discretized functional ae 


element is terms of the unknown by 


integration is performed over the element domain Finally t the total 

- functional is obtained by collecting the contributions from all discretized func-_ _ 


The stationary solution of the functional IT can be derived from the condition _ 
: is that the variation of the functional I(y,z,,) with respect to the unknown 
nodal parameters vanishes. That is 


in in which N = the number cof nodes the. cross 


fi formulation wi will be ap to two bar « cross sections. Fii First, a 
homogeneous thin curved bar of rectangular cross section will be used so that 
the calculated stresses can be compared with an analytical solution available 

from plane elasticity theory. The second example is chosen to be a curved 
_ Homogeneous Thin Curved Bar of Rectangular Cross Section.—The thin rect-— 

“ angular curved bar of Fig. 2(a) has a | in. (0.0254 m) width, 10 in. (0. 254 od 
depth, and a 20 i in. (0.508 508 cross — 
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CIRCULARLY CURVED BARS 
smo an 


is to an internal shear in the radial of 10, 000 Ib (44.48 
89 at its centroid. The discretized model of the cross section is shown in a : ; 
‘Fig. 2(b) and partial computer results are presented in Table 1. The stress - 
‘results obtained by the present formulation are compared in Table 2 and Fig. : 
i: * with those obtained from the plane elasticity theory. The direct shear stresses - 


B.- on the present formulation are the average values across the bar width. — 


(0. 254 m) 


4 i 


+ 


F CURVED BAR 


— 


FIG. 2.—Thin Curved Bar Subjected to Transverse Shear: (a) Dimensions; (b) ica ze 


CENTER 


complicated cross-sectional configuration. Choose the curved bar section = 


Fig. 4. It is made of steel with a modulus of elasticity and Poisson’s ratio = 5 


* of 30 x to" psi (206. 85 GN/ m’*) and 0.3, seapoctively. The section is subjected 


= demonstrates the applicability of the method of analysis to a more 
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—Shear Stresses due to Shear Forces (in pounds per square inch) on panes 


TABLE 1.— 
Curved Bar Section of Fig. 2 mon 


= — 
24.640E + 00 
—27.859E + 00 
2 859E + 00 
+e 
-22.034E + 00 
~21.587E + 


— 


+00 


—22.034E + 00 
18.585E +00 
—18.192E + 00 
> 18.192E + 00 
=18.585E + 00 
15.448E + 00 
+ 00 
+ 00 

=15.448E + 00 
+ 00 
-12.148E + 00 

12.148E + 00 
| + 00 
7 99. 201E — 01 
9ISE O01 

94.91SE — 01 


188 ~27.049E 01 


188 


on S81E 02 


in. = 0.0254 m; 


(-12.189E 01 


467E-O1 | 


rh 


va 


19.225E + 01 
19.225E + 01 
24.335E + 01 
 61.280E + 01 
$7.500E +01 
57.500E + 01 
 61.280E + 01 
90.825E + 01 
88.061E + 01 
88.061E +01 
90.825E + 01 
11.382E + 02 
‘11.192E + 02 


‘1L.192E + 02° 


+ 02 


| ‘13.105E + 020 


| 


15.087E + 02. 
15.087E + 02° 


| 087E + 02 


O87E 

+ 02 
1S. SO2E + 02 
4 457E + 02 
15.474E + . 
02 
+ 02 

IS.178E + 02 
15.294E + 02 
>: 15.294E + 02 


15.502E + 02. 
+020 
+02 


— 
4916 ST10 
i 
7 
— 0.688 
1313 
0.313 1,688 
» 
| 0313 2188 a 
OSE + 02 
0.813 
3.813] 52.2858 OF 
0.688 


CIRCULARLY CURVED BARS 


‘bor with Yourdary condiiions amd Tiyan, the 
‘TABLE 2. ._—Comparison of Shear Stresses (in pounds per s square inch) along Center 
Line of Thin Curved Bar of Fig. 2 Obtained by Present Formulation with those of 


Two- Dimensional Plane Stress Theory 


— 


217.80 244.69 
94.4 
Ne 
1,304.70, 
(1,429.40 
1508. 


‘DIRECT SHEAR STRESSES A ARE THE AVERAGE AVERAGE ‘VALUES: ACROSS 1 me BEAM WIDTH 


FIG. 3. —Direot Shear | Stress Distributions Gomer Line of of in. 


g Vor 
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— 1523.60 | 1478 § 
OS | 59.65 4071.08 
Note: 1 in. = 0.0254 m; 1 psi=6.895kN/m* 
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« FIG. 4 4.—Finite Element Representation of Curved Box Section 


4 as 


three- theory of elasticity. The assumption at any cross 

. 3 ‘section depend only on the internal forces at that particular cross section” 

; enables the three dimensional elasticity theory to be reduced to the determination 
of the stresses at a particular (two-dimensional) cross section. — the 


| 
7 
| ¥ | 

FIG. 5.—Direct Shear Stresses on Curved Box Section 

| to an internal shear force of 100,000 Ib (444.8 kN) in the z direction. The ; 

| finite element idealization of the cross section is shown in Fig. 4. The shear 7 
a stresses throughout the walls were computed. The average shear stresses across — 

| the wall thickness are shown in 
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GIRCULARLY CURVED BARS 1919 
{formulation can be used to compute the transverse shear stresses for a dis 
curved bar with any boundary conditions and arbitrary loading. T hus, the 
formulation requires only that the internal forces be known on the cross section __ 
| 
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Loaps on on Eaves oF Low Bu 


= the : assessment of wind loads on roofs, the ‘effect of eaves es (otherwise 


called canopies or overhangs) is sometimes overlooked. Windward eaves, 
however, may be loaded severely due to wind, since the deflected flow from , 
the windward wall gives rise to a pressure on the lower eave surface, _ Which 
‘reinforces the high suction of the upper eave surface right after "separation. 3 
_ Unfortunately, only a limited amount of information is available on the subject. - 
This is the reason why most Standards and Codes of Practice have poor 
; documentation of provisions selative to wind loads on eaves; in fact, they rarely 
_ Leutheusser (5) has published a _ comprehensive study on -mean wind loads 
on eaves using models of various geometries, but he performed all his cant - 
: under uniform wind conditions (no velocity variation with height). Thus, as - 
test conditions did not simulate natural wind effects, and therefore his results| 
_ may not be representative. Dreher and Cermak (4) have examined local pressures < 
on upper and lower eave surfaces it in a boundary layer wind tunnel we cre 


configuration only, but they have reported that pressures ‘underneath the o7 over- 
hangs can account for 50% of the total wind load acting on the overhangs = 

of a house . In contrast, Tieleman and Reinhold (10) have concluded that a 

presence, or absence, of eaves does not influence the magnitude of the wind — 

_ ‘pressures a great deal, i.e., the same pressure loads may apply on roofs with — 
or without overhangs. Finally, Best and Holmes (2) have examined local pressures /- 

on eaves of a low-rise building model and have recognized that combining the 

negative peak values on the top with the positive peak values under the eave 

_ is somewhat conservative, since peak values do not occur together. — ae 


=. The lack of pertinent information regarding wind loads on eave sections, 


‘Code and Standard Committee Officials to significantly - increase wind design a 
‘specifications for eaves. Ref. 6, for instance, states | that the Building Code 
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Age 


Japan increased the wind for the design of eaves in 1973. 
No research work, however, appears to justify such a decision. — ia SP ey 
bas _ Itis the purpose of this paper to describe the results of some recent experimental 


measurements is of n mean and fluctuating wind loads on eaves Of | low-1 -rise buildings 


vind 

with tributary a areas of structural elements has been demonstrated. All measure- “a 


ments have been carried out for a variety of representative low-rise building 
in turbulent flow conditions characteristic of the natural wind. Two 


terrain types and all critical wind directions have been considered. The analysis 


Of | the experimental data has finally led to some eee: for parse 
standards and Codes of Practice. 


4 


73 of Model with Large Eave Section ( (Root | Corner) 
a : elias section of about 80 ft (24 m) long, 8 ft (2.5 m) <a and about 
7 ft (2 m) high. Most of the fetch is required for the natural production of 


a a boundary layer which grows in a manner paralleling the atmospheric process 
_ under neutral conditions. A detailed description of the tunnel has been given 
in Ref. 3. Two terrain models have been used for the measurements reported _ 
_ here, representing open country and suburban (built-up) full-scale conditions, — 
- respectively. The velocity and longitudinal turbulence intensity profiles measured a, 
. at the test section without the model in position are representative of exposures 
A and B of the National Building Code of Canada (7) or C and B of the 4 
American ‘National Standards Institute (ANSI) Standard (1). 
_ The pressure measuring system used responds to pressure fluctuations on cl 
> the model of up to about 100 Hz with negligible attenuation or distortion. Several *. 
pages measurements can be made in parallel, each of which is sampled at 


+ 
| 
— 4 
_ Laboratory of the University of Western Ontario. This well-known tunnel a | 
> 
| 
— 4, ZZ 
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30 in time which a for input, the maximum 
and minimum values that occur, and computes the mean and root-mean-square > 

_ values. At the end of the sampling period, the measured pressures are converted = 
to pressure coefficients by dividing each by the reference dynamic pressure. 


The latter is measured in the free stream above the boundary _ where the 


at a scale of 1: 250) | a 1:12 ‘roof-sloped building, 80 ft (24 wide and 


ft, 24 ft, and 32 ft (5 m, 7.5 m, and 10 m) in full scale. Eaves of two different 


wide prt i 5 ft: and 101 ft (1.: 5 m and 3 m) in full scale were constructed 
“Wi 


: 3 (38 m) long. Three different eave heights were tested corresponding to 16 


: 


18+ 
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2.—Exploded Plan View of 3.—Exploded Plan View of 10-ft 
Eave Model Showing Pressure Eave Model Showing Pressure 


ing Locations (Note: 1 ft = 305 ing Locations 1 fe = 0. 305 m) 
= 
as extensions of the roof surface perimetrically. The eave sections were equipped © 


with pressure taps located | on both upper and lower surfaces. ‘Fig. 1 shows 

a picture from the top | of a corner eave section of one of the models. The 

holes for upper and lower surface points of measurement have been drilled — 

on the eave sections after drilling the path for the transmission of air into 

_ the tubing and pressure transducer. Because of the small thickness of the eave, 

- this path for each pressure tap had to be drilled from the edge of the eave 

at an angle of 45°. The thickness of the eave was subsequently | restored by | 

plugging up the small part of the path between pressure tap and eave perimeter. 
Figs. 2 and 3 show the pressure tap arrangement for both upper and lower — 


of the models with smaller and eaves, The former 
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model was equipped with 44 pressure. tops ps (22 « on pare surface), , whereas the "7 


latter had 58 taps (24 on the lower and 34 on the upper surface). As can . 
= the i sora taps have been | arranged in in such am manner or measurement 


_ close proximity of pressure taps on upper and lower surfaces serves s the perposs eo 
_of monitoring the total wind load acting on the eave at any instant. = oe 

a Experimental data of local pressure measurements carried out for the small = 
eave section [(eave width/building width) = 0.0625] have been summarized _ 
in Figs. 4 and 5 for upper and lower surfaces, respectively. Results are presented | 

in pressure coefficient form referenced to the dynamic velocity pressure at 
eave height. This height has been chosen because it was found that for low-rise 
this values & sensitive to to building height. 
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FIG. 4.—Typical Range of Worst Max, q » 5.—Typical Range of Worst Max, 

_ Min Mean, and Rms Local Eave Pree Min, Mean and Rms Local Eave Pressure — 
sure Coefficients (for all coy Coefficients (for all Azimuths)—Eave | 


q Eave Width = 5 ft (1.5 m), Upper Eave Eave ‘Width = =5ft (1.5 m), Lower Eave Sur- 


The measured pressure coefficients (referenced to gradient height) had thus 
to be converted to the new reference height, and this was done by using the 
_ Gaperimentelty determined variation of the wind speed with height. Table 1 | 
_ shows the relationship of wind speeds at eave height and gradient height for af 
the two ‘(terrain Pressure coefficients are defined as follows: 
_ in which Prax = the maximum instantaneous pressure measured over the sampling 7 
period; P xin = the minimum instantaneous pressure measured over the sampling 
period; aot = the temporal mean pressure; P,,,, = the root-mean- in-square pressure 
[(P— } 1/2 te = the dynamic velocity pressure associated 
with the mean wind speed, V,,,, at eave height; and all pressures are differential % 
pressures with respect to the static pressure at gradient height. a @ Snr € 24h 
The approximate range of local pressure coefficient values is caer 


™ 


4 
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5° intervals), two terrain 
roughnesses and three building eave heights. By virtue of symmetry, only one | 
quarter of the perimeter eave is necessary for the most critical pressure coefficient 7 
“values. Zones | and 2 appearing on the overhangs have been determined somewhat | 
subjectively and their boundaries should be considered very approximate. Fig. 
_ 4 shows that zone | is also larger for the maximum instantaneous pressure 


coef ficients than for the minimum mean and root-mean-square pressure coeffi- a 


the more uniform distribution of the lower eave surface pressure coefficients _ 
(Fig. 5) in comparison to the upper eave surface coefficient distribution (Fig. 

TABLE 1.—Ratios of Wind Speeds at Eave Height and Gradient tot 


te acting on the upper eave surface. It is interesting to note, however, : : 


ft=0.305m) | Open-country exposure Built-up exposure 


e's 


SSS 
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FIG. 6.—Typical Range of Worst Mex, FIG. 7.—Typical Range of Worst Max, | 
_ Min, Mean and Rms Local Eave Pressure Min, Mean and Rms Local Eave Pressure ; 
Coefficients (for all Azimuthe)—Eave 


Pe Coefficients (for all Azimuths)—Eave 
Width=10 ft (3 m), Upper Eave Surfaces Width = 10 ft (3m), Lower Eave Surfaces 


(8). This difference appears to be due primarily to the pressure taps being 


Lo ‘slightly farther from the edges than in the no-eave case. Obviously, some flow — 


differences also exist. Nevertheless, the present results provide representative 
— local loads on small eave sections around the roof. Mean pressure coefficients _ 
are positive for lower eave surfaces and their values are somewhat lower than os 


absolute values of mean pressure coefficie per surfaces. Also, 


ig 
i 
| 
4+). The small dependence of coefficients on building height also, as previously _ 
_ _ In general, results have somewhat lower values than were found for roof = | 
edges and corners of similar buildings without eaves exposed to the same terrain = 
the 


the rms local coefficients are for the lower eave 
for the the higher heights and the rougher exposure. baw 


This does not mean, n, however, that | pressures. are higher i in a the rougher exposure. ae 
_ In contrast, pressures are generally lower for higher terrain roughnesses because a 
of the decreased dynamic velocity pressure at eave height in which all pressure 

_ Figs. 6 and 7 show ranges of upper and lower eave surface local pressure — 
coefficients in a similar fashion for the larger eave section [(eave width / building 
— 0.125]. The number of areas for which the range of minimum and * 
mean pressure coefficients for eave upper surface is given (Fig. 6) has been 
ewe due to the rapid reduction of local pressures with distance from the 
edge. Maximum and rms values are quite uniform, so the number of areas 
has been minimized for them. Mean pressure coefficients for the outer eave 


surfaces. appear somewhat higher than in the case of the smaller eave section. 
This may be, again, due to the closer proximity to the edge of the pressure 

- taps used in the case of larger eave section or to alterations of the flow due 7 
to the big eave. Generally, the pattern of local pressures is similar for the 


‘ TABLE 2. —Extreme Local thai! Coefficients for Upper and Lower | Eave: Surfaces 


Building “Upper surface Lower Reference 
+1.62 Present study 


‘two different eave widths. For lower eave surfaces, the pressure coefficients — 

- appear a little smaller than the respective coefficient values for the smallez 
eave section. The distribution of these pressure coefficients: is, however, 
uniform as in the case of the smaller eave. 
_ Table 2 compares the extreme local pressure coefficients reported from the 

measurements on upper and lower eave surfaces. Only studies conducted in 

_ simulated atmospheric boundary layers have been considered for the — 

Despite the fact that there are geometrical and other minor terrain roughness _ 

differences | between these studies, data show a fair agreement. The higher values, ; 

of Ref. 2, however, may be attributed to the higher ‘‘eave width over building _ 


width’’ ratio. Nevertheless, more data are required to increase confidence in 


) 

| 
angle | 
§ 
| 
— 
| 
[am == smaill and large eave width. Vertical load coefficients were measured over different 
1 eave areas either separately or simultaneously on both surfaces by using the < 7 


pneumatic mmanifolding_ alone or in combination the 


Of eave sections was carried out “‘on-line’’ with a digital computer sampling 
“| one of: (1) matepantent pressures; (2) independently manifolded upper and lower 


by using both sides of a pressure transducer simultaneously. pauls  sviostts. ‘of 
_ The component loads were weighted by appropriate tributary area factors 
determined for each pressure tap as the percentage of the total area over which 
the pressure measured had been considered to be wisn 


: FIG. 8.—Eave Areas Examined (Note: Att ft = 0. 305 m) 


Typical results of maximum, mean, and minimum vertical | force coefficients 
for a half-end eave of a building exposed to open country terrain are presented 7 
_ in Fig. 9. Data are given for three eave heights, two eave widths, and all azimuths ~ 
tested. As expected, critical loads occur for 45° and 90° azimuths for the half-end d 
_ eave considered. The little dependence on height of the coefficients is also 
"apparent. The loads appear consistently higher for the smaller r eave Wi width. Similar 
trends have also been found for the built- -up (suburban) exposure. 
_ Comparison of extreme vertical force coefficients with the difference of local — 
coefficients obtained for points on the upper and lower eave surfaces (results — 
_ ees in Figs. 4-7) indicates some alleviation for area loading due to— 
the beneficial the designer’s point view) | lack of ‘Spatial 


| | 
a 


of local pressures es acting inside an area -: the upper or lower eave seutace. 


(1.5 m) exposed in “open country terrain. The solid lines of Fig. 10 represent 
the effective values measured, and these should be used for the design loads 


= that peaks of upper and lower surface area loads occur simultaneously, — 
4 matter when i in the record the peaks actually occurred. The | were found 


st 


by taking into account the true loads is in the “order of 20% for then most critical 
3 Similar results appear in Fig. 11 for the building with a 10-ft (3.0-m) wide 
_ perimeter eave. In this case, the alleviation is even higher—approx 30% for 
the most critical azimuths (0° or 90°). It should also be noted that this alleviation 
would be even higher if the true effective loads acting on eave sections were 
_ compared with the individual upper and lower surface area loads based on 
the envelope of extreme local pressure measurements; the latter used to be 
the common practice until very recently. This i is shown in Fig. 12 whose data 


= 
—suctions acting on the end-bay eave, interior bay eave and half-end eave (as ee . 
_d ned in io 8) of a 6 m) hich hnildine with an eave width of § ft 
be 
| 
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INSTANTANEOUS SPATIAL AVERAGING 


Instantaneous and Noninstantaneous Spatial Averaging—Eave Height = = 16 ft (5 
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FIG. 11.- —Peak Suction Coefficients on Eave Sections of Various Sizes by Using 


bs Instantaneous and nine gece ‘Spatial Averaging—Eave Height = 16 ft (5 


7 j 


eave corner areas and spatial averaging over these areas have been carried 
_ out. Loads have been obtained for the upper eave surfaces of both 5 ft and 
; 10 ft (1.5 m and 3.0 m) eaves, whereas total (lower-upper surface) corner loads 
have been measured for the larger eave width case only. This set of measurements 
is for an -open- country terrain. The two upper diagrams of Fig. 12 indicate 
the extreme suction coefficients measured on the upper eave surface, their 
arithmetic average (dashed line) and their effective values (solid line) which 
_have been determined by considering the nonsyachsonization of the individual 
peak values through pneumatic averaging. In both cases, the peak 


j 


~ 


PT4 


Line) and Noninstantaneous (Dashed Line) Spatial — Height = 
(5 m), Exposure = Open Country 
= 
coefficients appear very similar (somewhat higher for | the smaller eave) and 
a significant load reduction is apparent when spatial averaging is taken into | 
consideration. The lower diagram of Fig. 12 shows the effective suction coeffi- 
cients for the simultaneous wind action on upper and lower corner surfaces ‘ 
for the large eave section in a ‘similar fashion. In addition to the Snood 
due to spatial averaging, previously described, it is very interesting to. = 
i" . that the maximum total uplift coefficients are not significantly higher = 
the suction coefficients loading the upper eave surface; in fact, they have — 
‘somewhat smaller values for the critical azimuths. This beneficial result can ~ 
be enoed to the low positive pressure or suction acting underneath the eave | 


19300 
extreme € comer sections 
| 
| 

| 
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a corner ond to favorable correlation of loads acting on ‘inti eave surfaces. 
All of the preceding shows the amount of conservatism involved when the 
i upper bound load coefficient, which is derived by simple subtraction of worst 
case data measured independently for upper and lower surfaces, is considered. q 7 
F ‘Unfortunately, no comparison of the present results with other studies can 
be made, because, to the writer’s best knowledge, there is no other information : 


a simplified chart 
of pressure coefficients appropriate for. incorporation in a a design standard or 4 
a Code of Practice. This task has been undertaken previously for wind loads + 
- on low-rise buildings (9) without special consideration to eaves except some 
specifications for canopy corners. The main assumptions mi made codification 
of eave loading data as follows: 
4 Gust effects have been retained directly in the pressut pressure cosfficients; 
"consequently, only peak pressure coefficients are — in n the simplified d 
2. Eave height and cave width have not been considered as critical parameters -*- r 
‘the determination of effective pressure coefficients design purposes. 


tion (height and width) with the heaviest loading measured has been taken into a: 


b Terrain effects were significant ‘enough toe encourage the maintenance of 


the code model proposed has been based entirely on “the more conservative 
on Effective pressure coefficients have been determined as functions of the 
respective tributary area for which they have been measured, since this has 
been found to be a strong parameter for the load definition. Data for cave 
corners have not been considered since they are limited in quantity and a 


suggest extremely low values: in “comparison with other areas. However, if 
additional data support these 1 results, it is strongly recommended to include 
separate specifications for these areas in the chart. - 
5. For each configuration examined, the worst values of measured ne 
coefficients—regardless of wind direction—have been considered. To avoid the 
significant « conservatism involved in the design for such ‘‘worst “cases,” the 
designated coefficient has been established at 80% of the - maximum value. 


Arguments for the adoption of the 80th percentile le level have been pee oe 


Based on the assumptions stated, the codified eave pressure coefficients are 
peenante’ in Fig. 13. The effective pressure coefficients given in the chart are 
total coefficients appropriate to be used for the design « of eave sections s against 
: wind loads acting on both upper and lower surfaces. These values have oon 
- based on the bulk of experimental data for all configurations examined and 
do not reflect only th the ¢ representative ca cases presented in in this ran. The maximum > 


| 


uare foot) is given by the equation: _ 


0.0 00256 Vy 7, = mean hourly wind speed at eave height, 


_ miles per hour; and C, = the ‘designated effective pressure coefficient. Note 
that, by virtue of assumption 3 above, the designated c coefficients are appropriate — 


for open country terrain conditions, but they will give conservative wind | pressures 
when they apply to anotherexposure. 
“?_ It is recommended that the pressure . coefficients of Fig. 13 be wet | for wind 
loads on perimeter eaves of low-rise buildings whose eave height does not exceed | 
60 ft (18 m), (eave width/building width) < 0.12, and roof angle does not 
exceed 10° (8). The proposed height limit which is approximately double the 
maximum height of model tested, is justified from additional tests carried out ‘ 
- and not described in this paper. For higher heights or roof slopes (or both) 


FIG. 13.—Codified Pressure Coefficients for Maximum Loads on Seve Sections 


(Open-Cc Country Exposure) (Note: 1 1 ft? = 0. .093 m) 


been described in this paper. The most ‘significant ‘conslusions- can 


ow i. Mindente eave sections are very sensitive to wind loading. Suction loads 
_ dominate on the upper eave surface; pressure loads appear on the lower eave 
& surface ¥ with the exception of the overhang corner for which low positive or 
’ _ 2. Appropriate evaluation of eave wind loads cannot be made by considering _ 
the local pressures acting on various points of upper and lower eave surfaces — 

as used to be common practice in the past. Instead, the effective loading of | 
eave sections has has to be « determined nined by monitoring the statistics of the total 


> 
: es of building only) additional data are required for design specifications. 
| 
{| 
¢ 


3. A codified chart of design pressure coefficients has been suggested. ‘The ¢ 
design pressures are appropriate for perimeter eave sections of low-rise buildings 
not exceeding 60 ft (18 m) in height and 10° in roof angle. The eave width 1 
- should be smaller than = of the building width. Additional data are reed 
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= dynamic velocity pressure at eave height; ter 
V pat = mean wind speed at eave height; 


@ = roof angle; and 


rs Pun Ne 
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_ The optimization of nodal coordinates in discrete structures introduces several | 
of complexity above those of "Optimizing the element nt proportions v with 
the geometry held fixed. Major contributors to this complexity are: (1) Increased 7 
number of decision variables; (2) the different degree of nonlinearity in the 
- numerical behavior of these variables; and (3) the potential change in topology — 
_ under both fixed and variable geometry optimization. These difficulties prevented | 
effective utilization of available geometric optimization 


con far exceed those obtained assuming the geometry to be fixed. 
_ Various methods were developed by researchers to overcome the difficulties 

cited above. The majority of these methods are based on alternating optimization = 
& in two design rites the sizing variables subspace s and the coordinate variables | 


(13) treated s sizing and using a 
sequence of linear programming (SLP) approach. Most recently, Imai (4) and 
Schmit adapted the multiplier method (10,14), also treating all variables simulta- — 
neously, and overcame some of the aforementioned numerical difficulties. 
While the tools for numerical solution of complex involving 


optimization, from the designer’ s ‘point of view, still 1 remain. n. The ‘process is 
complex and costly, , even with the present decreasing costs of hardware and : 
central processing unit (CPU) time. One of the main reasons for the high cost 7 : 
is the numerous analyses which are usually rn in the optimization process. F 


improving intuitive, “experienced- -based For such purposes, ‘simplified 
processes, yielding approximate solutions, are much needed. This explains the — vd 
relative success of the stress-ratio method, otherwise known as the fully stressed _ 
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iterative design me the SLP, and the ‘optimality criterion techniques. ve ’ 
_ This paper proposes a two-phase sequence of processes addressing the above 
problems of structural geometry optimization. The method used is heuristic 
but simple and economically attractive. It is shown that, for a single-loading 
condition, satisfactory solutions can be obtained by fi first completing the optimiza- 
tion, considering only constant stress constraints and assuming the geometry 
to be fixed. This simple process, which usually produces a statically determinate 
truss, also simplifies the optimization of the nodal coordinates performed in 
the second phase. If the optimal design under multiple loads is treated as some 
combination of single-load designs (5,6,12,15,18,22), the final truss can be formed oi 
= and its ‘optimization refined, using more accurate methods. A similar process 
_ was suggested by Dobbs and Felton (1) in 1969. At that time, many of today’ * “a 
computer programs and methods were not to 
accuracy of approximate optimization We ot 


conceptual cast in the following 


minimize WiA, Q) = = Zpl L, 


in which which W = the weight of the truss; A. and L = the vectors of m member 


cross-sectional areas and lengths, respectively; Q = the vector of nodal coordi- 
mates; and p = the density of the material. ‘The set of constraint functions, 
g, are behavioral constraints, such as equilibrium, , compatibility and constitutive 
relations, | and side constraints restricting A and Q to vary between specifi ed 
upper and lower bounds. Ref. (17) is cited here as a representative of works 
which Qisassumedtobeconstamt,. 
_ As pointed out in the Introduction, it is advantageous in common practice 
not to optimize all variables simultancously. This is normally done by perf orming 


“assuming the § geometry, or Q, to be fixed; and the coordinate space, i in which 
Q is allowed to vary. The process in all these techniques is iterative and 
pr at alternating between the A and Q subspaces until convergence is 


_ be accompanied or followed aye ny proportioning. This necessity is due 

to changes in the equilibrium equations and member forces resulting from the 

- change in coordinates. Variation of node coordinates is usually based on a 
Search along a certain direction, € indicated by the particular gradient approach 7 


Soaaeees number of resizing and 1 re- e-analyses must be performed, requiring 
= or hundreds of analyses, which may be quite costly in a sizable problem. — 


Vy A representative method of geometric optimization is the technique used by ; 
Vanderpl 4). Il 
an aats (24). A for SADT calle sa 


| 
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Teached. € S€quence in which the Optimization 1s Carried Out in the two sudspaces 
“_——_ . It is important to note that while member sizes may be optimized under _~—_if 
| 


of analyses is s substantially reduced. However, a trial- and-e -error r process must Pi 
. performed to find the correct AQ whenever the standard values prove to 
g too large. Because of this problem and because of the internal linearization, 
“the program suffers from numerical instability and quite often stagnates. This 
(Spillers’) program, however, provides the user maximum control and flexibility. 
‘The stress ratio method » is used to update the member sizes, and it can be 
performed between geometric iterations in any number or as often as desired. -? 
Lipson and Gwin (8) also used the stress ratio method for optimization in 
‘the sizing subspace. Stiffness changes due to displacement constraints are 
performed using a similar scaling method. By using the complex method 3) 
in the geometry subspace, the gradient mates are completely avoided at the 


of the two- -phase geometric optimization were investigated. In the first sequence, 


- ‘Within the research program reported in this paper, two different a 


the geometry is varied with no attempts to optimize member sizes. In the second 

€ sequence, the optimization in the sizing subspace is first completed using = 

/: stress ratio method with fixed geometry and then the coordinates are — 

while adjusting the areas, A, in each iteration until convergence is reached. 
_Vanderplaats’ (23, 25) and Spillers and Kountouris’ (21) programs were used 


of techniques which consider all variables simultaneously, 7 in Ref. 4, for 


in n the ‘preceding section. . While the sequences investigated both lack the — 


they do offersome advantages. 
first sequence, which concentrates on the geometry, essential 


results were obtained. For this type of biel any general nonlinear programming 
methods and network theory are more appropriate. These methods were not 


in the case of multiple loadings. In the remaining sections, the second red Ms 


is proposed as a practical expeditious method of geometric optimization. An 
example is first presented for illustration and discussion. The conclusions drawn 
subsequently were based on several such examples. The results compared © 
satisfactorily with those obtained by more exact methods. = | 


general-purpose nonlinear optimization program CONMIN (23) whenever the iA 
= coordinates are modified; thereupon the member proportions are updated and a | 
7 the truss is reanalyzed. Size optimization is performed under fixed geometry 
: using the stress ratio method or a gradient method on the cross-section area 
arriation of this method is the process used by Spillers (19), where the 
_ geometry is varied in steps, AQ, obtained from a set of linearized Kuhn-Tucker — 
| 
this effort: on wever is sj d quite ing 
4 


was minimized ender the given load. Nodes | and 4 were constrained to move 
only in the vertical direction. The results of the optimization using Spillers’ 
and Vanderpiaats’ programs, _ discussed in the preceding sections, are shown 
in Figs. l(c) and 1(b), respectively. Bars which reached minimal cross sections, 
but did not totally vanish, are shown in dashed lines. In both cases the initial 
cross sections were all assumed to be equal, and the geometric optimization © 
began after a few stress ratio iterations without changing the initial geometry. | 
Since most of the members below the main diagonal vanished, the simpler 
truss shown in Fig. 1(d) was optimized; the results are shown in Figs. l(e) — 
and \(f), which correspond to Figs. 1(b) and 1(c), ‘respectively. 
Applying the proposed sequence of operations, the given truss in Fig. l(a) 
or 1(d) was first optimized using the stress ratio method and holding the oot 


> 


pay (62,29) ‘ 


FIG. 1—Example 1; E = 30 x 10° (psi); = 1 (Ibs /in.” ); allowable stress 


10°(psi); Loads in 10° (Ibs); numbers in indicate coordinate values 


‘fixed. This resulted in the truss shown in Fig. 1(g). This truss is fully stressed 
and will be shown later to be a ‘combination of two statically determinate 
configurations. Optimization of the geometry of the truss in Fig. 1(g) yielded 
the trusses in Figs. 1(h) and 1(i), which correspond to Figs. 1(b)-1(/), respectively. 
Note that in all cases the final locations of the loaded nodes are almost identical, 
while all other nodes are ‘‘floating.’’ Furthermore, if the kinks between nodes a 
d and 9 and nodes 4 and 9 were ae soso several other members would P. 
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GEOMETRIC OPTIMIZATION | 


- Optimizing the topology, or the connectivity of a Parone Structure, has been | a 
_ and still is a major stumbling block of structural optimization. While methods 
: for mixed-integer programming have been developed, the subject is still too 
complex for practical engineering applications. Two widely used methods of 
topological optimization are linear programming (2,11,20) and iterative stress a ; 
% ratio design (20). Using these methods, a set of redundants is eliminated from 
J a larger admissible set of bars connecting all admissible nodes. Most known 
_ geometric optimization programs, however, assume the topology of the structure 
to be fixed; some do not even allow for vanishing cross sections. Such programs = 
‘may be quite useful for optimizing a known topology. In preliminary design 
Ws the configuration problem (geometry and topology) is critical. In this section, - 
the role of statically determinate topologies is discussed as a potential, heuristic -_ 
basis for determining truss configurations. 
a Weight minimization Of stress- s-constrained trusses, under a single loading | 


toa ca feasible solution. If more than one such solution is optimal, a fully 
stressed indeterminate truss of equivalent weight may result (20). It follows 
that if, hypothetically, all determinate topologies are exhaustively a 
at least one of them will lead to the optimal topology and geometry 0 of the 
‘minimum weight 
e Figs. 2(a), 2(b), and 2(c) show ayes of the determinate topologies of the 
truss shown in Fig. 1(d). The results of geometric optimization of these trusses 
care shown in Figs. 2(d), 2(e), and af ), respectively. These configurations should 
‘be compared to those of Fig. 1 and to the theoretical optima in Figs. 2(g), 
2. and 2(i). It should be remarked that the fully stressed indeterminate solution, 7 
obtained by the stress ratio iterative method, shown in Fig. 1(g), is a combination z 
of and has equal weight to the two optimal d determinate trusses ; of Figs. 2(a) bf 
q - and 2(c). Furthermore, only one of these two trusses, the one shown in Fig. 
2(c), leads to the global optimum shown in Figs. 2(h) and 2(i). 
: ie Note in particular, that the solution obtained using Vanderplaats’ program, — 
- shown in Fig. 1(b), approaches the global solution in Fig. 2(h) while the solution re 
in Fig. ‘(c) obtained using Spillers’ program approaches the local optimum shown 
It is evident from the preceding example(s) that the initial iosidtingy’' and the ‘e 
method of optimization bear quite heavily on the optimum solution. It is also 
significant, however, that the coordinates of the loaded nodes converge in all _ 
cases to the true optimum, which suggests that this optimum is unique, oe 
the given loading and with constant stress constraints. Heuristicly, if this 
“uniqueness is tentatively assumed, the optimal coordinates of the loaded nodes — 


_ may be obtained by optimizing the geometry of any statically determinate set 
of members connecting these nodes to the supports. The topology may then > 
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=F,L, = 


sbog 


: _ in which F, = the force in member, i, is constant for determinate fixed-geometry 


configurations a and depends only on the loads, P (including the reactions), and 
the coordinates | points of application. Consider Michell’s 


FIG. 2.—Example 2; p = 1 (Ibs/in.*); allowable stress = = 22 x 10° (psi); loads in 
10° (ee): Numbers in “parentheses indicate coordinate values {in.); We weight 


“in which and = materiel constants; ads al 


‘For minimum V under fixed geometry, Michell’s showed that it is sufficient | 
to minimize Z since in this case C is constant. If the geometry is allowed | 
to vary, the contributions of both Z and C to the volume must be further 


minimized. Thus, if the weight reaches a minimum with the geometry fixed, 


under a single loading cont 
He showed that the quantity, C, in the 
om 
| 4 
i” (4) 
ruriher FeCUction May Oe ODlained DY Varying the potential Tunction, ©, 
ae by varying the coordinates. These notions are discussed in more detail in Ref. | 


"GEOMETRIC OPTIMIZATION 


‘Theoretically, only if W(A,Q) is convex and continuous, is there hope ‘for 


; a “reaching a global minimum, regardless of the optimization sequence. In problems 
with buckling constraints, the feasible design : space may be disjoint, as demon- 
ak. strated by Imai and Schmit (4), so that the general problem is nonconvex. | 
For many problems in practice, however, convexity can be assumed and 
7 _ advantageous sequences followed. The tremendous economy of fixed-geometry 
7 optimization, followed by geometric optimization of statically determinate trusses, 
“makes the proposed two-phase sequence very attractive, compared with the 
‘commonly used alternating “sequence. The number of analyses is ‘substantially 
-teduced. The solution may be subsequently modified using an exact method 
in a limited number of additional analyses including other constraints such as 
Euler buckling and allowable displacements. This, of course, assumes that the 
approximate and the exact solutions do not differ fundamentally, 
As shown | in Figs. 1(b) and l(c), the results of the iterative optimization — 
process, ‘continuously alternating between the coordinate and ; area design s sub- 
a “spaces may yield ambiguous results because of the topology. In large structures, 
this could be misleading. In optimizing determinate configurations, this deficiency 
is improved, and an experienced designer can extrapolate the final results and 
exercise judgement in providing ‘Tedundancy and buckling res resistance to 
-. } ‘It should be pointed out that, in this paper, only a single-loading condition — 
1 is considered. Multiple loadings, however, can be treated as combinations of 
single loadings (12,15,18,22). The stress-ratio method or programming procedures — 
- can, in such cases, be also effectively used in the fixed geometry phase. The 
problem of combining optimal geometries in the case of two-loading conditions — 
has been treated (6,12). a methodology f for case needs 


Concusions 
A two-phase sequence of optimization is proposed as a heuristic method for 


truss configurations in preliminary design. In the first phase, an 


optimal determinate configuration is obtained with the initial geometry held © 
‘fixed. In the second phase, the geometry is optimized. The advantages of this 
_ method lie in the substantial reduction in the number of analyses and the ined 
it offers the designer. It is also suggested that through the proposed technique 
te may be obtained a reasonable initial geometry that allows o one subsequently — s 
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C = defined in text, Eq. 
F = member force vector; 
= constraint functions; 
a = member subscript; Batoland (a « 
volume of truss material; 
weight: of truss material; hg © 
defined i in text, Eq.4; 
material constant; brs he 
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ol typieel antcumon em Pip. 7 
Cyciic SHEAR TRANSFER AND TYPE OF INTERFACE” 
; In the 1 1960’s Birkeland ant’ Birkeland (1) and Mast ( (4) developed a philosophy 
of connection design in which cracks are assumed to have occurred at disadvanta- 
-geous locations within the region of the connection. Reinforcement is then 
designed to transfer shear, normal force moment across these cracks \ 
the connection is loaded. No dependence is placed on the tensile strength of 
the concrete. They proposed that the reinforcement needed to transfer shear 
= across the cracks be designed using the “‘Shear-Friction’’ method of design. __ _ 
Provisions for the shear transfer reinforcement design u using the shear- ‘friction - 
method were subsequently included in the / American Concrete Institute, (ACI), — 7 
; Building Code (2). ‘iT hese provisions \ were based. on the test data obtained in 
- monotonic loading tests of specimens made from normal weight natural aggregates. 7 
Subsequent tests (7) showed that the shear transfer strength of lightweight concrete a 
under monotonic load i is inferior to that of normal weight concrete of the same 
compressive strength. ‘It was, , therefore, Bp er (7) that the ‘Shear transfer 7 


respectively, a as s 0. 75 and 0.85 times . the : shear teal ¢ strength of normal weight 
concrete of the same compressive strength and having the same reinforcement. — 


= S paper reports a study of the influence os cyclic shear transfer tahentier 
of aggregate type and of type of interface, i.e., a crack in monolithic concrete 
or an interface between concrete cast at different times. The primary purpose 
of the: study was to determine how sh shear transfer strength ‘equations based 7 
on monotonic loading tests should be modified when they are used in the desigz 


_ Note. —Discussion open until March I, 1982, To extend the closing date one month, 
_a written request must be filed with the Manager of Technical and Professional Publications, - a 
= Manuscript was submitted for review for possible publication on July 24, 1980. 5 
This paper is part of the Journal of the Structural Division, Proceedings of the American 
an of Civil Engineers, ©ASCE, Vol. 107, No. ST10, October, 1981. ISSN 0044- 
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— _ reversal of loads in earthquakes had not been verified experimentally, use has ae 
: ; been made of these provisions in seismic design. An example of this was the 7 


Scope. 
_ specimen of each pair was loaded monotonically to failure while the other was 
Subjected increasing, , cyclically reversing shear. 
— 


of Type of condition 
CNB composite | normal weight | bond broken uncracked 


Atom 
Shear Plane 
Closed stirrup pest 


‘Shear 


6 bar: és rat 


Fi A 


FIG. 1.—Typical Specimen (1 in. = 25.4 mm) 
In addition to the type of loading, the variables included in the tests were: 


_ Type of specimen, i.e., monolithic, or composite with the interface in the 
_ shear plane; (2) type of concrete, normal weight or lightweight; (3) condition - 
of the interface in the composite specimens, either bonded and initially cracked q 
or without bond and not subject to the initial cracking process before the shear 
test; and ( (4) the and strength of the ‘reinforcement thes shear 


q 
) 
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plane. The program of tests is set out out in Table 1 od 


Test Specimens.—Details of a typical specimen are shown in Fig. 1. It is — 
designed to be gripped by friction on faces “A” as ‘Shown i in The 
350i in.” (32,260 mm”) shear plane is subject to shear | without moment. ¥ Pato 


; wrapped around longitudinal reinforcement so as to ensure positive anchorage 

on both sides of the shear plane. Additional reinforcement was provided to 

prevent failure of the specimen away from the shear plane (Fig. 

The composite “specimens were cast in two stages, the interface between 

had two concretes lying i in the sheer Plane. The first cast part of these specimens 


hole rams 


FIG. 2.—Arrangements for Test 


s three old second part was cast. was 1s deliberately 


as 
roughened to an amplitude of 1/4 in. (6 mm) to conform to the requirements 
of the ACI Code (2). In Series CB dy apres effort was prc to obtain 


bond breaker of soft : soap and talc before the “second concrete | was Ss cast to 


prevent bond atthe interface, 
3 _ ‘The concrete mix proportions ; are are shown i in n Table 2. The normal weight concrete 


concrete. They were cured in the forms under polythene sheets, until the ll 
of test. The lightweight concrete specimens were moist- cured for seven ‘days 


| a 
q 
- 
| 
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TABLE 2. -—Concrete Mix Proportions be 


‘Fine aggregate | 1.81 2. 
Coarse aggregate 2. |. |. 
“Mix aged seven days at test; “Be four at 
“Type Ill portland 
a I portland cement. 


TABLE 3.—Properties of Monolithic Specimens 


yield forcement para-— 
vin pounds point (number meter pf,, 


f,, in kips in pounds 


Specimen square per square | per square 


number | inch stirrups) inch 


‘| 


three # ag 
| 6016 _three 3 


8 | ee #3 
ML3M cn 314 three#3 | 678 
ML3C | 4240 three #3 | 686 


“Average dry density 95 
Note: pei = 6.89 kPa; 


‘The normal weight concrete was made from a glacial outwash gravel of 5/8 : 
in. (16 mm) maximum size. The lightweight concrete was made from a ‘‘coated’”’ 


aggregate (Materialite) of f 1/2 in. in. (13 13 mm) maximum size. 


lized for studies of 
he concrete shall be 
) 
inkips 
MN2 wo#3 | 4722 | 389 
4 MN2 two # 466 38.6 
41.6 
57.1 
MLIM | 4140 | 66.0 | three#2 | 39% 25.8 
7 MLIC_ 66.0 three #2_ 396 a Bs 
(37.4 


_ The properties of the he monolithic and nate specimens are shown in Tables 
Arrangements.—The specimens were tested in the specially 
_ two-part frame shown in Fig. 2. Opposite sides of the specimen were attached 
to the two parts of the frame by gripping plates. The clamping force was provided — 
by post-tensioned high strength steel rods passing through the plates and through 
oversize holes in the ‘specimen. ° . Two soft rubber O- rings | were used to center 


‘TABLE 4. —Properties of Composite Specimens 


f,, in kips in 
per square two leg per square 
$tirrups 


a 


‘wor 
three # 3° 


three #3 by 
three #3 | 952 


be three #3 | 


Note: 1 psi = 6.89 kPa; 


- the ‘rods in the oversize holes to prevent the rods from bearing directly | on 
‘The shearing forces were provided by the diagonally opposed ‘Pairs of 60 

‘kip (266,880 N) capacity hydraulic center hole rams X and Y in Fig. 2. The 

loads were applied through load cells so that the applied shear could be monitored c 

continuously. Each linear assembly of rams, pull rods, and load cell was attached be 

to both parts of the frame by a spherical bearing at one end and a ‘ae 


q 
= 
Specimer 
number 
 CB2M | 6200 | 53.4 
a 6030 52.6 = 64 478 
| 6050 | 
| 
— 


at the other. This the | we vo parts ¢ of the to as 
Peers under load, and prevented the frame from restraining the deformation | 
of the specimen. Oil pressure was supplied to the rams through a four- -way . 
valve, by a Riehle ‘‘Pendomatic’’ pumping and measuring unit. oma bes eda 
-* Both the slip, (or relative motion parallel to the shear plane of the two halves »> 
~ the specimen), and separation, (or relative motion normal to the shear plane 
of the two halves of the specimen), were measured continuously using linear ‘ } 
differential transformers as the sensing elements of slip and separation gages 
to reference points embedded i in the face the specimen. 


—Testinc Procepures 
nitial Cracking.—All the specimens except t ose of Series were cracke 
In 1 Cracking.—All th hose of S CNB ked 
‘in the shear plane before being subject to shear loading. The crack was produced 
4 by appiying line loads to the back and front faces of the specimen along the 
tine of the ‘shear plane. The dilation of the Specimen normal to the shear pe 


reference frame surrounding the specimen. Loeding wa was continued until an 

- average dilation of 0.013 in. (0.33 mm) was obtained. When the line loads were =» 
removed, a residual dilation of close to 0.01 in. (0.25 mm) remained. This was > 
the average initial width of crack in the shear plane before shear was ‘as applied = 


to have occurred when the shear could not be increased further and both slip 

Those specimens having specimen numbers ending in C were subjected to 

_ the following program of cyclic loading. The shear was first continuously increased _ 

to 50% of the calculated shear strength V_(calc.) under monotonic loading, — 

after which it was reduced to zero. Shear of opposite sign was then applied \ 
to the specimen, again being increased to 50% of the calculated shear transfer _ 

; strength, before being reduced to zero to complete the first cycle of loading. © 

_ Each specimen was subject to 10 such cycles of loading. The maximum positive . 

and negative shears were then increased | by 8% of the calculated shear strength, 

i.e. to 58% of V _(calc.) for the next five loading cycles. After each succeeding 
five cycles of load the maximum shears were increased by the same amount. 

- This process was continued until failure of the specimen occurred. Failure was 

4 considered to have occurred when the shear could not be increased to the 


planned maximum shear for that cycle. 


In calculating V,(calc.), ACI strength reduction factor was taken as 1. 1.0, 


since e material strengths and specimen dimensions v were known accurately. a 
Ov (calc. ) _ based on the following equations previously proposed (6,7): __ 


(cale. = (0. + 400) 4. Ib (calc.) = (0. O.8pf, + 2.76)A,. ) 


Sel, lb [V_(calc.) = (0. f, + 1.38)A. 


| 
| 

| 
ig | 


but not more than 0.2 ft. A, nor 800 A. 4. 1b (. (5.5: 2 


in p= / Ay= area of shear transfer r reinforcement « crossing 


the shear plane, in.” (mm7); A. - area of shear plane, in.’ (mm_’); =! yield 
strength of shear transfer reinforcement, 


“Monotonic Loading Tests.— load aa 


a essentially similar for all specimens. Slip occurred along the preformed crack | 
- from the commencement of loading, at a progressively wl rate, as may 


cyclic 


ML2M , monotonic 


Lightweight 

Concrete , (Initial 

(Ini ially cracked) 


3.—Typical Shear-Slip Curves for tor Cyclic Loading the 
fai seen in een ad 3. The slip along the crack was accompanied by separation _ 
across the crack. In general, the separation was much smaller in lightweight % 
concrete specimens than in normal weight concrete specimens. __ specimens, 
P be tension cracks occurred adjacent to the shear plane in all specimens, 
‘more cracks occurring in the more heavily reinforced specimens 
_ specimens and in the bonded composite specimens, these cracks ' were ¢ uniformly 
distributed along the shear plane. In the composite specimens in which the 
bond at the interface was deliberately broken, the cracks tended to concentrate 
in the vicinity of the shear transfer reinforcement. In these specimens also, 
a few cracks perpendicular to the shear plane occurred near failure, at at the 


q 
ini 
_ 


of spalling adjacent to the shear apy in most spalling 
na was most pronounced in the specimens in which the bond at the interface was” 


broken. In these specimens the spalling o primarily in the vicinity 


changed 


Cyclic Loading Tests.—In these tests, ies response of the specimen 


as the number of cycles of and the level of 


‘Intermediate. 


cycle 


— iustrated in Figs. 4 and 5, which show typical shear-slip and slip-separation __ 
curves at different stages of loading, for both = ‘monolithic and composite _ 
_ The separation n values are the change i in semediies due to application of the La 
thane. To obtain the total separation, the initial crack width or the width of 
any gap at the interface due to use of the bond meen) o- must be added. a 
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0 CYCLIC TRANSF 

Response to the first — of loading i is characterized by a gradual reduction 


and by retention of almost all of the slip caused by the maximum shear until 
the shear reduces to about half its maximum value Fig. 4a). 
* Response to succeeding cycles of loading is characterized by a low ‘shear 
stiffness at low values and a gradual increase in shear stiffness with increase 
_ in shear in both positive and negative directions. As the number of load cycles =s J 
increases, the shear stiffness at low shears decreases and the increase in stiff: ness 
with increase in | shear becomes | greater. . This causes the -Shear-slip curve 
a complete loading cycle to assume a progressively more pinched appearance — 
as the number of loading cycles increases. This was particularly marked in 
the lightweight concrete specimens =e in the composite specimens without | 


cycles 
approaching — 


— 


FIG. 5. pre Nate of Typical Shear-Slip Relationships for Cyclic and Monotonic 
_ Loading; Initially Cracked Monolithic ithic Specimens | (1 | kip = = 4, 4,450 | N; Tin in. = 25.4 mm) © 


» 


bond at the interface, and was probably due to the smoother crack faces in : 
The slip at maximum shear increased slightly each cycle, for about the first — 
- cycles. Thereafter, for the normal weight monolithic and bonded composite __ 
_ specimens, the specimens responded in a stable manner Fig. (>). zat a 
_ For maximum shears of +0.5V,(calc.) the maximum slip was about +0. 007 a 
in. (0.18 mm) for the monolithic specimens and for the bonded composite — 
specimens, and was about +0. 015 in. . (0. 38 mm) for the composite specimens a 
_ The slip at maximum shear and the shape of the shear-slip curve remained 
essentially the same for a given maximum shear, until the maximum shear reached — by 
about 90% of the load at which failure occurred. At and above this load, the 
7 4 slip at maximum shear increased with each cycle, and by progressively increasing 
, amounts. The characteristic Shape of the shear-slip curve also changed in that 


as Lhe Shear’ iS increased in DOLD positive an negative cirections a 
| 
| 


after increasing as the shear increased, the shear stiffness then decreased again 
as the maximum shear was approached (Fig. 4). Failure occurred when the 
shear stiffness under increasing load reduced to zero, after) which slip in increased 5 


rapidly even though the load was decreasing. ot 


without bond at the interface, the slip at maximum shear increased by varying 


amounts in all load cycles. * ovitiegg de 


load due to shear acting in ‘one direction rather than in the other. This: was 
_ presumably due to the crack faces having an irregular profile. The general shape _— 


cycles of loading. _ The pend 
at ‘zero “shear usually stabilized at from 0.003-0.005 in. (0.08-0.13 mm) 
less than the initial value. After this had occurred, ‘the slip-separation curve 
‘stabilized, its form being similar to that shown in Fig. 5. The exact shape 
of the slip-separation curve would be a function of the actual profile of the 
crack faces in any particular case. As the maximum shear per cycle increases, 
by the slip and separation at maximum shear increased, but the separation 
at zero shear remained approximately constant until shortly before failure. 
"However, in the last few cycles before failure the separation at zero shear 
_ increased with each cycle of load, as shown in Fig. 5. This increase in separation 
at zero shear occurred in the same cycles in which the shear stiffness started 
_ to decrease when approaching maximum shear, and the slips: at | maximum shear 
increased significantly with each cycle load. 
In the case of the Series CNB specimens in . which bond at the interface 
_was broken, but the specimens were not cracked before testing, the separation — 
at zero shear did not decrease initially, but subsequent behavior was similar 
that of the specimens in the other test series. 
Diagonal tension cracks occurred in all specimens subjected to cyclic loading. 
Two sets of cracks occurred, both inclined at 45° to the shear plane, and at 
right angles to one another. The two sets of cracks opened and closed alternately, 
as the direction of the shear was reversed. Until close to failure, the cracks — 
were almost invisible at zero shear. They were fairly uniformly distributed along : 
o shear plane in the case of the monolithic specimens and the composite 
specimens in which the bond at the interface was not broken. In the case 
of those specimens in which the bond at the interface was broken, the diagonal _ 
— cracks tended to concentrate in the vicinity of the shear transfer _ 
r reinforcement. In these specimens a few cracks also occurred perpendicular 
the shear aligned with the shear transfer reinforcement. Gu 


Specimen failure was characterized by both slip and separation increasing oo 

with the load carried by the specimen either held constant or decreasing. =~ 
Failure Data.—In Tables 5 and 6 the following data are given for each ae 
_ tested: (1) Ultimate shear transfer "strength; (2) number of cycles of load to 


In the case of the lightweight concrete specimens and the composite specimens ~g 


_ The shear- slip curves were in general approximately symmetrical about both ie 


of the shear- -slip curve was unaltered, but the relative motion of the two halves - 


OI the specimen became Centered On a position different trom the Starting position. 
> 
q 
q 
4 
f the concrete : w 


Specimen "strength, cycles to failure, failure,® (monotonic 


(63 | — 0.031 


MLIC 


“Values of separation do not include initial crack width. ait 


at 


cycles to” ‘tallure, ilure,” 
failure i strength) 


0.035 0. 0260 


0.021 0.022 > 


_ 
TABLE 6.—Test Results for Composite Specimens 
— 
timate =| Number leparation | (Cyclic 
Series CB (Bonded, Initi ack od) 

| 

CNB3C 
"Values of separau © ck Wid 


q 


“which (3) ‘slip and separation at at ultimate | 
(4) ratio of ultimate strength of specimen subject to cyclic loading to that of 
companion specimen subject to monotonic loading, 

i. ‘The ultimate shear i is defined as the maximum shear carried by the specimen 
a during the test. In a a cyclic loading test this may be either the maximum shear — 
age! 

to which the specimen has been cycled, or the maximum | shear reached when = 
£ maximum shear was being increased at the end of a group of cycles of _ 


loading toaconstant maximum shear, 


a For cyclic loading, the values of slip and separation at failure listed in Table » 
‘5 are the averages of the values measured at the 
shear in he ast in the ant load cycle before failure. 
‘Test Resutrs ap: all 
ks To facilitate the Subsequent analyses of specimen slimes the ways in which 


= in the reinforcement crossing the crack, as a result of separation — 


of the rough crack forces when slip occurs, | 
2. By direct bearing of small asperities projecting from the faces of the crack. 

= By dowel action of the reinforcement crossing the crack; i.e. direct resistance 

of the bars to shearing action at the crack, 


_ Behavior of Monotonically Loaded Specimens. —At low values of applied shear, ; 


: of asperities on the faces of the crack. With little or no separation occurring, 
the reinforcement crossing the crack is not stretched significantly. Consequently 

_ little or no compression exists across the crack and hence no significant frictional _ 

~ resistance to shear can be developed. Also, with near zero slip between the 
a faces, shear resistance due to dowel action will be > very small. 1% tar veil 

As the applied shear is increased, increasing numbers of the asperities are — 

crushed and the intensity of bearing pressure on the remaining asperities increases 

_ ‘more rapidly than the applied shear. Deformations, therefore, increase more 
rapidly, i. e., the shear stiffness decreases; 
As slip and separation increase, an increasing amount of shear resistance 
will be provided by dowel ; action and by friction between the crack faces. — 
However, at shears approaching ultimate, the tensile stress in the reinforcement 
will approach its yield point, if secondary failure of the concrete near the i 
_ plane is prevented. The “‘dowel action’’ contribution to shear resistance will 
consequently reduce, and, at failure, it is believed that almost all the shear | 
Borge will derive from direct bearing of asperities on the crack faces ary ‘ 
- friction between the crack faces due to the large tensile force developed q 

- the reinforcement. The proportion of shear resisted by friction will increase 
as the total yield strength A,f, of the reinforcement increases. 

Dowel action probably contributes little to shear transfer strength when the 

yield of the shear transfer reinforcement is developed. This hypothesis 


qi 
| 
4 
| 
4 
shear is very Migh and litle OF No separalion OCCurs. indicates 
t at is stave mo e shear resistance am the 
4 
4 
4 
i 


to a shear plane, the shear transfer strength is increased by the same amount 
: 2 as would result from providing additional reinforcement with a yield strength — aa 
equal to the external compression force (5). If significant dowel action occurred _ 
_ after bar yielding one would expect a greater increase in shear strength | from 
a of the extra reinforcement than from the action of the ——_ 
s Because of the difference in concrete compressive strengths, the behavior 
- strengths of the normal-weight concrete and lightweight concrete monolithic 


wilt. 


Composite , rough interface, 


Pow no Weald, Initially cracked 
ods sc 
0.02 004 0.04 } 0.06 0.0 8 010 O12 O16 
FIG. 6. ._—Comparison of Typical Shear- Slip Relationships for Cyclic and 
Loading; Composite Specimens (1 kip = 4,450 N;1in.=25.4mm) 
- specimens cannot be directly compared. The behavior of the lightweight concrete 
_ specimens under monotonic loading was, however, generally similar to that 


be De 


a ‘loading the behavior and strength of the initially cracked composite —- 
with good bond at the interface (Series CB), was very similar to that of the 
initially cracked, monolithic normal-weight concrete specimens (Series MN). 

i is in with previous te tests 


682M, monotonic 
| 
| 
| 
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In cartier monotonic loading oni (8) of composite specimens in which wie 

at the interface was prevented, behavior and strength were found to be similar 

to that of comparable initially cracked, bonded composite specimens. In the 

_ ‘current study, the behavior of specimens CNB3M and CB3M confirmed this 4 
finding, but specimen CNB2M did not behave as well as specimen CB2M. 
The reason for this disparity is not known. 
- Behavior of | Cyclicallly Loaded Specimens.- —On first loading, behavior is as 
described for monotonically loaded specimens subject to low loads, with most 
_ of the shear resistance being due to direct bearing of asperities on the faces 

ala _ When the shear is reduced to zero, the slip does not reduce until the applied : 
shear is less than the net restoring force due to the combined effects of elastic 
deformation of the reinforcing bars and of the concrete, and any friction between 
the crack faces. When the applied shear is zero, some slip remains, corresponding _ 

to whatever frictional force continues to act between the crack faces. oe 


7 a _ As the shear is increased in the reverse direction for the first time, the 


mechanisms by which shear is resisted are the same as on first loading. ‘The — 
ing 
in load cycles i is different. In onder for 
the asperities to be brought into bearing, a slip must occur almost equal to 
the maximum slip that has previously occurred. Prior to this, shear resistance 
is developed only by dowel action of the reinforcement and by any small friction 
between the crack faces. The shear r stiffness i is consequently mt much lower initially 
than in in the first cycle of loading Fig. 4a). When the asperities come into 
bearing, further slip requires deformation of the asperities. The resistance of 4 
the asperities to deformation results in a sharp increase in resistance to —, 
“and an increasing tangent shear stiffness until the maximum shear is —_ es. 
_ On removing the shear, the response and underlying mechanism of behavior 7 
a P- similar to that in the first cycle. When the shear is reversed in a, 


the response and mechanism of behavior a are as described. 
_ With each cycle of load, the surfaces in contact are abraided and become 
- smoother, further reducing the frictional resistance to shear at low shears. This _ 
is reflected in a further reduction in the shear stiffness at low shears, which 
about 20% the effective stiffness at maximum shear in the same load 
_ As the m maximum num shear is increased, the m maximum slip an _ separations become ; 
Se This leads to the development of larger dowel forces and tensile stresses 


in the reinforcement. Consequently, an increasing fraction of the shear is — 
friction between the crack faces and by dowel action, 
_A further change in behavior occurs | in i the | last few loading cycles before 


— shear stiffness approaching maximum shear Starts to decrease, and the 

- and separation at maximum shear increase significantly in each cycle. This’ 

behavior is in large part due to local crushing of the crack faces. Some of © 

the mortar particles produced by this become trapped between the faces of 

_ crack, wedging it open, and also acting like ‘‘ball bearings’’ when the crack 

9 move relative to one another. Progressive yielding of the reinforcing da 
under the combination of direct tension, and of bending and shear resulting E 


- 
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It can be seen in Figs. 3(a) and 6(a) and Tables 5 and 6, that the behavior — 
and strength of the initially cracked composite specimens with good bond at - 

the roughened interface (Series CB) was very similar to that of the initially re 
cracked, monolithic, normal-weight concrete specimens (Series MN). (In Figs. 
3 and 6, the slip plotted for cyclic loading is the numerical average of the 
slip occurring at the maximum positive and negative shears in each loading — 
cycle.) In both cases, “until significant deterioration of the crack faces began 


was close. to | the ‘slip ovcursiog g at same shear in in the ¢ companion monotonic 
- loading test. In both Series MN and Series CB, the ratio of the strength under 
_ cyclic loading to the strength under monotonic loading averaged a little over 
_ _ Fig. 3(b) shows an example of the progressive deterioration of initially cracked, i 
monolithic, lightweight c concrete specimens subjected to cyclic load. At all levels: 
of load, the slip under cyclic loading was greater and the separation was less 
than that which occurred under monotonic loading at the same shear. It can 
s seen in Table 5 that the separations at failure in the lightweight concrete _ “9 
- specimens were very much less than in the comparable normal weight concrete : 
F specimens. This difference in behavior was noted in earlier (7) monotonic loading — 
tests. It was then postulated that this is because of the difference in minor — 
roughness of the crack faces in the two cases, due to the absence of sand — 
particles in the lightweight concrete. (It was pointed out that the magnitude 
_of the slip and separation at failure was of the same order as the size of S 

“| In the normal weight concrete the crack passes along the interface between 

_ the cement paste and the aggregate particles, because the bond strength between 

: the cement paste and the aggregate particles is less than the tensile strength — 

_ of the aggregate particles. The minor roughness of the crack faces, therefore, 
reflects the size and shape of the sand particles. In all lightweight concrete, 
the bond strength between ‘the cement t paste and the aggregate is greater than 
the tensile strength of the | aggregate. _ In this case the crack passes through — 

_ the aggregate particles and both ‘the minor and major roughness are less than _ 
_ Stable shear-slip behavior probably depends upon the interlocking of asperities 
on the crack faces. This would not be possible to the same extent in the case 
of the smoother crack faces of the lightweight concrete, as in the case of the 
‘Tougher crack faces of the normal weight concrete. The shear transfer strength _ 

_ of cracked, monolithic, lightweight concrete is less than that of cracked, monolithic 

- normal weight concrete under both monotonic and cyclic loading. As a result, | 
bs ratio of the shear transfer strength under cyclic loading to the shear = 
strength under monotonic loading is about 0.8 for both concretes, = © 
In the cyclic loading tests, the behavior of the unbonded composite specimens, 
deteriorated much more rapidly | than did that of the initially cracked, bonded 

7 ‘composite specimens. It can be seen in Fig. 6(b) that, in the unbonded composite 
specimens, the slip at maximum shear never stabilized for a given range of 
_ as was the case in the cracked, bonded composite specimens and ae 4 

‘monolithic, normal weight concrete » specimens. This This difference in behavior m 


g 
| 
f 
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the cracked, composite epecimens, which had good bond at the interface, 
it is likely that the crack propagated along and close to the formed interface. 
At cracking, sand particles would be exposed to create minor roughness conditions _ 

approaching those for a crack in monolithic concrete. This minor roughness 

in the cracked, bonded composite specimens was evidently able to resist the 
repeated abrasion of cyclic loading as effectively | as was the minor wahete 
in the cracked monolithic concrete, since in both cases the reduction in shear i 
transfer strength due to cyclic loading was about the same. 
a In the case of the composite specimens pape: bond, the reduction in shear _ 


- formed ee equivalent of a fine crack at the interface, and the major roughness 

was probably the same as for the cracked, bonded specimen. However, the 
- minor roughness would only be that formed incidentally in the mortar, when 
tee interface was a at the time the first half of the specimen was 


“consisting of exposed sand grains in the cracked, bonded composite specimens 


- occur perpendicular to the shear plane at the location of the reinforcement. 
This indicates: that for these specimens, larger dowel forces were developed 


~ composite specimens, where such crackirg did not occur. The greater dowel | 
forces probably resulted from the combination of lower axial tension in the 
_ bars and the occurrence of larger slips in the shear plane, , both are to oe 
q tn Fig. 7, average values of damping factor B for each test series are a 
against percentage of cycles to cause failure. The damping factor calculated 
Le 


Area under ‘‘skeleton curve”’ 
in which A, and A, are as shown in Fig. 7. No attempt was made to draw 
- skeleton curves for the shear-slip hysteresis loops obtained in these tests. Instead, 7 
_ it was assumed that the skeleton curve divides the hysteresis loop into two 


: be due to a difference in the minor roughness in the two types of composite _ “ 
j 
q 
: | | nly in the Case of the composite specimens without don id Splitting Cracks 
7 
. a bounded by LMNOPQL in Fig. 7. The damping factor — 
sé in this way has a maximum possible value of 0.318. The significance §& 
1 _ Of the individual numerical values is questionable, but it is thought that the 4 


1961 
ends i in to variation of the wehees are of interest and reflect the observed 


In all cases the damping factor dropped rapidly in the first few cycles of _ . 


loading to 50% of V_,(calc.), reflecting the greater amount of energy dissipated Ca 7 


_ in the initial shearing off of asperities, as compared to later cycles. Similarly, _ 
in the last few cycles before failure, the damping factor increased rapidly in 


Light-weight concrete « 


Rog > 0 


i cracked 


(b) Con weight cc concrete 
Percent of cycles to cause se failure 
Type 
all cases, due to the dissipation of energy as as the crack faces were increasingly — 
severely damaged approaching failure. me) Tomaso et 
_ In the interim, the damping factor decreased only a ‘small amount in the 
eases of the cracked, monolithic, normal weight concrete and the cracked, 
_ composite concrete. This setae to the the period during which the shear- slip 
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hysteresis stable for these ty types of specimen. after 
_ initial tendency to exhibit the same damping behavior, the damping fect 
increased significantly i in the cases of the cracked, monolithic, lightweight concrete 
‘and the composite concrete specimens without bond at the interface. This increase - 

‘in damping corresponds to the increasing slips at constant maximum shear, , 
noted in these specimens at relatively early stages in their loading. eS) er) : 
_ Ultimate Strength. —Comparing the measured shear transfer strengths listed 
in Tables 5 and 6 with the calculated strengths using Eqs. | and 2, which © > 
are listed in Tables 3 and 4, it can be seen that in all cases, Eqs. | and 2 — 
yielded a conservative estimate of ultimate shear transfer strength in in the monotonic = 
These equations have previously been shown (6,7,8) to predict shear mater 
‘strength quite closely, when the shear plane was horizontal at the time of —— 
_ the concrete. In these tests, the shear plane was vertical at the time of casting 
the concrete, and this may be the reason for the higher strengths obtained. 
— Itis probable that the circumstances leading to the higher than average values 
monotonic shear t transfer "Strength in these tests, would also lead to 


appropriate to base conclusions concerning the effect of cyclic loading on sane 
transfer strength, on the relative strengths of companion pairs of —— 


maximum shear was not t studied. However, because of the stable behavior of | 
the specimens until just below the failure load in the present tests, it is sa 
that the strength under randomly varying maximum shear would not differ 
significantly from the strength obtained in this study. ‘However, it should be 
noted that the shear stiffness will be much reduced at lower maximum loads 
per cycle than the maximum. previously reached, if the maximum shear per 
cycle varies randomly. 

‘The following conclusions ‘tons to shear transfer across an initially cracked 
shear plane or unbonded interface in composite concrete, when secondary failure . 
of nearby concrete due to dowel splitting effects or diagonal tension cracks 
is prevented by the provision of suitable reinforcement. __ 


It is fecommended that ‘shear transfer Strength under cyclically reversing 


for the following ca cases: a a shear plane i in n normal weight o or r lightweight r monolithic — t 
- concrete and a rough interface between concretes cast at different times, when 
good bond has been obtained at the interface. 


a 2. If the interface is roughened as specified in Section 11. 1.9 9 of ACI 318- 7 7 
_ (2) and if good bond is obtained at the interface, then after cracking shear 
_ transfer behavior under both monotonic and cyclic loading will be essentially 
- the same as in the case of shear transfer across a crack in monolithic concrete. _ 


_ 3. If bond at the interface between concretes cast at different times is sis destroyed, 7 


= 


and the calculated strength using Eqs. | and 2. | 
| Only one loading history was used in this study, involving a progressivel —_ 
| 
4 
| 
=. 
> | 
| 
a 


shear transfer behavior under cyclic loading deteriorates rapidly and the shea 
transfer strength is only. about 0.6 of the shear transfer strength under monotonic 


strength, of interface roughness and bond between concretes cast at different 
ite times; every precaution should be taken to obtain adequate roughness and good 
a... bond at interfaces that ¢ are likely to be subjected to cyclically reversing load, 
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The following symbols are used in this 

= total area of shear transfer reinforcement esenting & the shear plane; 

area under hysteresis loop ‘‘skeleton curve’ 

= area bounded by aes in Fig. 


= 


= compressive strength concrete on 12 in. 
sof, yield strength of reinforcement; 

) = calculated shear transfer strength; in sone 
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Unrate Ca CAPACITY OF CRCULAR 
By Joseph. A. Yura,’ M. ASCE, Nicholas Zettlemoyer,’ 


cross section provides for buoyancy needed during installation in the ocean 
environment. Tubular members are also used in many truss type structures | 
which require long slender compression members, since the tubular cross section | 
exhibits a high strength- to-weight ratio. 
_ In most instances, the connections of tubular members are Pome by | on 
penetration welding of the carefully contoured ends of the branch members 
to the continuous chord of the truss. Typical simple tubular joints are shown _ 
in Fig. 1 along with some of the geometric terms used in defining joint behavior. z 
The word “‘simple’’ implies no stiffening devices and, in the case of K joints 


an The first widely used design approach for s simple tubular j which 

in 1973 in the American Petroleum Institute’s Recommended Practice (API RP | 

_ 2A), was based on a punching shear stress concept (2). The branch axial load — 
and moments were converted to an average applied shear stress in the chord | 
wall which attempts to punch out an area of chord around the perimeter of | 

4 the branch-chord intersection. The punching shear resistance was — 
_ from a lower bound fit to the results from 38 ultimate strength experiments 
available at that time (11). The experimental ultimate capacity was used because 

Of the significant post- yield test strength and the inadequacy of theoretical 


atic -plastic solutions. A safety fi factor of 1.8 was incorporated in the ee) 


When more experimental evidence became available, the single shear 
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: _ Circular tubular members are used extensively in offshore structures because = ff 
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_ resistance fi as was no longer a lower bound, so the API concept was was expanded 
_in 1977 to include adjustments for type-of-loading and type-of-joint. The = 
son of the more complex capacity equations with test data showed rome 
. improvement in reliability, but scatter was still significant (12). ae 
The objective of this paper is to establish improved ultimate € capacity 


development benefits from a carefully screened data base, "which “existed 
’ of early 1979. Some of the data are new, unpublished results from the oe 

of Texas, including data for out-of-plane moment loading on the branch (25). 
_ After preliminary discussion of the data base and the equation format, iis 
_ paper : systematically addresses axially loaded T, Y, , and DT joints, , axially loaded 
K joints, and moment loaded joints. In each section there is a review of i 
ormance characteristics, a of the new and an 


«FIG. 1 —Typical Simple Tubular Joints: 
of each equation’s accuracy relative to the data base. The prediction accuracy © 
_of the equations in the 1973 and 1977 editions of the American Petroleum Institute’s 
(API) Recommended Practice (API RP 2A) (2), with the factor of safety of 
1.8 removed, is also _— for the same data base. While other f published capacity 


exist (7, ant space limitations have excluded them from 


_ Data Base Validity. —Comparisons of joint capacity test results, with the 
- predictions of existing capacity equations, generally show extensive scatter 
(11,12). There are a number of likely sources of the scatter. First, many of 
"the tests were conducted on very small joints—the chord diameters were less 
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CIRCULAR TUBULAR JOINTS 


and local behavior, particularly where fracture is involved. Second, for many — = 
. the material properties were not measured and only minimum specifie ed 
“yield strengths | are given. Yet, it is common for the actual yield strength in 
steel pipe to be more than 20% greater than the specified value. Third, much mt 
of the test documentation lacks deformation data. If deformations are large, 7 
®@ capacity may be controlled by deflection rather than load. Fourth, occasionally 
resulted from member rather joint in some cases x 


8k An evaluation of the test data was undertaken to provide a sound experimental — ae 

basis for tubular joint capacity equations. Only those joints with chord diameters o 
larger than 3. in. (140 mm) and measured material Properties were considered 


selected does permit a significant ‘number of data to be retained. Each test 
_ report was examined for deformation data and a capacity limit other than joint | 
failure, such as general cross section yielding of the member(s). These and 
any duplicate data were eliminated from the base. Also, a few ccmqueedin 
T results were discounted ‘since short chord lengths gave al abn oe 
to the boundary conditions. 
— In the screened data base, ‘the minimum load from three potential failure 
criteria has been used to define the useful capacity of each joint. One criterion 
is the maximum load attained during the test. This is the limit recorded in 
_ virtually all test reports. However, in the axial tension situation, the load associated _ 
by with first crack i is considered a more appropriate limit. In laboratory specimens 
the branch is usually overdesigned to ‘eliminate branch failure. While this 
- overdesign is believed to have little or no effect on the first crack load, it 
_ results in unrealistic reserve strength beyond the crack’s first appearance eed : 
‘The third failure criterion incorporated in the data base is the load at an | 
~ excessive deformation limit. For axial loading, it is assumed that a branch member z 
reaches its practical deformation limit when the strain along | its entire length 
is four times the yield strain (1,24). It is further assumed that the comparable — 
practical limit of local joint deformation at each branch end is twice the — 
yield deformation. Thus, the deformation I limit, A, may be 


in which L = branch len th; and E = Young’ s modulus. A member length 

- of 30 times its diameter (a typical upper limit for offshore structures) is used 


to evaluate the deformation limit. However, ‘for joint and load types 


plot for the joint is relatively flat c over a a broad r range of Soleeatien: a 
Equation Format.—lIt is difficult for one formula to accurately predict the 
ultimate capacity of tubular joints for a large range of joint geometries and 
loading configurations. Among other things, the post-yield behavior varies. By 2s 
introdpcing more than one equation it is possible to improve both the simplicity 


and | the accuracy. y. The authors have found that separation of the data into — 


«Six groups is advantageous and does not result in an unmanageable set of equations. > 


so 
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API RP 2A makes use of the ‘punching shear concept to predict capacity — 


g 11,12). While the concept has proven useful, it does not properly reflect 
the actual mode of failure at a joint. Several recent equation recommendations 
have avoided the failure mode issue by describing the capacity in terms of _ 


. ultimate branch force or moment (7,16,20). The branch loading format has = 
: selected for the new equation development. Use of branch force and monet 


should simplify the design procedure, since member end forces and moments 

are the common products of a structural analysis, 

au Performance Characteristics. —Typical load-deflection behavior of T ‘and DT “ 

tubular joints subjected to axial compression and axial tension is shown in 


BRANCH ‘ 
(4st CRACK 


Man 


Fig. 2 ie a ‘companion set t of et of four tests -_ matied member sizes aie 
of the tension specimens is characterized by a increasing 


> 
— ’ 
1 
«FAG. 2.—Load-Deflection Behavior of Axially Loaded T and DT Joints (21) 
> 
"distortion of the cross section (ovaling) until a first crack has formed. [No rd . 
1 | 7 first crack data are actually reported (20); the crack loads shown are for illustration _ 
— 


only.) The branch load contin continues i increase gross separation occurs. => 
The DT joint behavior in tension is similar to the T joint, but both the stiffness 
Compression loading produces the joint capacity. Failure is usually 
* associated with yielding, buckling, and gross distortion of the chord wall. Again, Bi 
the DT joint is weaker than the T joint, except possibly at B values near 1.0. 


The initial stiffness of the joint is similar for both compression and tension 
loading. For purposes, , the theoretical elastic Tesponse of the branch 


4 The joint ‘stiffaess is significantly less than the member ‘stiffness, | especially 
_ Compression Data and Capacity Equations. _The published compression data 
- for T, Y, and DT joints, which are considered reliable, are presented in Table — 
1. The ultimate test load, Pr, is nondimensionalized by the factor T” F,, which» 
7 is related to the plastic bending strength | per unit length of | plate. ‘This factor 
7 is common to many suggested ultimate strength formulations for circular ~— 

_ rectangular tubular connections, especially those based on a yield-line analysis — 
(7,13, 19). Unless otherwise noted, the branches intersect the chord at a 90° — 

The T and > 4 data are plotted in Fig. 3. Even though the range of y (defined 

in Fig. 1) is quite large, the nondimensionalized ultimate load shows good 4 
correlation to B alone. The Y tests show an average 18% increase in strength | 
above the companion T tests. This increase is close to the 21% indicated by 

_ the K, term in the API formulas (2) in which K, = (1 + sin 6)/2 sin lee 

he z and Y data are available for B up to 0. 7. At larger values of B, some =) 

= were conducted (10), but the short chord lengths affected the results O. 

Ss actual structures, T and Y joint failures at the B = 1.0 condition are unlikely © 

because the chord member will fail due to the high bending moments introduced 
the branch force. Thus, the following could ‘i used to represent 


P= 


§ the limited number and scope of Y data. Also, 50, tension te: tests to be  ahiiaee 
dater suggest that the K, adjustment is nonconservative. 
It is interesting to reflect on the constant 3.4 in Eq. 2. As B approaches 
€ the branch represents more of a point load on wel local chord surface, 
which might be viewed as a flat plate in the relative sense. A yield-line poor ral 
of a square plate with a perpendicular point load gives a constant of 4.0, 
independent of the plate planar dimension. Thus, the limiting value of 3.4 seems 
reasonable with reference to the yield-line approach. 
The DT data from Table | are plotted in Fig. 4 and cover the full practical 
‘Tange of y and B. At low values of B, the ultimate load is similar to that 
developed in T joints. However, in the mid range of 8, DT joints show a — 
capacity. At high B, the increases significantly because the 
load i is transferred between the branches, primarily by membrane action rather 
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TABLE 1. _—Data T, Y, and DT Joints | in Compression 


or. 
CB-70-0.2" 
-10-0.4° 
CB-100-0.2" 


0.721 | 0.721 | 1.011 


0.832 | 0.832 
0.802 | 0.802 
| 46.5 . 0.663 | 0.663 


CS-40-0.2° 
CS-40-0. 4° 
CS-40-0.6" 
CS-40-1.0° 
CS-70-0.2" 
 CS-70-0.4° 
CS-70-0.6" 
CS-70-1.0° 
 CS-100-0.2"* 
 CS-100-0.4° 
CG-1° 


4 
= 
in | in | in 7 | Col | Cok | 8/ 
mile | mile | | in | | umn | umn | Equa 
te [mega | | | sin 8/ | 8/ | tion 
Test | ters ters ters cals | y 6 ifs, | bat. 
— | (3) | 4 | | | | | ao | ay | 
139.8 | 101.6 | 6.5 | 323 | 10.8] 0.727 | 17.26 | 0.884] 0.884 | 1.003 
' 139.8 | 101.6 | 6.5 | 323 | 10.8] 0.727 | 18.70 | 0.957 | 0.957 | 1.086 
if 164.5 | 42.7 | 4.7 | 440 | 17.5 | 0.258 | 8.65 | 1.121 | 1.121 | 1.037 7 
a : 164.5 | 76.3 | 4.7 | 440 | 17.5 | 0.462 | 11.97 | 0.870 | 0.870 | 0.980 
319.5 | 60.5 35.5 | 0.190 | 7.05 | 1.017 | 1.017 | 1.008 wa 
ae 319.5 | 139.8 | 45 | 410 | 35.5 | 0.440 | 12.23 | 0.761 | 0.761 | 1.043 
| 455.7 | 89.1 49 390 46.5 | 0.195 | 6.67 | 0.856 | 0.856 | 0.936 
a 455.7 | 165.2 | 4.9 | 390 | 46.5 | 0.362 | 10.41 
| 1645 | 763 
CB’-70-0.4" | 319.5 | 139.8 
 CB’-100-0.4" | 455.7 | 165.2 
165.2 | 42.7 | 4.7 | 480 | 17.5 | 0.260 | 6.84 | 0.894 | 0.894 | 1.013 
_— 165.2 | 76.3 | 4.7 | 480 | 17.5 | 0.460 | 9.61 | 0.701 | 0.949 | 1.022 
; 165.2 | 114.3 | 4.7 | 480 | 17.5 | 0.690 | 12.93 | 0.617 | 0.881 | 1.020 
| 165.2 | 165.2 | 4.7 480 | 17.5 | 1.000 | 30.49 | 0.572 | 0.817 | 1.033 
318.5 | 60.5 36.2 | 0.190 | 6.01 | 0.858 | 0.858 | 1.023 
) | 318.5 | 139.8 420 36.2 | 0.435 | 10.01 | 0.618 | 0.820 | 1.099 
318.5 | 165.2 35.4 | 0.520 | 11.29 | 0.594 | 0.833 | 1.113 
: 318.5 | 3185 | 4.5 | 410 | 35.4] 1.000 | 37.04 | 0.563 | 0.804 | 1.255 
457.2 | 89.1 390 | 47.6 | 0.195 | 5.54 | 0.709 | 0.709 | 0.933 — 
457.2 | 165.2 | 4.8 390 47.6 | 0.361 | 8.90 | 0.616 | 0.740 | 1.100 a. 
139.8 | 139.8 | 323 | 10.3] 1.000 | 32.37 | 0.712 | 1.017 | 1.097) 
CG-2* 139.8 | 114.3 | 65 | 323° | 10.3} 0.818 | 17.99 | 0.755 | 1.079 | 1.113 
CG-3* | 139.8 | 101.6 | 6.5 | 323 | 10.3] 0.727 | 14.46 | 0.752 | 1.039 | 1.075 
139.8 | 76.3 | 65 | 323 | 10.3 | 0.546 | 11.37 | 0.823 | 1.167 1.083 
139.8 | 486 | 65 | 323 | 10.3| 0.348 | 8.42] 0.956] 1.126| 1.0 
| 190.1 | 48.3 | 4.69 310 | 20.3 | 0.254 | 6.94 | 0.881 | 0.881 | 1.036 
2° ree 193.7 | 48.3 | 6.50 | 330 | 14.9] 0.250 | 6.81 | 0.964 | 0.964 | 1.024 
193.7 | 48.3 | 9.39 280 | 11.0 | 0.250 | 6.97 | 1.081 | 1.081 | 1.048 
188.9 | 101.6 | 4.65 | 310 | 20.3 | 0.538 | 11.28 | 0.676 | 0.955 | 1.085 
193.7 | 101.6 | 6.50 | 330 | 149] 0.520 | 10.78 | 0.734 | 1.030 | 1.061 
193.7 | 101.6 | 9.30 | 280 | 10.4| 0.520 | 11.37 | 0.861 | 1.210 
» 7 | 490.0 | 159.0 | 4.56 | 310 | 20.8| 0.837 | 19.20 | 0.620 | 0.885 | 1.164 _ 
193.7 | 159.0 | 6.50} 333 | 14.9] 0.820 | 16.00 | 0.598 | 0.854 | 0.985 
4 9¢ Pee 193.7 | 159.0 | 9.35 | 280 | 10.4 | 0.820 | 18.23 | 0.759 | 1.084 1122 
340 | 15.9 | 0.522 | 10.97 | 0.729 | 1.025 | 1.077 
: || | | 340 | 15.9 | 0.522 | 10.63 | 0.707 0.992 | 1.044 
| 


ee ~ Note: For column 9, the mean of ‘column 9is 0. 779, the standard deviation | 0.144, and the e coefficient 
of variation 18.5%. For column 10, the mean of column 10 is 0.923, the standard deviation 0.136, q 


- and the coefficient of variation is 14.7%. For column 11, the mean of column 11 is | _ - standard — 
0.076, and the coefficient of variationis7.1%. 


“Deformation controlled. 


_*From Ref. 18. pa 
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4 suggests that no significant effect on joint ade. saint 
_ is clarified in Fig. 5 where the data from Table | a are organized as a function 
of y. At B less than 0.6, a variation of y from 10-36 shows no effect. At 
higher B, there is evidence of aslighteffect. 


| lower bound equation: in Fig. 4 is 


q 


in which Q, = = 1.0 for B ‘< 0.6; and Q, = = 0.3/B(1 — 58/6) for Bp > 6. The 
sin @ term is used in Eq. 3 even though no X (or DY) joints are reported - 
in Table 1. Use of the DT equation for X joints is inferred from the relationship 


of Tand Y data 


= 


FIG. 5. of on or Strength 


Data and Capacity data which provide information 
on the formation of the first crack are limited. Table 2 Gocuments the 19 tests | 
which were found to satisfy the reliable data conditions. The 16 tests from _ 
the large JSSC program (20) are are not included in the table because the first 
crack loads were not reported. The test capacity, P,, is the first crack load © - 
in all cases, although = many instances the load at the deformation limit was | 
_ Fig. 6 plots the test load versus B. At best the data can be described as 
scattered. Various attempts were made to organize the data in a tighter a arrange-_ 
ment using other test variables, such as +, but little improvement was observed. | 
Actually, the scatter cannot be reduced significantly because of the variation 7 of 
in the first crack } oad itself. The four data between ‘points | a and b represent 


fix 


= 
10 


CIRCULAR TUBULAR JOINTS 


os surface cracks, or hints thickness cracks, 2 are reported. Also, fracture charac- 


323.9] 60.3 | 6. 5 | 0. | 0.855 | 0.890 
323.9 | 101.6 5 | 0. 0.668 | 0.695 
323.9 | 273.1 4 ; 0.403 | 0.420 


0.748 | 0.748 


variation 39. 8%. For colema mean 411, the standard deviation 0. 603, and the 


a *From Ref. 6. bas, ’ 


— 


A lower bound of the data follows the same trend as the compression T 


data (Fig. 3). Tests with the largest chord diameters are particularly close to 
the lower bound. The limited DT and Y data also follow the lower bound. 
A comparison of Beale’ s companion T and Y data G) indicates th the additional 


1973 
almost identical specimens, yet the first crack load varied by 60%. When failures 
_ 3 are due to fracture of the material, considerable scatter should be expected — 
Test ters | ters | ters | cals |_ ‘ai 73 77 
(7) | @) (9) (10) | (11) 
Tit 323.9 | 73.0 | 12.7 | 283 | 12.3 | 0.218 | 6.82 | 1.173 | 0.838 | 0.906 
| 323.9 | 73.0 | | 25.5 | 0.214 | 17.56 | 2.471 | 1.765 2.352 
7 - 406.4 | 88.9 | 64 | 283 | 31.5 | 0.209 | 15.61 | 2.111 | 1.508 218 3 
323.9 | 141.3 64 | 25.5 | 0.425 | 14.83 | 1.051 | 0.751 q 
= 219.1 | 141.3 16.8 | 0.643 | 24.20 | 1.278 | 0.917 
as 323.9 | 141.3 | 6.4 | 283 | 25.5 | 0.425 | 16.78 | 1.189 | 0.849 1.462, ; 
| 323.9 | 141.3 | 64 | 283 | 25.5 | 0.425 | 15.22 | 1.078 | 0.770 
| 219.1 | 141.3 | 15.5 | 0.645 | 39.4 | 2.124 | 1.526 | 2.517 
| 141.3 15.5 | 0.645 | 39.1 | 2.108 | 1.514 | 2.497 
219.1 | 141.3 «15.5 | 0.645 | 33.6 | 1.811 | 1.301 | 2.146 
{ 323.9 | 60.3 64 290 | 25.5 | 0.173 | 5.71 | 0.994 | 0.710 0.854 
323.9 | 101.6 64 (290 | 25.5 | 0.302} 9.14 | 0.991 | 0.651 | 1.002 
323.9 | 273.1 | 6.4 | 290 | 25.5 | 0.840 | 19.05 | 0.573 | 0.488 0.984 
087 
ose 
| Note: For column 9, the mean is 1.252, the standard deviation 0.593, and the coefficient of variation - : 
a 
— — 


= LOADATIstCRACK 


02 
6.—Tension of T, and DT Jo 


_ tests on matched 90° and 30° DT specimens (15) also suggest the API K, term 
overestimates the increase in joint capacity. Novikov’s tests are not includ ed 


in Table 2 because no first crack loads or material Properties were — 


experimental of the specimens are to 
Eqs. 2 and 3 and to API predictions in Table | . Histograms of the results 
are given (26). For these tests, the 1973 API equation typically overestimates. a 
_ the true capacity. The adjustments made in the 1977 API edition still result __ 
in significant scatter, but the ratios of test/theory are closer to the ideal value 
of 1.0. It can be seen that Eqs. 2 and a provide a a good evaluation of the 
lower bound joint capacity with reduced scatter. 
_ The statistical summary given in Table 2 for the tension data shows that 
all methods of strength prediction result in considerable scatter, while the 
: coefficients of variation are 40% or greater. Once again, the 1977 API approach 


provides peep close to the data mean, while Eq. 2 represents essentially > 


"Performance Characteristics.— —Various studies on K joints have indicated that 


nad 
bea) 
| 
FORFOUR 
mes 4 | eal = or 
| 


CK-100-0.2 


TABLE 3.—Data Summary: Axially Loaded K Joints 
« 
te | | | mega- sin | 9/ | 9/ |Equa- 
Pg pas- | API | API-| tion 
Test ters | ters ters | cals y | B 73 | 77) | 
| | | | ® [ao | ay] a2 | 
 Bouwkamp, = 90°; 0,=45°(Ref. 4,5) 
323.9 | 168.3 | 6.4] 290 | 25.0 | 0.520 | 0.236 | 36.90 |2.150 | 2.039 | 1.725 
168.3 | 60.3 | 5.6] 395 | 14.6 | 0.358 | 0.186 | 20.60 |2.049 1.756] 1.223 
5.6| 395 | 19.2 | 0.406 | 0.024 | 20.60 |1.664 | 1.299] 1.041 
4.8} 425 | 28.1 | 0.326 | 0.324 | 14.58 |1.308 |1.201] 0.986 
7.1] 290 | 11.3 | 0.358 | 0.186 | 18.22 |1.957 | 1.677] 1.082 
| 12.7| 494 | 20.0 | 0.431 | 0.905 | 16.70 | 1.040 | 1.040 a _ 
BF 508.0 1219.0 | 12.7] 494 | 20.0 | 0.431 | 0.905 | 16.70 |1.040 | 1.040 | 1.339 
508.0 | 324.0 | 12.7 20.0 | 0.638 | 0.154 | 25.92 | 1.086 | 0.988 0.996 
10 508.0 | 324.0 12.7] 494 | 20.0 | 0.638 | 0.154 | 27.31 }1.145 | 1.041 
508.0 | 324.0 | 12.7] 276 | 20.0 | 0.638 | 0.154 | 34.79 |1.458 | 1.326 | 1.337 
12 $08.0 | 324.0 |12.7} 276 | 20.0 | 0.638 | 0.154 | 33.54 |1.406 | 1.279 | 1.289 
165.5 | 42.7 | 4.6 | 17.8 | 0.258 | 0.730 | 9.34 | 1.214 | 1.214] 0.925 
318.4] 60.5 | 4.4] 412 | 36.2 | 0.190 | 1.425 | 10.82 | 1.544 | 1.544 | 1.544 
| 456.9 | 89.1 | 4.9] 402 | 46.5 | 0.195 | 1.357] 9.14 |1.178 | 1.178 | 1.286 
458.9] 165.2 | 4.9] 402 | 46.8 | 0.360 | 0.182 | 24.38 |1.700 }1.454) 1.439 
TK-400.2 165.5 | 42.7 | 4.7] 490 | 17.6 | 0.258 | 2.460 | 9.28 |1.003 | 1.003] 1.118 
TK-40-0.4 165.9 | 76.3 | 4.7] 490 | 17.6 | 0.460 | 0.760 | 15.68 |0.951 | 0.951 | 1.084 
TK-40-0.6 | 165.2] 114.3 | 490 | 17.6 | 0.692 | 0.031 | 31.05 | 1.223 | 0.995 | 1.057 
TK-70-0.2 | 318.4] 60.5 | 4.5] 422 | 35.4 | 0.190] 3.849] 8.93 | 1.063 | 1.063 | 1.274 
STK-70-0.4 321.4] 139.8 | 4.5] 422 35.7 | 0.435 | 0.855 | 16.24 | 0.842 | 0.842 | 1.247 
 TK-70-0.6 317.7 165.2 | 4.5) 422 | 35.3 | 0.520 | 0.509 | 21.11 J0.919] 
456.9 | 89.1] 451 | 45.7 | 0.195 | 3.714 | 7.38 [0.793 | 0.793 | 1.039 
q *7K-100-0.4 | 458.9 | 165.2 | 5.0| 432 | 45.1 | 0.360] 1.364 | 13.44 |0.785 | 0.785 | 1.313 
165.2 | 76.3 | 1.6] 348 | 51.6 | 0.461 | 0.131 | 38.20 | 1.673 | 1.419 | 1.853, 
165.2 | 76.3 | 2.3] 288 | 35.9 | 0.461 | 0.131 | 28.05 | 1.370 | 1.162 | 1.361 
| 65.2 | 60:5 | 23] 338 | 35.9 | 0.366 | 0.165 | 23.52 | 1.447 | 1.227] 1.362 
& 165.2 | 48.6 | 2.3] 340 | 35.9 | 0.294 | 0.205 | 23.00 |1.762 1.494] 1.565 
065.2 | 76.3 | 2.3] 294 | 35.9 | 0.461 | 0.131 | 35.56 | 1.737 1.473] 1.725 
165.2 | 76.3 | 3.2] 314 | 24.8 | 0.461 | 0.131 | 30.04 |1.639 | 1.390] 1.457 
165.2] 76.3 | 4.5] 348 | 18.4 | 0.461 | 0.131 | 28.33 | 1.691 1.434) 1.374 
165.2 | 76.3 | 278 | 13.8 | 0.461 | 0.131 | 24.78 | 1.612 | 1.367] 1.202 
165.2 | 76.3 | 6.0] 282 | 13.8 | 0.461 | 0.131 | 22.17 | 1.443 | 1.223 107% 
: . Bigs 165.2} 76.3 | 1.6] 358 | 51.6 | 0.461 | 0.131 | 42.92 |2.270 | 1.925 | 2.082 _ 
- Reage f 165.2 | 76.3 | 2.3] 339 | 35.9 | 0.461 | 0.131 | 41.80 |2.470 | 2.095 | 2.033 
don 165.2} 60.5 | 342 | 35.9 | 0.366 | 0.165 | 31.31 }2.325 1.972] 1.813 
165.2 | 48.6 | 2.3] 270 | 35.9 | 0.294 | 0.205 | 32.72 |3.025 |2.565| 2.226 
165.2] 76.3 | 2.3] 285 | 35.9 | 0.461 | 0.131 | 35.34 | 2.084 | 1.767 | 1.715 
165.2] 76.3 | 3.2] 263 | 25.8 | 0.461 | 0.131 | 39.97 |2.602 | 2.207 1939 


323.9 | 60.3 | 6.5 15.4 | 0.173 | 1.225 oats 1.853 | 1.853 | 1. 
323.9 | 101.6 | 6.5 15.4 | 0.302 | 0.237 | 13.91°] 1.614 | 1.395 | 0.945 


= - Note: For column 10, the mean is 1.521, the standard deviation 0.497, and the coefficient of variation 
& 7%. For column 11, the mean is 1.355, the standard deviation 0.393, and the coefficient of variation 


in tension branch. «4 


the gap, g, between the intersecting branches has Py influence on 
- the capacity, as do the variables which affect single branch joints (16,23). If — 
_ the gap is large relative to the members, the joint performs like a single branch 
joint. As the gap approaches zero, the overall joint strength is increased because 
the bending stiffness of the chord wall between the branches is increased. For 
_ large gaps, the joint capacity is usually controlled by plastic bending and buckling - 
of the chord wall in the vacinity of the compression branch. For small ati iy 
the joint strength depends on conditions i in the gap. 

Many axial K tests have had « one 90° branch loaded in compression, a am 7 
a 45° branch loaded in tension. In these cases the gap was in tension, and ,- 
failure was governed in part by fracture of the chord or welds in this region. 


_ Therefore, considerable scatter in the results could have been expected. a 


_ Nakajima experimental program included companion tests which had the direction 
branch loads reversed (14). Generally, the value of ‘sin 6 was slightly 
lowe er when the 45° branch was in compression. In this case, compressive ve chord — 

stresses: developed in the gap ape. More importantly, the need for a K, 
‘termasin API wasnot supported. 

‘a Data and | Capacity Equation.—A summary of the experimental K data i is shown 

in Table 3. Except for two tests (25), the branches in each specimen ; are the 

- same size. 8. refers: to the angle of the compression branch | and 6, to the 

tension branch. The joint failure loads listed in Table 3 are those on the 
compression branch (maximum load or load at the deformation limit), except 
for the TK joints (three branches). In these tests the single branch in tension © 
failed by developing a crack. we 
The effect of the gap ratio, g/d, on joint strength is shown in Fig. 7. The © 

P,, are ‘nondimensionalized by Eq. 2, the capacity prediction 


1976 
16 165.2 | 76.3 | 6.0] 360 | 13.8 | 0.461 | 0.131 | 20.81 | 1.635 1.386] 1.010 
165.2 | 76.3 | 6.0] 292 | 13.8 | 0.461 | 0.131 | 22.56 |1.772 | 1.503] 1.095 
$07.2 | 326.4 352 22.8 | 0.643 | 0.117 | 28.03 | 1.364 | 1.181 L052 
507.2 | 326.4 22.8 | 0.643 | 0.117 | 29.67 | 1.437 | 1.250 
Zimmermann, = 60°;8,= 60°(Ref.27) 
419.0 | 168.3 | 10.0] 340 | 21.0 | 0.402 | 0.284 | 15.32 |1.129 | 1.054 
419.0 | 168.3 | 10.0] 340 | 21.0 | 0.402 | 0.284 | 15.42 | 1.136 | 1.061 | 0.888 
419.0 | 168.3 | 10.0] 340 | 21.0 | 0.402 | 0.284 | 15.02 |1.107| 1.033] 0.865 
{ 419.0 | 168.3 | 10.0] 239 | 21.0 | 0.402 | 0.284 | 17.28 | 1.273 | 1.189 0995 
| 
> | 
| 


= and Y joints. Other researchers ha’ have related the joint capacity to g/D ‘rather - 


than g/d, but an analysis of both approaches indicates that g/d provides a __ 
a better organization of the test data. For gaps greater than the branch diameter, ¥ 


"NAKAJIMA 


ALL OTHERS 


Fi iG. 7. —Effect of of Gap on Axially K 


as a T or Y Eg. 2. For gaps the ‘bound capacity 


on he ities oe of the various capacity equations for K joints is “ 
"given in Table 3. It is seen that Eq. 4 gives a good prediction of the lower _ 
bound capacity, while the 1973 and 1977 API equations show greater scatter. i 


Since the of data at g/ d< < 0. 2 tends to ‘the 
e 4t is suggested that the 17 Nakajima tests be neglected when anaes the 


accuracy of the various equations, 


| 

P, = G.4 + 198) | 
| 

plane 


“bending (IPB) out- bending (OPB) moment is in Fig. 8 (20). 
Typically, the joint subjected to in-plane bending reaches a maximum load at 
moderate deformation, with actual failure due to plastic bending and sae 
7 of the chord wall on the compressive side of the branch, and fracture through © 
the chord wall on the tension side. For out- -of-plane bending, the chord — 
distorts locally, resulting in lower joint stiffness and strength. Fracture may 
occur, but only after excessive deformation. The practical deformation limit | 
is defined as the joint rotation corresponding to the angle at the end of a simply 
supported beam, which is uniformly loaded to four times the first-yield - 
A beam length of 30 times the branch diameter is used to establish the limit. 
Out-of-plane results are apt to be more sensitive to 6 than in-plane cei, . 
_ At high B values the branch transfers load to the chord primarily through membrane 
action in the chord, as to wel Improved capacity and stiffness 
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ong FIG. 8.—Load- Deflection Behavior of Tubular Joints Subjected to Bending a 

are expected, just as for the DT compression specimens. On the other hand, 

_ the mode of load transfer for in-plane bending load is not altered significantly 


ip Plane Bending Data and Capacity Equation.—In tests used to study the 7 


bending capacity, the branch is normally loaded at its end as a cantilever. Since 
the moment varies along the length of the branch, the useful strength must 
_ be related to a specific location, such h as the maximum moment which | occurs — 
at the chord centerline, M., or the moment at the chord surface, M,, where 
the failure actually occurs. In subsequent tables both moments are listed for 
: each test cited, but M, is used to develop the capacity formula since it provides 
a slightly better organization of the data. Capacity estimates based on M, are a 
conservative for situations in which as analysis 
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A2-X-90"° 24 

- Note: For column 10, the mean is 3.212, the standard deviation 0.577, and the coefficient of variation 
18%. For column 11, the mean is 1.627, the standard deviation 0.185, and the coefficient of variation : 


Ref. 25. ka ahe 


TEST 
TEST 


in- Bending Loading of T ond K Joints 


TABLE 4.—Data Summary: In-Plane Bending 
in in in qn Col- | Col- | umn 
q li- li- | mega- sin 8/ 8/ |Equa 
me- | me- | me-| | 6/ | 6/ | APL | API-| tion 
ters | ters | ters cals dT’F, aT’F, 77 | 
219.1 | 71.6] 6.3] 314 | 17.4 | 0.327] 9.24] 9.64 | 3.818 | 1.909 200 
298.5 101.6 | 7.2] 294 | 20.7 | 0.340] 9.24] 9.79 | 3.485 | 1.743 | 1.171 
; 7 7* | 219.1. | 101.6 | 5.5 | 305 | 19.9 | 0.464 | 12.45 | 12.99 | 3.482 | 1.741 
219.1 101.6 | 8.4 13.0 | 0.464 | 9.81 | 10.24 | 3.118 | 1.559] 1.006 
219.1 | 101.6 | 10.0] 368 | 11.0 | 0.464} 9.34] 9.75 | 3.121 | 1.560} 0.956 
219.1 | 139.7 | 6.0] 314 | 18.3 | 0.638 | 16.34 | 17.06 | 3.396 | 1.698 1316 
12" 219.1 | 139.7 | 8.8] 422 | 12.4 | 0.638 | 12.90} 13.46 | 3.013] 1.507} 1.039 
| 298.5 | 193.7 | 7.3] 296 | 20.4 | 0.649 | 17.48 | 18.53 | 3.448 | 1.724 1408 &g 
219.1] 177.8 | 5.9] 314 | 18.6 | 0.812 | 20.85 | 21.77 | 2.978 | 1.489 | 1.384 
219.1] 177.8 | 8.6] 422 | 12.7 | 0.812 | 17.68 | 18.46 | 2.832 | 1.416 1174 
— B-70-0.2° 318.5 | 60.5| 4.4] 441 | 362] 0.190] 6.45] 8.02] 3682/1841] 1.150 | 
318.5 | 139.8] 4.4] 441 | 36.2 | 0.439 | 12.48 | 14.81 | 3.083 | 1.542 1329 
| 457.2] 89.1] 4.8] 402 | 47.6 | 0.195 | 7.37] 9.52 | 3.776 | 1.888) 1.297 
—-B-100-0.4" 457.2 | 165.2] 4.8] 402 | 47.6 | 0.361 | 11.79] 14.53 | 3.263 1.631] 1.437) 
= 
— 


ecrosen 


5.—Data Out- of- Plane Bending 


| 457.2 
4 


~ 


$07.2 | 455.9 
$07.2 | 326.4 
507.2 | 455.9 


al Note: All tests are deformation controlled except BL-100-0.4. For column 10, the mean is 1. 470, 

_ the standard deviation 0.353, and the coefficient of variation 24%. For column 11, the mean is 1. “i 
the standard deviation 0.302, and the coefficient of variation 19.4%. “ead column 12, the y 
1.171, the standard deviation 0.179, and the coefficient of variation 15. 3%. ont 


(344760, 


| mil- | mil- | mil- in Ms, | M. | 8/ 
li- li- | mega- sin sin 8/ 8/ Equa- 
| me- | me-| pas- | | APL | API | tion 
ters | ters | ters | cals B |d7°F,| aT°F,| 73 77 
(1) (2) | 3) | 4 | (6) | (7) | (8) | (9) | (10) | (11) | (12) 
— ; 
165.2 | 76.3 | 4.5] 471 | 18.4 | 0.462 | 5.46] 6.08 | 1.570 | 1.570 
q BL-70-0.2* | 318.5 | 60.5 | 4.4 36.2 | 0.190 | 4.10] 5.10 | 2.340 | 2.340 | 1.084 
7  -BL-70-0.4" 318.5 | 139.8] 4.4] 441 | 36.2 | 0.439 6.56 | 1.366 | 1.366 | 1.068 
 BL-100-0.2"_ | 457.2 | 89.1 | 4.8) 402 | 47.6 0.195 | 4.12] 5.32} 2.101 | 2.111 
165.2 | 4.8] 402 | 47.6 | 0.361 4.36 5.37 | 1.207 | 1.207] 0.9200 
171.5 |11.1| 352 | 22.8] 0.338 | 4.39| 5.31 | 1.618 | 1.618 | 0.952 4 
7 171.5 | 11.1] 352 | 22.8 | 0.338 | 4.83] 5.84 | 1.780 | 1.780 1.047 | 
326.4} 11.1] 352 | 22.8 | 0.644 } 7.29) 8.81 1.403 | 1.403 | 1.147 
326.4 | 11.1] 352 | 22.8 | 0.644 | 7.86] 9.50 | 1.513 | 1.513 | 1.236 
455.9 | 11.1 252 22.8 | 0.899 | 13.40] 14.79 | 1.398 | 1.398 ol : 
| 455.9 | 11.1] 352 | 22.8 | 0.899 | 13.45] 14.86 | 1.403 | 1.403 | 1.305 
- 4 507.2 | 455.9 te 22.8 | 0.899 | 14.96 | 18.70 | 1.200 | 1.595 | 1.452 
— 507.2 | 455.9 352 | 22.8 | 0.899 | 14.90] 18.62 | 1.195 | 1.588 | 1.446 
352 | 22.8 | 0.644] 5.83] 7.04 | 1.122] 1.122 0.917 
352 22.8 | 0.899 | 12.31] 15.39 | 0.987] 1.312] 1.195 
7 352 | 22.8 | 0.644 | 8.32] 10.06 | 1.602 | 1.602 | 1.309 , 
352 | 22.8 | 0.899 | 14.67] 18.33 | 1.177] 1.564 
| 
) § 
| 


loads on the two branches of the K joint are listed individually and were applied | . 
Ses: that the branch moments opposed each other. The useful moment capacity 
at the chord surface is plotted in Fig. 9 as a function of 8. There is significant 
scatter, but the lower bound follows the same trend as the axial compocusion a 


Eq. 5 implies that the in- -plane n moment can be converted to an equivalent axial 

load by dividing the moment by 0.8d, which is close to 0.785d, the distance 
between the resultant tension and compressive forces on a tubular cross section 
with an elastic bending stress distribution. Gibstein suggested that the ultimate — 


‘moment is related to ott in addition to B (8). However, the JSSC data do 


4 


- Out-of-Plane Bending Data and | Capacity Equation. —Test data ain 17 branches 

_ Subjected to out-of-plane bending loads are summarized in Table 5 « The data 

cover a wide range of y and £ values, , and include tests on T, Y, , and K joint i 


The Yura data provide tests which established a test 


In the K joint tests the branch shear loads were in the same direction. The 7 
- maximum | load in each branch at failure is recorded in Table 5 and shown 
4 in Fig. 10. In test C2-1, both branches failed simultaneously, with P,, equal - 
to three times Boe . This test produced the lowest failure loads, as shown by 7 ; 
the open data points. However, in test C2-2 the branches did not fail simulta-_ 
I; neously, with P,, equal to one and one-half times P,,. Each branch showed 
an improved load capacity over the C2-1 
___ Despite the aforementioned K joint differences, it appears one equation can 
reasonably predict the lower bound joint capacity. The plot in Fig. 10 shows 


“J ‘The comparisons of test and theory listed in Tables 4 and 5 show that the qi 
API equations predict a lower bound to the data but are very conservative. 
; Eqs. 5 and 6 also predict a lower bound, but with reduced scatter and ee 
A reliable data base of tests tubular joints has been established and 
summarized. This base excludes data without measured material properties, a 
wae 
properties, deformation information, or fracture information. The base 
- 


a 
| 
q 
NOL Uibdstein trend. so € simple lower oun piven > 1S 
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‘also tests ats with « chord diameters. less than 5.5 in. (140 since scaling 
questions exist for smaller sizes. val soe 4 lobo ow? adr 
- _ The screened data base has been used to develop improved lower bound - 
_ capacity formulas for a variety of joint and load configurations. These new _ 
equations « are summarized in Fig. 11. One formula is recommended for both — 
tension and compression in the branch of a K joint. This formula is identical 
to that for T and Y joints in tension and compression, except for the additional 
gap factor. Branch moment capacities can be converted | to equivalent axial 
load capacities by dividing by 80% of the branch 


"experimental results. The new equ equations ‘predict ‘the lower bound 
with less scatter than the API equations. The 1977 API equation sororities ; 
-Toughly the mean bel the T axial data, but is lower bound for joints with applied — : 


for the new equations, , the bad fi API PI equation shows distanct improvement over | 
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BENDING 


OUT-OF -PLANE 
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FIG. —Ultim 
There are a number of questions which ‘should bee addressed prior to design 
implementation of the new equations. First, one or more safety factors should — 
be established. With t the | new w equations, one common safety factor | becomes: 


this case lower bound) with a relatively constant degree of scatter. Second, 

loadi interaction should be considered. Some sort of interaction equation is 

to predict joint capacity under combined branch loading. Also, the effect of 
& load in the chord should be addressed (11,22). Until reliable test data are generated > 


4 American Institute of Steel Construction, Specification for the Design, Fabrication 
Erection of Structural Steal for Buildings, 7th ed., 1969. 
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_ VIBRATIONS OF ‘CABLE-REINFORCED 


By O. G. Vinogradov,’ D. J. Malcolm,’ A. M. ASCE eldest) 
Pure, ron 57 and P. G. Glockner,’ M. ASCE ad? it 


_ Due to the high flexibility of an air-supported structure, its equilibrium shape E 


and associated external loading can be found only by considering the interaction © 


: between the structure and the loads, e.g., aeroelasticity in the case of wind 


~ Joad . As a step toward solving this complex interaction problem, the dynamic — 


characteristics of air- supported structures under fixed external loads should be — 


determined; progress towards this first goal has been limited (7), but such an 

- analysis will itself be useful in design by showing the sensitivity of dynamic — 
characteristics to different static loads and internal pressure. 
_ The free and forced vibrations of an initially stressed membrane shell of — 
revolution have been studied by Leonard (3) who found that the magnitude 
of the initial stress has a significant effect on the natural frequencies. A wind 

gust analysis of a cable reinforcing net for an air- r-supported « dome is , presented 
(6) in which the wind load has been considered as a static part plus small 
random fluctuations, and the spatial distribution of the wind pressure has a 
taken in accordance with data from tests on rigid domes (4). Ref. 61 is an interesting 
attempt to simulate the complex loading of the flexible structure. a al biome 396508 

‘The present paper deals with the dynamic characteristics of cable-reinforced _ 
inflatable caps and the influence of internal pre pressure and static loads on ree — 
characteristics. The results presented are based on a numerical method, with | 

_ corresponding computer code, developed for the nonlinear static and small forced — 

oscillation analysis of arbitrary cable-reinforced air-supported structures. 

+ The inflatable cable-reinforced structure is modeled by a cable network with 
the action - of the membrane neglected, except for transmitting the pressure — 
load to the cables. The pressure load is considered as distributed along cable — 4 
elements while all other loads, including inertia forces in the dynamic problem, 7 
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- concentrated at nodes. It is assumed that cab cables con b can be divided into segments divided into segments — 
such a way that each | is loaded by a ‘uniform pressure line- load lying in 


of elements sufficiently to approximate any pressure loading. The equilibrium — 
of each segment under the uniform plane pressure line-load is described by | 
simple nonlinear algebraic equation. Another advantage i is that the finite element 
matrix 
The cables are considered a as s elastic elements. ‘undergoing large 
It is assumed that the change in shape of the structure does not influence | 
the magnitude of the external pressure load. Under these assumptions the 


geometric nonlinear static problem is reduced to the determination of the 
equilibrium: shape of the structure under the given load. . In the Program, 


an additional ‘condition of “virtual shapes” i is introduced. 
_ When the equilibrium shape of the structure, and the corresponding g stiffness 
matrix are known, the small, forced oscillations are investigated by transforming a 
_ the static stiffness matrix into the dynamic complex stiffness matrix, including 
inertia terms from concentrated masses. are! 
| ‘Then method of assembling | the overall structural stiffness matrix at each cycle 
is based on that used in the SAP IV computer program (1), in which the a. 
of non-zero elements in the overall structure stiffness matrix is noted and is 
used in the solution of the linear equations, 
_ The numerical analysis has been applied to a a spherical membrane cap y comprised 
_ of two meridional reinforcing cables lying in mutually orthogonal | vertical planes: 
and having 3 degrees of freedom, and a spherical cap comprising 14 reinforcing 
cables lying along geodesic lines and having 147 degrees of freedom. Strong 
interdependence has been found between internal pressure, static load, and 
_ Some of the results presented here w were ediscussed at the 1980 Spring Convention — 7 
‘The cable segment shown in Fig. 1 is restrained against displacement at ends 
q and 2, and is in equilibrium under the plane uniform pressure, p, and end 


4 
| 
> 


orces H, = H,=H 


‘ which L, = unstretched length, and EA = axial stiffness of t of the cable. . 
follows from: the geometry ‘of Fig. 


‘in which L = the span of the cable, the following equation may be obtained oa 
‘or the unknown 1 angle, asban ile 

8+ecos@=—-—€ .. 
wea 

in which e = L,/L and E = 0.5 p L,/EA. The stiffness of the element in _ 
the deformed configuration is found a considering changes in end forces due 


to small end and results its (Appendix I 


4 


q 
dost 


_ The solution of the nonlinear static problem i is found by the sietation method. 
a Using the linearized structure stiffness matrix and residual nodal load forces 
as external forces, the coordinates for each node at step j + | are — ert al 
in which 1 u’, = the coordinate increments or displacements; se aie eet 


tic elongation, AL,, of a segment is 
| 
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The computation is continued until the residual force vector insignifi- 

cant. The process may fail to converge if the difference between the initial” 

and final geometry is too large at any one step. If this occurs, then a value , 

of ’ = a < 1 should be chosen in Eq. 8 so that displacements will be | small 

at this step. Applying this method of relaxation coefficient adjustment to cables — 

subjected to pressure loading is equivalent to requiring convexity of the cable _ 

segment at each step, i.e., e’ > 1. If this condition is not satisfied, \’ must 
be reduced and the step repeated. This procedure ensures convergence of oe 
iterative process for any initial geometry with unstretched cables. 


vue ‘od 


« equations describing static equilibrium, and by : adding inertia forces 

in which 8 = = y/m: = logarithmic decrement of the material; m = mass; 
os = circular frequency; and v = amplitude of vibration. For quasi-elastic materials © 


in a relatively narrow frequency range, the real and imaginary terms of the | 
y w frequency rang ginary 


_ For quasi- nani materials with _— 8, , and after separation of real and imaginary 


- The real part of 6* is, ‘therefore, expressed by Eq. 4 where ‘the value « of 6 

is known from the static analysis, and the imaginary part of 6* is given | through - 
Eq. 14. The substitution @ = 6* transforms the static stiffness matrix, K (Eq. 7 - 

5), into the complex stiffness matrix, K*, which is transformed into the dynamic 
complex stiffness matrix by adding the inertia stiffness of the concentrated = 
masses, —mw’, to corresponding diagonal elements. External dynamic forces _ 


are also considered as complex in order to describe the phase se tse) 


1 
) 
rum shape and 
assume all masses concentrated atthe nodes. 
a ‘The algebraic system describing small vibrations is obtained, according to 
} 
| 
7 
1g 


_ The FORTRAN computer program consists of two parts: (1) The solution a 

: of the nonlinear static problem; and (2) the determination of amplitudes and ¥ 

phases of small vibrations. The switch to the second part occurs when th “a 

The effectiveness of the program depends mainly on the « 
its static part because it contains two kinds of iterative processes: (1) The 


ee of Eq. 4 for the cable segment; and (2) finding the equilibrium shape 


FIG. 2.—Two-Cable Spherical (EA 6 KN = 4. 6 > 
4 is solved by the method with an initial value of 


a which i is on a half-sine-wave approximation and ensures ‘convergence 


within 5% of the correct shape after2-4 cycles, 
SS initial geometry of the structure was taken as that of the structure, under — 
small internal uniform pressure, without concentrated loads. To ensure conver. 


were e needed fora value of 8, = = 0. 25 x 107°. = 
4g Execution times on the University of Colne s Honeywell DPS level 2 , 
Computer System were | sec/iteration for the structure with 60 cable segments 
Fig. 7). Calculations of and of vibeations for same structure 


| 
q 8 
| 
: 


numerically. Results for the first ‘structure | were with 
exact results, while some of the results for the second structure were compared 
_ with corresponding data for the static case (5), 


Z Spherical Membrane Cap with Two Reinforcing Cables.—If the change i in shape 


are proportional to the square root of the pressure 


Because the segment and, thus, the overall stiffness matrix is directly Proportional — 


to the internal pressure, it follows that the eigenvalues of this matrix are also _ 
proportional to the internal pressure and, therefore, natural 
are proportional to the square root of this variable. ; 


VERTICAL 
FREQUENCIES 


HORIZONTAL 


1 
LOGARITHM OF DIMENSIONLESS PRESSURE 


3. —Log Natural Versus Log | Internal Pressure for 


Two-Cable Spherical C a 

_ The results presented in Fig. 3 confirm this relationship for the two-cable 
~ reinforced spherical cap over a wide range of pressures. However, for small 
values of p, when the displacements under the static load are large, and also g 

7 for large values of the pressure when elastic elongations are large, again resulting — 

‘in large displacements, the deviation of the real horizontal and vertical natural 
frequencies (solid lines), from those described by Eq. 17 (broken lines), becomes 
significant. The zones of deviation coincide with those where the deflections 


from the initially ‘Spherical Shape are (Fig. 3). ‘There is, therefore, a 


“proportional to the square root of a magnitude of internal pressure. The width 
of this range depends on the stiffness of the cables and decreases with decreasing — 

Fig. 4 shows the effect, to scale, of vertical concentrated static loads, such 
as snow and water, on the change in shape of the structure. The nse ll 7 
between the change of shape and change of resonant frequencies can be read - 

from ail 4, using points on the stiffness curve as reference points. As esi 
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be seen from the diagram, the relationship between force and natural frequencies _ 
_ is nonlinear. For relatively small changes in shape (u/R =< 0.15), it may be 
Z approximated by a straight line so that a decrease in natural frequencies becomes __ 
proportional to the displacements or the Yo Pine 
a Fig. 5 shows results for a horizontally applied static load. These results n 
demonstrate a reversed influence, in stiffness and 
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4" DISPLACEMENT RATIO, 
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VERTICAL 
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Spherical (pR = N, 7.55 Ibf) 


FORCE RATIO, F/pR 
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FIG. 5.—Vertical and Horizontal Natural Frequencies Versus Shape for Two- aes 

‘Under real loading conditions, the spectrum of wind load comprises low 

_ frequencies, which may be considered as an equivalent static pressure, as well 
a: as high frequencies. Under the action of slowly applied static forces, a structure» 
undergoes changes of shape and corresponding changes in its dynamic response» 

y _ which may differ significantly from the response in the undeformed configuration. — 

It is important, therefore, to check for any overlap between the spectrum of 

loading frequencies and that of Ge: structure in the displaced configuration. 


The case of wind pressure, static load acting 
vs shown in Fig. 6 (to scale). Note that the static concentrated load increases 
the effect of the static wind pressure on the dynamic response of the structure, ies ; ; 
because the stiffness of the structure depends on its shape and pressure. Under 
wind loading, both of these parameters are affected. The existence of two peaks 
% in the diagrams of Fig. 6 is a result of interaction between vertical and horizontal © a 
_ vibrations at corresponding resonant frequencies due to the unsymmetrical ha” 
Spherical Membrane Cap with 14 Reinforcing Cables. —It was found from 
“a previous static analysis (5) that the ‘‘best’’ cable configuration consists of cables _ : 
that follow lines of geodesic curvature. The cable- reinforced ae ia 


Pint #0.04 psi,284 Pa 


as 


« 
= 
2 
Qa 
= 
da 


f 


7 T he static load has been taken as the dead weight of cables applied at nodes. 


4 
Tet 
a FIG. 6.—Vibrations of Two-Cable Spherical Cap Loaded by Wind Pressure — cd 
rs to obtain the spectrum of vertical and horizontal vibrations. A unit harmonic _ 
— force is applied at the apex of the cap (Fig. 7) in the vertical and horizontal 
1 _ Fig. 8 shows the relationship between internal pressure and the first three i 


‘The ‘results. of the analysis (solid lines) do ‘not 
7 exactly to those predicted by Eq. 17 (dotted lines), the reason being the same 
as for the two-cable reinforced membrane. The smaller the internal pressure, a 
the larger the change in shape due to the static load. As the internal pressure — 
increases and the importance of dead weight diminishes, the the results approach ~ 
_ those given by Eq. 17. Deviation from those results will also occur if the clastic 
extension of ey is sufficient to change the undeformed oo 
We ELEVATION 


FIG. 7. —Fourteen-c Cable Reinforcement Spherical = 13.3 x kN = 

_& An example of the dynamic response of a spherical cap ) under an internal - 
pressure of 0.4 in. (96 Pa) of water, is presented in Fig. 9. It can be seen 
_ that natural frequencies in the vertical direction are close to each other od 
are distributed in groups. The resonances would not be so sharp if acoustic 
damping and corresponding added masses had been included. As a result of - 
relatively small material damping = 0. 01), “the interaction between vertical 
and horizontal vibrations at resonance have become significant. Satie: 

_ In order to estimate the influence of the static part of the wind load on 
the dynamic characteristics of the cap, an idealized wind pressure of +3.8 
‘Db/sq ft (182 Pa) plus corresponding to the leeward side, minus to the windward | 
side) was as applied to the structure with an an internal (of 4. 3 ft 
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FIG. 9.—Magnification Factor Versus Frequency for 14- Cable Spherical Cap Under ; 
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0.—Effect of Idealized Wind Pressure on 14- Cable Jo Spherical Cop 
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FIG. 11. —Vertical Magnification fone Versus 14-Cable | Spherical Cap 


8 © 22 24 2 2 
12. _—Horizontal Magnification Fac Factor V Versus for 14-Cable Spherical 


(206 Pa). The deflected shape is shown in nin Fig. ‘10. The ssaplivede versus frequency © 
curves are shown in Fi igs. 11 and 12 for the vertical and horizontal ae 
respectively. The oscillating unit load was applied at node 20 (Fig. 10). It -_— 
be seen that some nodes have resonance frequencies lower than that of node 
“Ing general, the response of the structure differs considerably from that without — 
wind load, - the reason being that wind load changes the “pressure distribution 
on the structure as a result of which some cable segments experience higher 
_ pressures, while on others the load is decreased. Thus cable forces are redistribut- — 
ed throughout the structure affecting the natural frequencies of all cable elements, 
though the change i in shape the be insignificant. 


The stiffness | matrix fo rt or uniform pressure- -loaded cable element may . a 
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= | 
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‘obtained small variations of the governing ‘nonlinear equilibrium 
= 
X 
4 H’ = Hi, + Hi, = ty 


2 


Tay are the projections of cable tension 


FIG. 13. 3.—On Derivation of ‘Finite Element Stiffness Matrix id 
AG. 


H at ends J and 2 in the x and y directions, respectively. In general, variation 


ues of X, Y, R, H,,, H H,, and H,, are considered, and 


= Hy + Hyy8Hyy = + 


(27) 


(26) 


i] R (28) 


= 
= 
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‘From 24-26 one obtains 


y(Hyy + Hay) = BH ~ Hx) + 
and from Eqs. 27- 29 it follows that ‘forees; 


: te may now be noted that if only the case Y = 0 is to i aa considered, then 


om Hy — Hy, = - 0, X = L, and Ay = Hay =| pL/2, and from Eqs. : 30 and 


Hy, Y (2) 
21, 22, 23 and the relation H, , = H sin @ are used, then Eq. 32 a 
6(1—esin6) 
¥ in which e = L,/L. To obtain the change in H,, due to an  4-displacoment, R 
only 8Y = = 0i is used and it follows from Eq. 33 that somenytt* 


4 substituting 8X = 0, one obtains = 


_ 


Using ‘Ret: 24-26 to remove 8H, ,. x? and a similar Procedure will 
result in the following influence of -only displacement ement 


=—tan@sY 
the € reciprocal th theorem emas 
8X, 


cos 8 


| 
- 
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sin®@ +cos 


—cos 6 sin 6 


it is shown that the air- supported structures are 
— strongly dependent on internal pressure and static loads, including equivalent — 
static wind pressure. For the relatively small static loads, e.g., snow and dead 
weight, the changes in natural frequencies can be considered as proportional 
to the square root of the internal pressure. The influence of an idealized static 
_ wind-type load is to decrease the fundamental frequency of the structure. 
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= cross-sectional area of cable; 
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complex modulus of elasticity; 


cable segment end forces; 

= finite element stiffness matrix; RA 
"unstretched length of cable Tat 


spanof cable segment; 
‘mass; 


= "internal pressure; pith’ aod ives on 
= radius of curvature of cable saint 
static displace ents; 
= of vibration; 
“local coordinates of ends of cable (i = = 2: 


logarithmic decrement of material; 

i increments of displacements of cable segment e1 ends (i = = 1, 2); 

increments of end forces of cable segment (i= 1,2); | 

angle (Fig. 1); res iret 
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complex angle; an and 4 
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of shear walls « or a braced steel frame in the center of the ‘building, connected 

to exterior columns by flexurally stiff Cantilovers at one or more levels, 
_ Shown diagramatically in Fig. \(a). 

To the best of the author’s knowledge, ‘this was first 

used by Barbacki in 1962 in the design of the 47-story Place Victoria Building 

in Montreal (1). It has since been the basis of design for many other tall buildings 

of a similar size, in both steel and reinforced concrete framing. == 8 


_ When lateral loading acts on the si structure, the outriggers and columns 
hat 


7 would have resulted i ina a free- -standing core Figs. 1(b) and \(c). Outrigger bracing 
is one of the most efficient and economic systems for controlling drift in tall 
“I buildings: the core and columns are obviously necessary components of the © 
_ structure while the additional elements, the outriggers, are usually incorporated 
within the plant levels and, therefore, ey only minimally with the usable 
The magnitude the reduction in drift core “moment. depends on 
relative flexural rigidities of the core, the outriggers, and the columns acting» g 
7 z axially about the centroid of the core. The amounts of reduction also depend = 
on the location of the outriggers within the height of the core. i —t™S 
a Previous studies of the problem have assumed the outriggers to be flexurally : 
a rigid. Taranath (2) investigated one-outrigger structures and showed that the 
optimum location of an outrigger is close to the mid-height of the building, 4 
(0.455 of the total height, from the top). McNabb and Muvdi (3) verified Taranath’s 
result for a one-outrigger structure and showed that the optimum locations a, 
- two outriggers are 0.312 and 0.685 of the total height from the top. ‘Stafford : 
Smith and Nwaka’s study (4) produced generalized results for the optimum 
 'Prof., Dept. of Civ. Engrg. and Applied Mechanics, Mc McGill Univ., MacDonald Engrg. 
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locations in mul in multi-outrigger structures for any number of outriggers. They showed © 
that if the outriggers are located at equidistant height intervals, but excluding, 
one at the top of the structure, a close to optimum reduction in drift is achieved. — 
_ Ineach of the above investigations, the assumption of flexurally rigid outriggers" 
was adopted to simplify the solution of the problem. It was noted, however, 
_ that the flexibility of the outriggers would have an influence on the total drift 
and the core moment as well as on the optimum location of the “outriggers— 
for the maximum reduction in drift. In this paper, the results of an investigation — 
of mullti- with flexible outriggers are presented (5). Generalized 


_ the magnitude of drift and reduction in core moment to be estimated. The 
_eauations developed | for the solution of the optimum locations a are too > complex 


i: a computer solution; therefore, graphical solutions are presented. These allow — 


. rapid manual solutions for multi-outrigger structures, within the constraints of 
the assumptions. 


-PREE CORE 


uF 


“FIG. 1. Outrigger- -Braced Structure; (b) to Loading: and (c) 


n practical structures, the ‘size of cores and columns usually reduce uy up the 
height of the building; thus, the results for a uniform structure, as assumed 
here, cannot be applied exactly. Therefore, the results of the study are useful _ 
primarily as guidelines for estimating member sizes in the preliminary stage 
_ of design. When the exact sizes and locations of the components have << 
. decided, a stiffness computer analysis sh should be used to provide a more p accurate 
‘see. _—The development of the analysis is based on the following 
-_ simplifying assumptions: (1) The structure behaves linear elastically; (2) only ‘ 
axial forces are induced in the columns; (3) the outriggers are rigidly attached 
to the core, , and the core is Ne ade cantilevered | from the base; © the on 


The in of the member ‘Property assumptions for a structure under 


— 
| 


"FREE' CORE B.M. 


Toth tik 


fom, OF he 
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2 2. Dimensions; External Loading Bending Moment 
(ec) Restraining Moment Due to Outrigger 1; (d) Restraining Moment Due to Outrigger — 
(e) Resultant Bending Moment in Core 


using the sizes of core, columns, and outrigger in - the wrens region, will give 
results acceptable for preliminary design purposes. 
Compatibility Equations for Two-Outrigger Structures.—A two-outrigger struc- 
- ture will be used to demonstrate the method of analysis. The analysis of single 
or other multi-outrigger structures can be developed easily from the two-outrigger ; 
case. Compatibility equations” will be developed for ‘the rotation of the 


| will depend upon the material Of construction and the features 
: - of design. Assumption 4 is less restrictive than it might, at first, appear. Although, ie 
ina practical structure, the inertia of the core and sectional areas of the columns ll a 
a reduce the height, the results of concern in a preliminary analysis, i.e., the a 
drift, the core base moment, and the column axial forces are dominated 
s _ by the structural properties in the lowest region. Therefore, a uniform analysis, _ _ 
bar 
q 
7 


. ; and outriggers at the outrigger levels Fig. 2(a). The bending moment diagram 


for the core consists of the free bending moment diagram due to the externally 
load Fig. 2(b) reduced by the outrigger: restraining moments, which 


and 2(d). The r resulting bending moment diagram i is, , therefore, a stepped version 


of the original free moment diagram Fig. 2(e),} 
The core rotations at levels | and 2 Fig. 2(a) can be shown the ~area 


| 
The corresponding expressions for 1 the rotations of the outriggers at locations 
1 and 2 will now be developed to establish compatibility with the core. The 
rotations of the inboard ends of an outrigger where it attaches to the core 
= consists of two components, one due to differential axial deformations of the 
columns and the other due to the ‘outrigger “bending under the action at its 


ends of the column forces. The components of rotation at 1 and 2 due to 
axial deformations of the columns are given by 


) 


while those to the the outrigger are by 


wet 


and Egs. 1b, 2b, and 3b for compatibility of 1 rotations at the level of the lower > 


— 
| 
— 
— 
~ 
) 
- Eqs. la, 2a, and 3a can be used to express the compatibility of the core and = 7 
Outrigger rotations at the level of the upper outrigger, thus — Vim 
> 
q 
| 


It should be noted that (£/), is the effective flex py tigidity < of the outrigger, 


= as though | its length extended from the - column to the centroid of the core. 
(EN), can be me from the actual flexural rigidity of the outrigger (EI’), by —_ | 


GGER ACTUAL OUTRIGGER EFFECTIVE is 

aid 


converting the flexural rigidity of a wide- | beam Fig. the 
rigidity of full- beam Fig. 3(b), 


.. 


full-span beam will be used the paper. Conversion 


from, or reconversion to the actual outrigger —" should be made using 


and M,S( 


The ‘solution of Eqs. 6a and b | give expressions for the 
moments that the outriggers apply to the core & bas > 


0 

q 
[S, + S(H X,)] + M,S(H X,)=——(H?- X3) a) 
which S and 5, are parameters as given 


ax = 
St. + S,SQH X,-X,)+S°(H- 


_ The resulting core moment at any level, at a distance x from the top Fig. 


in which, me<s< i, 0, and for X, 


0. 
7 The forces in the columns due to the bracing | action a1 are + o for 


— 


x 


VALUE OF 


wt 


4 4.—(a) Outrigger tone Levels & in Outrigger Structure; Outrigger d 
Optimum Levels in Two-Outrigger Structure; 
 <x< < X, and + (M, + - M,)/d for x > X,. Then the maximum bending moment 
in an outrigger is the product of the resultant column force at the end of the 
outrigger and the free length of the outrigger, 
_~Drift wee. —The drift at the ¢ top of the structure is obtained conveniently 


— 
| 
— 
i 


‘ 


7 
| 
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4. —Continued: (c) Outrigger Optimum Levels in Three- Outrigger (d) 
Outrigger Optimum Levels in Four-Outrigger Structure ott 
by the second moment-area a method. to the moment Fig. 
_ x? )+. +M x2) 


in which the first term on the right-hand side represents the free cantilever 


q 
id 
. 
7 


deflection under the full external loading, while the second term is the | drift 


reduction, A, due to the outrigger system. 


of optimum of the ‘outriggers to give ‘minimum drift is found 
maximizing the drift reduction, A, the second right-hand side term of Eq. 11. ~~ 
P Differentiating Eq. 11 with respect, first, to . X,, then, to X,, gives two 
simu neous equation 


xX, 


in which dM /aX,, aM, /dX,, , and dM,/dX,, dM,/dX, are the derivatives 4 
of M, and M, from Eqs. 9a and 5b, respectively, with respect to X, and X,. a 
——— substituting for M, and M, and their derivatives from Eqs <<. 
_ Eqs. 12a and 6 can be solved to obtain the optimum locations of the outriggers. 
The solution of Eqs. 12a and b is complex, and a computer must be used. ae 
In the complete expressions for Eqs. 12a and b, the structural properties 
were expressed initially by parameters S and S,, as given by Eqs. 7 and 8, 
_ "respectively. It was found possible to rewrite Eqs. 12a and b in terms of more 
= meaningful dimensionless structural parameters a and B, representing the core-to- 
column and the core-to- -outrigger inertial ratios, respectively, where 
=. 


an 

It was ‘then found possible to simplify Eqs. 12a and b further by combining 
a and into a single parameter w, as given 


a uniform outrigger structure with flexible outriggers. 
12a and b were solved to give the optimum | outrigger levels” 
x, for a range of values of the structural parameter w. These are plotted in 
Fig. 4(b). The restraining moments M, and M, and the top drift A, for a structure : 
with a | characteristic parameter and ‘its outriggers located at ‘the ‘optimum 


of x, and xX, ‘from Fig. 4(b) into Eqs. 9a and 5 and 11, respectively. It should 
be noted that outriggers located for optimum drift reduction do not give a maximum 
core base moment reduction. The core base moment progressively reduces as 


the “outriggers are are placed lower on the str structure. 


4 
| 

) 

| 


_ The solution of multi- -outrigger structures with three, four, or more  outriggers 

(5) leads to problems of equations Ser moments, drift and optimum outrigger 
levels corresponding t to Eqs. 9, and and respectively, for the 


ed equations 
to te written as follows. bas at “pt od) 


Restraining Moments. —These are expressed best in matrix form for sim 


S(H-X,) S,+ S(H- x) ‘S(H X,)  sur- x) 
“aes 


r which 1 n is the pe of outriggers. Eq. 16 requires the structural properties 


and levels of the outriggers X,-X, to be specified. Then the general expression 


In the region ‘between the top and the first outrigger, | the second —— the 
right-hand side of Eq. is taken t to be zero. g2 
= 0). 
Optimum of Outriggers. .—The required outrigger | X, to X 
in an n-outrigger structure, to optimize the reduction in drift are given by the 


M, 


| OUTRIGGER-BRACED STRUCTURES 2009 
Ff 
7 
M, 
x H x 
| 
| | = 
4 | 
| 
.. My; || H? X32 |-2,] x,M, | =0 a) 
M,,, M,,, By Mn H? - x; X,M,,. 


a which M, is the re restraining g moment to. i, and M, i is s the derivative 
q 
“Note est i and j both extend from I-n. Thus, by substituting expressions for + 
the restraining moments from Eq. 16, and their corresponding derivatives, Eq. 
19 can be solved to obtain the optimum outrigger locations X,-X,. These can 
be solved for ranges of values of the two structural parameters a and B, (Eqs. 
13 and 14) « or for a range of the single characteristic parameter w » (Eq. 15). 
Thus, curves can be drawn for a multi-outrigger structure giving the optimu 
& of the outriggers for a range of values of the parameter w. “tie = / 
_ Since the most usual number of outriggers in a structure is from 1-4, curves 7 : 
_ for the optimum outrigger levels in these structures are gven in Figs. 4a) q 


— 


to 4(d), respectively. 
Perror 


; performance of outrigger structures, with respect to variations in the size and 
Structural arrangements of the core, outriggers, and columns. These 


‘include the core-to-column and core- to-outrigger inertia ratios and the number 


Optimum Outrigger Levels. —Observation of the plots of outrigger optimum 

_ levels, Figs. 4(a)-4(d), the following, 

The more flexible the outriggers, ‘the other properties remaining constant, 
higher up the structure their optimum levels 2 are. 


“sensitive the structure is to the outrigger flexibility. 
3. For a constant value of core-to- “outrigger inertia ratio, i. e., B constant, 


decreasing values of column stiffness, i.e., a increasing, cause the optimum 
outrigger levels to move downwards towards the optimum levels for the case 


of flexurally rigid 
= 


Reduction in Core —A useful measure of the influence of the outrig outrigger 
system on the core moment is to express the reduction in the oun base moment 


termed a “moment reduction efficiency,” denoted by %. Fully composite 
behavior occurs when plane sections through the building before bending remain 
plane during bending. In the fully composite case > core moment reduction 


given by 


| 
4 | la 


in which S is given by Eq. 7. Then, expressing ‘the actual moment reduction 
at the base as a a percentage of of the my alte reduction taal ein 


a function - of the 


FIG. 5. —Moment Reduction | Efficiencies for Located Outriggers 


characteristic parameter in Fig. Using Fig. in reverse, the 


base moment can be obtained from the graph by 


In this ante mort reduction in the core moment has been considered as secondary © 
= to the reduction in drift. The, so called, ‘‘optimum outrigger levels” are calculated 


& the core base moment, expressed by Eq. 22 and Fig. 5, are e those corresponding — 


to the optimum outrigger locations. If drift is not the critical consideration, — 
: - additional reductions in core moments at the base and at the intermediate heights _ 


moment reductl re Witt — 
optimally located for maximum drift reduction is given as ® 
| 
| | 


: an outrigger placed close to the base. However, this would contribute little 
to the drift resistance. For outriggers located at other than optimum locations, 


4 the corresponding core moments and top drift can be calculated using Eqs. 


¥ Reduction in Drift.—The drift een of an outrigger : structure can also] 
€ usefully expressed as a ‘‘drift reduction efficiency,” A%. In this case, the _ 
rift reduction due to the outrigger system is expressed as a percentage ul 


the drift reduction that would occur if the system behaved fully compositely. 4 


FIG. 6. —Drift Reducti 


VALUE OF w 4 


Corresponding to to the fully composite moment Bes. 20 and 21, 
the fully composite drift reduction i is by - 


dineuntis as a function of the characteristic - parameter w. The curves can 
_ be used either to assess the efficiency of a proposed ctenctaras system or to 


determine the top. of the system using 


| 
— 
| 


Fig. 6 shows that the larger the ‘number of ot outriggers, “other 
_— constant, the greater the drift reduction. However, the “‘law of i. 
. _ returns” applies in that the increments of increase in efficiency Teduce for 
each additional outrigger. This effect of diminishing return —. that four 
= outriggers is probably the maximum number that can be justified. 


6 also indicates that outrigger ‘flexibility Teduces the drift efficiency 0 of 


The conclusions of th the can be summarized as follows. | 
1. The use of o to mobilize the axial stiffness of columns 
is an efficient way of reducing drift and core moment. a oe 
2 2. The behavior of an idealized multi-outrigger siracture can be aetenie 
s terms of a nondimensional characteristic parameter, w, which is a function = 
of the core-to-column system, and core- -to-outrigger, flexural rigidity ratios. _ 
is 3. ‘In any multi- ~outrigger | structure, there is a set of optimum levels of the . 
-outriggers to give minimum top drift. These optimum levels may be expressed 
in a uniform structure as a function of w. 
; 4. The larger the number of outriggers in a structure, other properties being 
“ —_ the oat the reduction i in drift. However, since each additional outrigger | 


of the outrigger system on drift and core moment. 
The influence of a ange in ‘outrigger enbility is larger for greater inertial 


a 
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wee moment of in inertia of core; 


= = actual moment of inertia of 
effective moment of inertia of outrigger; 


= moment in core at distance x from top; 
restraining moments due to outriggers 1,2,i; 


in core ¢ base moment for fully-composite structure; 


distance from topof structure; 
of outriggers 1, 2 etc from top of structure; 


“a non- 1-dimensional perameter Tepresenting flexural ratio 


non- parameter representing flexural Tigidity ratio 
é 
ol of core- -to- -outrigger; 


att 
duct drift in full te structure; ; 
reduction in top Tit in fully- composi es ructu 


= drift reduction efficiency; rage 


rotation of core or outrigger; and all 
characteristic parameter for outrigger s struc- 


tures with flexible cantilevers = B/ ‘120 + 00 


bee. 
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Willi William R. Spillers," M. ‘ASCE and Sohan an Singh? 


Shape dbliialosaéa (geometric optimization) of structures deals with situations 


in which the shape of a structure is not known a priori, but is to be determined - 
along with other parameters, using some sort of optimization algorithm. Since 
it is most common to optimize other structural parameters ke keeping the shape 
- fixed, the additional difficulty associated v with changing geometry makes | problems 


__ of shape optimization the most complex structural optimization problems encoun- 


7: _ The literature of shape optimization within the finite element method is rather _ 
" sparse. There is, first of all, a branch of it in which node locations do not 7 
but shape varies as some element thicknesses go to zero (13). (This 
uy 
.% is analogous to using an algorithm: to determine an optimal skeletal sebstrecture 


a given an unnecessarily complex initial layout.) Since the node locations are 
‘not determined directly by the analysis, this type of shape optimization is not 
the spirit of the present work. 

the approaches that have been used in continuum-type optimization all 
problems, there is little commonality. These methods Tange from the fersie 
in which Umetani and Hirai “postulate a a bone growth r mechanism based 
upon the existing state of stress to the quasicontinuous (1) approach of Cine 

: ~ and Haug who only resort to the finite element method after performing sensitivity a 
studies based in the calculus of variations. Somewhere in the middle lay the | 
york of Ramakrishnan and Francavilla (7) and the work of Dems and Mroz — 
“ (2) who | optimize specific boundary-related parameters. It is typical of the 
state-of-the-art of structural optimization that no uniform approach to shape 
optimization has yet emerged. it is typical of the state-of-the-art of shape © 


_ optimization that problem complexity obviates all but relatively simple applica- 


_.. present effort derives from the 1971 Ph.D. thesis of Friedland (4), which 


dealt with optimal shape for trusses. Friedland’s work was subsequently extended 
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to other « optimal shape for structures (8-10). The point is 
e that the discretization Process of the finite clement 1 method produces a structural > 


‘This similarity of form is used here as a basis for extending shape optimization 4 
5 “procedures developed for skeletal structures to continuum problems — 
modeled by the finiteelement method. 
_ This paper treats the most simple finite element problem of plane stress using 
4 constant stress triangular” element. It is hoped that the reader will regard 
‘ this effort as an attempt to demonstrate feasibility and forgive—at least for 
the time being—the shortcomings of this particular element. 
_ With regard to optimization itself, the paper adopts a simplification developed - ; 
in earlier work in which the problem is formulated in terms of ‘‘force’” variables | 
neglecting the constitutive equations, i.e +» the plate thickness is written in terms | 
of element stresses, which must, in turn, satisfy node equilibrium equations. 
i: argued earlier, this procedure is most simply explained in terms of plastic 
design, but may also be regarded as an “embedding” of an elastic design problem 
in a larger space that does not contain the constitutive equations. The point 
is that since the optimal design is statically determinate, the satisfaction of 
- constitutive equations can be argued a posteriori. The algorithm itself uses 
Lagrange multiplier method to convert an optimization problem toa problem 
“of solving a nonlinear system; then, Newton’s method is applied. The algorithm 7 
is simple in concept and relatively direct. At every step in the process of iterating, 
q “all the node ¢ coordinates (which are allowed to change under the problem statement) _ 


tan ‘ade 


lane 


It is, first of all, necessary to establish ‘the ‘notation of the finite element 


mato Following Zienkiewicz (14), the simple case of all triangular, plane 


Constitutive Equation o= 


oS which o = stress matrix; € = strain matrix; D = elasticity matrix; B 


shape function matrix; q° = (applied) node force matrix; v“ = element volume | 
matrix; and u = element displacement field. 


epi poh: (14) provides the detail of the finite | element ‘method for the 
case of triangular, constant strain, and plane stress elements as indicated below. 
Element Variables. —Associated with each triangular finite element i area 


stress matrix, o, a strain matrix, €, and an elasticity matrix, D, sothat = 


in which E s Modulus; and v = s Ratio. The ‘Strains 


are related to the element displacement field (u, v) 


a 
) 
} 
| 
a 
| 
| 
q 
| 


a 


yo ny ni big 903 


bw 
Variables.— 


d vector, q. so that 


. 


< 

hape ‘Function Matrix.—The contribution of any element nodes 
m to the equilibrium equations of the system can be written *) 

in which = thickness; and A = element area. the 1 matrix 

Bcanbe partitioned as 


in which B,= | | 

ati 


and, for example, 


+ 


The point i is ; that, for the type of element under discussion, these matrices ‘ 


‘The intention of this section is to set up an optimization problem in which 
the volume of material, written in terms of stresses, is to be minimized while 
_ Satisfying node equilibrium and varying some subset of node coordinates. While 
_ any well behaved stress-volume relationship c or. be used (8), for reasons of 
simplicity the Mises yield condition, 


will be used here to determine the element meaanae given the element stresses. 


2017 

a displacement 
a | 
| 
| 
adi : 
terms of these variables Eq. (7) can be written explicitly as 
| 
a 


id In Eq. 12 f, is the dita onant in uniaxial tension or compression. Furthermore, _ 
it should be noted that only a restricted set of node coordinates may be allowed ~ 


to vary during the optimization process since, if all nodes were allowed to — — 


move, the loads — to the supports giving an optimal design of zero 
optimization can now be formulated to find | unknown node 
and the element thicknesses which 


Miaimize ya 


volume). 


a For it is to a at 

: this point. As discussed earlier, the approach presented in this paper is motivated — 
by the truss problem where it has been quite successful. In the truss problem, — 
the ‘‘force”’ variables used are in fact bar forces. For the finite element method, _ 
the obvious choice of * ‘force’ variables, the member stress, leads to numerical 


n of th In; an attempt 
arrive at variables with the same as the truss problem, the 


That is, new variables F, N, K, A ‘that are analogous to the variobies 


Bb € are introduced. In terms of these new | new variables, , the optimization problem 


Optimality Conditions. —Proceeding. in a classical manner, it is now possible 
to replace the above optimization problem by a nonlinear system in the following 


manner. First of all the Lagrangian of Eqs. 17-19isformedas = ~— | 


in which 8 is a matrix of me mukigiers, and r is s the : set of inti 

“that i is allowed to vary. The Kuhn-Tucker (optimality) conditions are then derived “f 


21 


| 

J 

| 
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FINITE ELEMENT EXAMPLE 


It has been convenient in the past to solve systems, such as Eqs. . © 23, 
rl in two steps, first fixing the geometry and computing element thickness, then — 
\ fixing the thicknesses and computing new node coordinates. The case of fixed _ 
_ geometry has been described elsewhere (3), and it will only be noted here that 
‘z to optimize in this case simply requires iterating with any standard finite element * 
cy ‘program using (Poisson’s ratio) v = 1/2 and at each step using the yield condition a 
_ For the step in which the thicknesses are fixed and new node conedinates 
are to be computed, Eq. 23 is simply linearized in a Newtonian manner and 7 
ia used to compute these changes i in | <opeeanins. Some of the details of this: sep 
some from nonlinear it convenient to rewrite Eq. 
5= or F=0 (24) 
which “operator ‘refers. to the unknown node coordinates and 
_ Kis the geometric stiffness matrix. Weing Newton’s method, Eq. 24 is linearized 


= iti is, therefore, n necessary to ‘compute the gradient 1 matrix of both Vf and K 5 ~ 
in order to compute the changes in coordinates, dr, 
_ Some Gradient Terms.—As indicated in Eq. 20, the application of Newton & 
_method to the optimality criterion of Eq. 23 requires that some — 
¥ of finite element terms be computes. These Gexivetives are discussed below. oT 


= | 
Lol respect to the x and y coordinates of the 3 nodes can be written simply — ; i 
terms of the b’s and c’s of the finite element method 
Since f is defined by Eq. (9) it follows that 


From definition of the b’s and c’s (Eq. 10) it follows that 


from which the Hessian matrix of the function, f, follows 's directly. no ad at 1 he 7 


‘The geometric stiffness matrix and its gradient are somewhat more difficult — 7 
to deal with. It is first of all convenient to define a new set of coefficients Pico 


/2 1/2 
=c,A- 


/2 


Second also be computed as bis 


2 


It follows now from Eq. | 


= vector form . returns again when Kei is multiplied by da Sas 


— 
(30) 
Ox, Vv | 
| orthat may be partitioned as a A 
| 


FINITE ELEMENT EXAMPLE 
.+c%, F,)8,, + (6%, F,+ cf, Fo), 


_ It is clearly possible to take the gradient of Eq. 25 by Am there i is no 
convenient form to put it in, no further derivatives will be indicated here. - 
af Following Ref. 3, again two applications of the above optimization procedure | 
are presented here (Fig. 1). The first is simply a fixed ended plate (beam) 
_ whose lcwer right hand node is loaded vertically. In the second example, a 


hole which is allowed to move is added to the first — 


— 


ex 


Several « comments concerning numerical perfo performance are in order. First of 


_ all, there is the common problem that the step size, which emerges from Newton’s 

method, is sometimes too large. It is, however, a relatively simple matter to 
reduce the step size ¢ automatically when the new design does not have a reduced 
weight. A more difficult problem has to do with the aspect ratios of the finite - 
elements. As can be seen from the figures at the end of this paper, the aspect -_ 
ratio of some finite elements can degenerate quickly. The solution to this problem 
should be to lay out a new mesh when this degeneracy occurs. That has ie 
been done here, and is left t to fu future efforts. ‘Tables: and 2 indicate 


| 
| 
= 
| 
la 
@g “to the extent that step size difficulties increase with increasing iterations. Ong__ ; 
— _ ‘might expect this effect to disappear is relayout is used at every itecaticn. 


With to shapes that have evolved under (Fig. 2), 
it would be theoretically possible to construct a Michell truss to transmit the — 
ae efficiently to the supports (5). This truss would follow two << 
families of lines forming a Henky-Prandtl net. While paper is concerne< 
with finite element layout rather than truss — it should be noted the — pf ‘ 


Step size Volume, 


for step 2 cubic inches ubic inches 


0.05 6.830 
TABLE 2.—No Hole Case 


Volume, Volume, in 
cubic inches inches Ad) 
7 


step 
63360 
| 


final comment on graphics is in for the case 
fixed geometry, it is almost impossible to understand thickness distributions _ 

= without some sort of three-dimensional graphic presentation. When node locations _ 
are also allowed to change, it becomes doubly difficult to understand and, thus, lke 
monitor design generated by an algorithm. It is then, in fact, _ almost impossible Fe 

_ to deal properly with problems of this kind without on-line computer graphics _ 

&s displays. Fig. 2 attempts to convey something of the spirit of the er 


graphics environment but space limitations make that difficult i the ‘print 
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Conciupinc Remarks 
methods = 
Pro “Venkayya (12). These methods largely to fill the 
~ vacuum left by ‘direct methods of mathematical programming that, as a rule, * 
do not work when applied to large systems (6). As an adjunct, these methods _ 
_ work in structural terms; theorems about these methods are structural theorems. 
> _ From one point of view, the approach used here is direct and even say 
= plastic design problem is taken in its most simple terms and treated directly he 
using the finite element method and Lagrange multipliers. This approach certainly _ 
raises (and does not settle) the question of the relative efficiency of moving 
: possible nodes at each step versus perturbing only the boundary nodes. 
= _ Finally, this paper is clearly exploratory in nature. It will probably be some 
_ time before basic issues of more ¢ sophisticated finite element 1 models, , grid size 


of structural optimization a1 are settled i in the. context of the finite element method. 7 
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‘TECHNICAL Notes” 


yy To provide a ~ _— ASCE for publication of technical ideas. that have not t advanced, 
- yet, to the point where they warrant publication as a Proceedings paper in a Journal, the 
publication of Technical Notes was euthodand ~—s the Board of Direction on October 16-18, 


1967, under the following guidelines: 


and Professional Publications, ASCE, 345 East 47th Street, New York, N. Y., 10017, poe 
with a request by the author that it be considered as a Technical Note. satan opal til 

_ 2. The two copies will be sent to an appropriate Technical Division or Council for review. 

If the Division or Council approves the for publication, it be 


4, :- The technical publications staff will prepare the material for use it in the earliest poniite 


issue of the Journal, after proper coordination with the author. 
: _ §. Each Technical Note is not to exceed 4 pages in the Journal. As an approximation, — 
each i manuscript page of text, tables, or figures is the equivalent of one-half a Journal 

4 6. The Technical Notes will be grouped in a special section of each Journal. = o 
7. Information retrieval abstracts and key words will be unnecessary for Technical Notes. 

Technical Notes will not be included in 


«10. Technical Notes will be included in ASCE’s annual and cumulative ‘subject and ‘author 


‘The manuscripts for Technical Notes must meet the following requiremen’ 


" 1. Titles must have a heagh not exceeding 50 characters and spaces. wr 
_ 2. The author’s full name, Society membership grade, and a footnote reference stating present 
employment must appear on the first Page « of the manuscript. Authors need ‘not be Society 
3. The manuscript is to be submitted as an original copy (with two duplicates) that is typed — 
double-spaced on one side of 8-1 /2-in. (220-mm) by |1-in. (280-mm) white bond paper. 
4. All mathematics must be typewritten and special symbols must be properly identified. — 
‘The letter symbols used must be defined where they first appear, in figures or text, and arranged 
_ §. Standard definitions and symbols must be used. Reference must be made | to the lists 
_ published by the American National Standards Institute and to the Authors’ Guide to the 
7 6. Tables must be typed double-spaced (an original ribbon copy and two duplicate copies) 
; on one side of 8-1/2-in. (220-mm) by Il-in. (280-mm) paper. An explanation of each table 
appearim the text, 
“- 7. Figures must be drawn in black ink on one side of 8-1 /2-in. (220-mm) by 11-in. (280-mm) 
oe. Because figures will be reproduced with a width of between 3 in. (76 mm) to 4-1/2 
in. (110 mm), the lettering must be large enough to be legible at this width. Photographs must 
be submitted as glossy prints. Explanations and descriptions: must be made within the text 


8. References cited in text must be be typed at the end of the Technical Note in alphabetical 


9. Dual units, i.e., U.S. os followed by SI (International System) units in npeeeameene, 
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DETERMINATE TENSEGRIC SHELLS 
At tensegric shell consists of a net of cables and bars connected at both 
7 of their ends to the cable net noes, no bar is connected to the other bars. — 
Because of the lack of a continuous supporting rigid element the tensegric shell a 
_ is loose and does not have a definite shape before erection. The shell is stabilized 
to the required shape by prestressing. Prestressing, of the tensegric shell can ; 


+ 
‘some cables. The prestressing gives the shell its shape and introduces tension _ 


in the cables and compression in the bars. When the shell is in its prestressed - 
space no bar touches the other bars. By disconnecting some bars or cables 

the shell is disassembled and can be transferred to another place for re- -erection. 

_ A model of a tensegric shell has been presented already (2). 
: Tensegric shells were ‘developed from tensegric bodies (1), fundamentally the 

_ tensegric shells are similar to the well known one skin pneumatic shells. In 
the case of a pneumatic shell the nonpermeable envelope of the shell is stabilized | 


= the shell gets its proper shape when it is prestressed by pressurized air, 
whereas in the tensegric shell the cable net, which is related to the pneumatic 4 


- shell envelope, is stabilized and the shell gets its proper shape when it is prestressed 

by lengthening the properly arranged bars connected to the net nodes. The 
advantage of the tensegric shell is that no special envelope and no constant — 
pumping is required and there are no problems considering openings as in the a 


In this. work the behavior ofa determinate tensegric shell when it it is is prestressed 


t the first stage of designing a tensegric shell the nature of the tensegric © 


net to be used in the structure must be fixed. Various nets and possible — 


arrangements of bars suitable for tensegric shells were described earlier (2). — 
_ Of particular interest are the regular tensegric shells in which at each node 
4 the same number of cables are joined to the bar. In regular determinate tensegric 
_ shells, where the number of bars, cables, and the reactions along the shell 
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a written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on January 15, 
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edge is equal to three times the number of nodes, there are five cables joined & 
to the bar at each node. An example of a net to be used in determinate tensegric | 
After the nature of the net has been fixed it must be fitted to the geometry 
4g the specific shell under design and the length of the different bars and cables 


to 


must be found appropriately. In Fig. 2 the tensegric net shown in Fig. | is 
_ fitted into a spherical determinate tensegric shell. This is the initial geometrical 
configuration of the tensegric shell. Because the shell has no continuous rigid 
supporting elements, constructing a tensegric shell with bars and cables of the 
length analyzed from the initial geometry leads to a loose structure without 
a definite shape which is of no use. In order to stabilize the structure and 
to give it a definite shape the structure must be ce ie iigubal 
it 


_ 


| 
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ee FIG. 1.—Net and Bars Arrangement in Case of Determinate Tensegric Shell = g 
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some of the cables. The “prestressing changes. the geometry of the tensegric 
_ Shell from the initial geometry to the prestressed geometry and introduces inner 
7 forces to the cables and bars of the structure. Only tensegric shells in which 
prestressing induces tensile forces in the cables and compression in the bars 


; 7 are suitable. A tensegric shell achieves a stable geometry when equilibrium 


_- between the inner forces i in “the cables and bars and the reactions along the 
shell edges due to the prestressing at all the nodes is established. The inner 7 
forces at every node of the shell should satisfy equilibrium of forces in the — 


equilibrium equations from all the shell nodes in a matrice form yield oo 


‘in which Pa = the vector of the > inner forces in in the cables and bars | of the 
-tensegric shell and the ‘reactions along the shell edges; and A = the matrice 
related to the geometrical configuration of the prestressed tensegric shell. It 
is convenient to assume that during the prestressing the geometry of the tensegric 
Shell is changed from the initial geometry to the Prestressed geometry. A node 


Matrice A is a a function of x,, u,, and 1 w, of all 
- Shell nodes. Eq. | has nontrivial solutions only if the rank of A is less than = 
the number of unknowns, in the case of a determinate pe shell A is 
a determinant and implies 


in which A° = the matrice consid 

shell and it is a function of a and z, of the shell nodes pe A’ = the : 
linear matrice of the nodal displacements u,, v,, and w,. Considering Eq. 3 


andi ignoring second order terms Eq. 2 takes the folowing form: — 


z which cs = the co-factor i, j in determinant A°; and a, = the element i,j 
“of A’ matrice. Eq. 4 indicates a linear relationship between the nodal displace- 
_ ments, u,, v,, and w,. The displacements of the shell nodes during the prestressing 
are restricted by the fact that the length of some cables and bars is not changed 
- during the prestressing. If the length of a member k connecting nodes i and — 
d is not changed during the prestressing; the displacement of nodes i and j, 

in the case of small nodal displacements have to 
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a In the case where the nodal displacements are small, matrice A can be written q 
= 
— 
q Ul spi 
It seems most convenient to impose some restrictions on the distortion of 


seeking for a minimum general distortion which implies that: —" 
+ v2 + w4) = Min ad af atittes allt 


HOTS 


the summation in Eq. 6 is on all the moveable shell nodes. el 


cnuitidiie the linear Eqs. 5 and 6, Eq. 6 yields to an appropriate set of linear a 


equations from which the nodal due to the 


| forces, in the bars and cables of ‘the shell an 


_£, applied on it can be found by using, Eq.’ 


-_of the fact that the geometry of the tensegric shell was achieved by prestressing, 
_ only a few loading cases are related to a consistent set of the linear equations 
given by Eq. 7. In the general loading case this set of the linear equations 
is not consistent; a fact which indicates that the tensegric shell cannot sustain | 
the external loads having the geometry achieved by the prestressing. In these 
loading cases the geometry of the shell changes ‘until equilibrium i is established 


7 in which B = the matrice related to the new geometry of the tensegric shell. ‘ 
It is assumed that during the change in the geometry of the shell a node, i, . 

located at Xs Ves and z, in the shell prestressed geometry has been moved 7 

- z u,, v,, and w,. In the case of a small nodal displacement the inner — 


, in member k connecting nodes | i and ji is given by Eq. 9: ee ee a 


= =p. + fis fi=— > ,) Kv, 1 v 


‘in which p, = the inner prestressing force; E, and A = the elasticity” and 
cross-section; and 7 = the length of member ki in the prestressed > Fon 


In the | case of a small nodal displacement £ B takes the form of: ae “Al aed 


in which =a linear function of the ‘nodal displacements. 


Eqs. 10 and and ignoring second order effects takes 
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The method of proper arrangement of the bars and cables in the case of 
a determinate tensegric shell were described. It was shown that to stabilize 

The condition for a stable tensegric shell was ‘formulated and the or 
of analysing the shape and the inner forces of the prestressed tensegric shell 
bY It was shown that most types of external loads when applied on the determinate — => 
tensegric shell will cause a considerable geometrical distortion of the shell. 
The method of analyzing these distortions i is considered i in the work. —_ 


—Revenences 
‘Emmerich, D. G., “Morphology « et ‘Aujourdhui, ol. 160, 
2. . Vilnay, De. , “Structures made of infinite regular ten ric nets,” ud JASS Bulletin No. — 


Vol. XVIII-1, ‘Apr., 1977, pp. 


which is a set of linear equations of the nodal displacements of the temsegric 
— ¢€ prestressing is to the extent that proper tensile forces are left in the J 
cables under the expected external loads. q 
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anyone who has significant comments or questions regarding the content of the paper / technical 
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ments, on the free vertical vibration of suspension bridges. The first statement _ 

4 of the paper talks about the First Tacoma Suspension Bridge which failed, 

- in 1940, by the first antisymmetric torsional mode of vibration. Actually, “a 
great majority of authors solve the suspension bridge as a plane structure loaded _ 
symmetrically about its longitudinal center line. However, a few authors (12, 13,14) 7 

a solve the suspension bridge as a three-dimensional structure considering the ‘ 


torsional stiffness of the deck and making allowance for asymmetric loadings. | 
_ Although the torsional stiffness of the bridge deck is taken into consideration 
in his analysis (Eq. 11), the author did not try to offer any analytical or even 
experimental results of torsional vibration, as if the structure is a plane structure _ 
- with no torsional stiffness. In other words, the author did not make use of | 
the torsional stiffness of the bridge and his numerical example does not show 
value. Rather, the author. presented a separate paper (3) for free torsional 


in which (f/1) is the sag-span ratio of the cable. 

However, W, defined by the author should be the ‘‘average’’ one 

of the cable per unit span length (and Lana per unit cable length) : | since | the Jf 
cable profile is parabolic under dead load. 
2. In his sixth assumption (on p. 2054) the author states that the suspenders . | 
- s ‘are assumed to remain vertical during the vibration,’ * while the example (Fig. 

l(a)] shows some inclination, 

ow The ‘integral of Eq. 4 is given | by T (15) 


while it is given we the writer (13) in a much simpler form a wa fot 
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sec not cos as stated by the author on p. 2054). 
Also the integral of f Eq. 6 is given by the writer as 


x = 


4. Eq. may ‘alternatively be w written as 


Ix 
in which d’v/dx? = curvature of the bridge deck = ck = M/EI: d/dx (dv/dx) 
12. Fukuda, T., Proc. Japan - 
Society of Civil Engineers, No. 242, Tokyo, ‘Japan, Oct., 1975. 
ms 13. Hegab, H.1.A., “‘Statics and Dynamics of Suspension Bridges with Notes on Suspension ~ 
vail Cables and Nets,”’ thesis presented to the University of Tasmania, at Hobart, Australia, 
in 1977, in partial fulfillment of the requirements for the degree of Doctor of Philosophy. q 
14. Sih, N., ‘‘Torsion Analysis for Suspension Bridges,’’ Journal of the Structural Division , 
"ASCE, Vol. 83, No. ST6, Proc. Paper 1431, Nov., 1957, 1431-1-1431-81. 
15. Timoshenko, S., “‘The Stiffness of Suspension Bridges,” Transactions, ASCE, Vol. 


Errects or Fire on REINFORCED CONCRETE MEMBERS sit 


4 Discussion by I. D. B Bennetts” and LR. Thomas* nas” 
_ This paper has raised a number of important considerations with regard to 
the prediction of the behavior of concrete structures in fire. The authors have 
ably presented the advantages of analytical methods of estimating fire resistance 
the use of traditional fire resistance tests. 
os The various approaches which have been used to estimate the fire resistance 
of concrete members are summarized in Table 2 and have been numbered from — 
ts one to ten in the Table. Approach | is the standard fire test. Approaches 2, 
_ “November, 1980, by B. Ellingwood and J. R. Shavers (Proc. Paper 15802), —t™ 
Research Officer, The Broken Hill Proprietary Co., Ltd., Melbourne Research Labs., 
273 Wellington Rd., Clayton, Victoria 3168, Australia. 
“Sr. Research Officer, The Broken Hill Co., Ltd., Melbourne Research 
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DISCUSSION 


‘he remaining» 
approaches (5-8) a are re analytical methods chad sonaesa the behavior ofa — 
in real fire conditions. To do this they combine prediction of the fire environment 

> with the response of the member in that environment. These methods es 
_ require information on the fire compartment (e.g., fire load, ventilation, emissiv- 

“ The modelling 0 of the structural Tesponse ranges from simple flexural strength 


approaches (2-5, and 10) to relationships for steel 
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"TEMPERATURE Aa 
We believe that the authors are correct in suggesting that care should be 
—_—- in using some of the simpler design methods, as it is possible that 


temperatures on the concrete compression zone. Other problems s which could 
negate the validity of some of these simpler methods are those of spalling, — : 
_ interaction of shear and bending, (25), bond failure (27), and longitudinal splitting 
J along the beam web. This latter problem has been noted by European workers 
is to be due to the differential thermal expansion across the 


Tempera compression failure may occur due to the influence of elevated — 


— 


: 3, and 4 are semianalytical methods based on experimentally determined tempera-_ 
; = ture distributions. Approaches 9 and 10 are analytical methods which attempt os 
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FIG. 13.—Differential Thermal Exp of Concrete Beams 
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4 Factors relevant t __ Approach Number 


Standard time temper- 
ature curve in stan- 
dard furnace 
environment 
Fire load and ventila- 
tion requirements of 


Thermal conductivity 
concrete and steel} 
_ Specific heat capacity 
of concrete (and 
Experimentally deter- 
_ mined isotherms as 
function of time 


=. 


Failure criterion 


Member failure 
 terion—failure to’ 

satisfy equilibrium 

Member failure 
terion—maximum 


deflection or deflec- 
tion rate 


— 


Required mechanical characteristics 


_ Fire resistance test of 

full size member 
Critical temperature of 

steel (based on pre- 
vious assessments 

using 8 and 9) } 

Reinforcement: 
Influence of temper- 
ature on proof or 
ultimate stress 
4 Influence of temper- 


ature on stress- 


(instantaneous) 
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Thermal 
characteristics 
Thermal 
behavior 
Influence of t temper- 
ta ature on concrete 
crushing strength 
Influence of temper- 
a ature on concrete 
Stress- strain rela- 
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tionship (instan- 
taneous) pun 
expansion 

characteristics 
of full constitu- 
tive equations tak- 
ing into account 
Stress history, 

time and tempera- | 
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webs of concrete beams (Fig. 13) with the resulting stress field causing the 
_ cracking. If the simpler flexural strength approaches are to be used then safeguards — 
must ¢ exist t to to ensure that alternative failure 1 mechanisms e not occur. Some — 
progress along these lines has already been made. 
_ . As the authors have pointed out, fire resistance tests are expensive. In addition, 
b: it is only practical to test the simplest of structures as it is difficult to simulate 
the real boundary conditions. On the other hand, full scale fire tests indicate 
the level of importance of spalling, shear, or bond failure, and furthermore _ 
4 such tests do not require detailed knowledge of thermal or heat transfer properties. : 
| Although there are limitations and problems with the standard fire test, it 
seems to the writers, that there are also difficulties with the use of analytical 
methods, particularly the more sophisticated methods. The calculation of temper- 7 
atures throughout the member cross section requires the use of reliable thermal — _ 
E properties. Iti is well known that i it is + difficult to determine experimentally the 


is variable (Fig. 14) and appears to ) be not only a function of concrete type : 
but also testing technique (17). If the concrete under consideration is prone “ 
to a . or if the concrete cover tends to crack and fall away as the sumed 
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of the difficulties m mentioned above for heat flow ana ow concrete 
members can be overcome by ‘‘calibrating”’ the theoretical model (Approaches - 
9 and 10) against the results of fire tests (Approach 1). Thus, one can choose . 
thermal properties for the concrete which seem to result in good correlation 
between and theoretically determined temperatures. Tt tis difficult 
r 
4 A piindens with the use of analytical techniques to predict temperatures within 
a member is the availability of sufficient data to properly model the heat transfer 
. across the boundary of the member. This depends on both convection and a 
‘Tadiation: and the writers question the accuracy with which the = 


_ The use of the most sophisticated structural analysis technique ( (Approach - 

- 8): relies on empirically based constitutive relationships for both steel and concrete : 

(18). The determination of such relationships for concrete involves a very large 
e experimental programme and the writers question the applicability of these 


constitutive relationships to other concretes. s. Approaches 5- 10 assume that perfect 
_ bond exists between concrete and steel and do not predict the member shear- 
‘moment or shear strength. For some structures therefore, care must be exercised 
using these methods to predict fire resistance, 
_ The engineering goal must be to predict the response of structures in the 
"teal fire environment. This combines prediction of the thermal environment, 
prediction of the ‘thermal response of the structure to that environment, and 
prediction of the structural behavior. This is the intention of methods 5-8. 
The writers would appreciate the authors comments on how methods 5-8 can 
be validated with regard to both the thermal and structural response in the 
real fire situation. Does sufficient data exist to allow this to be achieved? 
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‘Space Tru Truss S STupIEs WITH ForcE LIMITING Devices 


by E. A. Smith,? M. ASCE al 
- The authors present ‘some interesting experimental results that can increase’ 
- our knowledge about space truss behavior and give some insights into the 
difficulties and potentialities of Force Limiting Devices. This writer has also” 
been studying space truss behavior and has a number of comments and questions 
_ This writer agrees with the authors’ observations about the effect of brittle- -type } 
- post. buckling characteristics of transition length struts on the behavior of the _ 
total system, which can also be of brittle-type. Lack of fit will prestress the 
_ members, , some favorably, and advance or delay the achievement of the member 
ri peak load, the buckling load, thus causing the premature initiation of failure — 
of the total system. However the authors’ results seemed to indicate that even 
i accounting for lack of fit, none of the monitored members achieved the buckling 
) loads obtained in the member tests. Kahn and Hanson have presented results _ 
- (15) on tests of mild steel transition length struts that also fail to reach the 
This writer believes that one reason for this phenomenon is the sensitivity 
_ of the peak load of the member to slight variations in the pre- and post- -buckling — 
_ behavior. Rosen and Schmit have shown (16) that the effective axial stiffness 
of a pre- -bent member is givenby = 
tom 
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er; ¢, = the maximum nll 
initial eccentricity; r = the radius of gyration; om p = the ratio of wed applied 
load to the buckling load, 
Schmidt has derived a similar expression for loaded 

_ members. ‘The stiffness begins to decrease rapidly when p is greater than about 


0. 8. If inelasticity occurs, then the effective buckling load decreases" causing P 


- additional increase in p and decrease in stiffness. The effective axial stiffness 

therefore is quite sensitive to the effect of the initial eccentricities as the peak 

load is approached. Because of the steepness of the post-buckling response 
aq curve, the peak load at the intersection of the pre and post buckling responses 
‘ botham and Hanson have showed (14) that it is possible to considerably reduce an 
the percent deviation of the maximum compressive bend from the analytical =e J 
solution of a member by accounting for initial curvatures, ‘ 
_ Although the authors’ stiffening of the joints would reduce the effect of 
jackknifing, it would not eliminate initial eccentricities that can lead to the 
above described imperfection : sensitivity. In addition, it should be noted that 


~ buckling of adjacent members. For example, in the final test, #10, members 


the stiffening of joints did not appear to allow unbuckled members to attenuate 7 


9, 1, 10, and 3 all failed at the same external load cho only 9 and 10° 


4 were rigidly connected to lower stressed members. 


| 


_ While the authors’ claim to have tested a ‘ “statically determinate”” subassem- 
blage, without the assumption of ‘symmetry, the assemblage is clearly both 
externally and internally indeterminate. If the reaction points are unequally stiff, 
or if initial slack existed in the supports, prestress will occur. How did the 
authors support the structure to ensure that all the reactions were equal and 
thus not contributing to prestress and do they have any initial force measurements 
for the determinate Such data _ be = 


14. Higginbotham, A. B., and Hanson, R. D., “Axial Hysteretic Behavior of Steel 
Masten, ”* Journal of the Structural Division, ASCE, Vol. 102, No. ST7, Proc. Paper 
15. Kahn, L. F., and Hanson, R. D., “Inelastic Cycles of Axially Loaded Steel Members,” 


ies of the tga Division, ASCE, Vol. 102, No. ST5, Proc. Paper 12111, 


(16. Rosen, A., and Schmit, L. A., “Design Oriented Analysis of Imperfect Truss 

International Journal for Numerical Methods in 


| 
} 
— 
| 


Couumn B BASE with AxiAL Loaps AND Mowests* 


The be congratulated for this very fine study an and 

_ to the understanding of the effects of axial and undirectional bending moments 
on the concrete, base plates, and properly designed anchor bolts; using two 
stress grades of wen bolts of various diameters, and various degrees of 

a Due to the type of -“through- -bolt” anchor, of course, the: effect of an 
-end- anchorage failure could not be produced. Do the authors know at approxi- 
‘mately what percent of bolt working loads the bond break-down occurred? — 
nee less credence should be given to anchor bolt bond allowance. 

_ Two directional (or diagonal) bending was not included. Are the two preferred 7 
— of analysis practical or indeed safe approaches for this problem? ‘te 

Did the authors consider the possible benefit, eg. , greater base fixity, that 
al occur from use of thicker base plates? The writer believes that the plate 
deformations, then, would not be too inconsistent with the “‘strain compatibility _ 
approach,”’ (section at bottom of base plate remaining plane after paren 
An aniivels of the 13 test conditions show that (in theory) the edge pressure 


¢@ the concrete will never be greater than 0.56 f/, cf which is well within oe 


permissible of A.LS.C. Specifications, when using. strain compatibility. 
While it is generally true that failure load, based on anchor bolt achieving its’ 
full working tension, is decreased about 63% (see author’s conclusion number 
_ 3) however, for test specimens 3 and 4 values of P,,, increased to 13.9 kips : 
and 10.7 kips, or 20% and 87% increments respectively (1 kip = 4.45 kN) 
when compared to values in Table 3, column 3. 
; Perhaps an additional measure to provide greater base fixity of columns — 
be to keep anchor bolts shorter while providing end-anchorage of bolt to carry 
_ _ Cannot strain compatability be considered consistent within the stress range — 
working loads, as opposed to failure loads? doc 
‘The authors’ method « of measuring concrete bearing area after the test seems 
infallible. However, as. as they indicate, the use of reinforcement in the menor 
should further reduce the chance of concrete failures. de 
= Finally, does the large concrete stress of specimen 4, Table 6, represent — 
_a “*‘break-down’’ of the ultimate strength method for that range of small mee 
bolt size and large eccentricity? or 46° echalaue. would by 
yy The writer looks forward to the possibilities of more rigid guidelines to column 2 
_ base design, but would also plead for a further look at the practicality of the 


strain compatibility concept in such areas as maximum base fixity and biaxial 
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How ca an we achieve a designer s relationship between A.I.S.C. 
— 
ome then of die the 
ON LATERAL BUCKLING OF BeaMs* 


Discussion by Narayanan Subramanian’ 


_ ‘The writer would like to appreciate the authors for writing an interesting 
a paper, which provides ‘information about the lateral buckling of beams with = 

ome. However, the writer a like to make the following observations. 


1. The authors muna the use of the formula suggested by Timoshenko for 


"predicting the lateral torsional buckling load of the cantilever beam, after making — 
allowances for the openings. ‘This formula, Eq. applicable. only in 
: elastic range, , ie. & , the formula is applicable only if the quantity 17/AL is very 
small (5). The authors have chosen the beams, such that the value of 17/AL 
ranges between 7.6 x 10 ~* and 1 x 10~*. But in practice, the value of 7? /hL 


may be greater than those considered by the authors and thus the buckling 


in Ea ‘should be re ad “Sab. 


i a very narrow rectangular cross section in which t/h = a small quantity. (The 
authors have considered beams having ¢/h values ranging from 0.06- -0.08). 

Theoretically buckling may occur also in the case where t/h is not very ‘small 

but where the length L is very large. In this case a large deflection in the 

plane of the web will be produced before lateral buckling occurs, and this 

-— Geflection should be considered in the derivation of the differential equations — 
of equilibrium. It has been shown that in these cases 3 the constant value of in 


b/h. Dinnik and Federhofer have’ shown that for values of b/h equal t to o 1/10, 
: /5 and 1/3, the value of +y will be equal to 4.085, 4.324, and 5.030, respectively 
_ 4, The formula given in Eq. | is applicable only when the load, P, is applied 
at the centroid of the cross section. It has been shown that the loads applied | 
at the level of the top flange at points free to twist represent a particularly a 
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DISCUSSION 
severe condition for ©). It will be useful to euusider this case 
_ since in many practical cases, the loads will be applied at the top flange. por oe es 


6. Nethercot, D D.A., ‘‘The Effective Lengths s of Cantilevers Governed by Lateral Buckling,”” 
Structural Engineer, Vol. 51, No. 5, May, 1973, pp. 161-168. 


"the behavior of a tall structure ‘subject te to wind-induced lock-in excitation. The 
"writer has obtained some of the results outlined by the authors, e.g., an increase 
in the reduced velocity results in a change in the nature of the crosswind excitation. 
_ At lower velocities the forcing function exhibits randomness; however, as the | ae 
reduced velocity is increased, 
random n noise. Furthermore, ‘studies conducted by the writer have indicated = 
that there is an apparent decrease in the net damping, as the reduced velocity — 
is increased. It would be interesting to learn if the authors tried to quantify 
change in the net damping in their study. 
_ The authors have arranged the equation of motion such that all the motion — 
y cua terms are placed on the right hand side of the equation. ‘This has x 
resulted in different spectral curves as the reduced velocity fis. increased (Fig. : 
4) . In Eq. 2, if appropriate adjustments for the change in the net damping ~ i 
= of the system at various reduced velocities are made, all the curves in — 
} — . 4 may collapse on one curve. A simplified expression for the amplitude © 
aerodynamic as a function of reduced velocity, can 


aerodynamic damping, even when the motion induced effects are 
_ dominant. A similar approach has been reported in Refs. 22 and 23. Pid e chy 
= The model proposed by the authors is essentially based on a forced vibration 

: theory. The uniqueness of the forcing function is not necessarily insured, since — 

it is obtained by an inverse approach. An accurate assessment of the system “oq 

properties is foremost for the success of this technique. It would be v very beneficial a, 
for the readers to know whether the model was free to oscillate in the two 
: ~ principal axes of motion or the alongwind direction was restricted. Experiments | 
conducted by the writer have indicated that the behavior of a model is quite 
i January, 1980, by Kenny C. S. Kwok and William H. Melbourne (Proc. Paper 15976). 
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"different i in these configurations. results have that 
the redistribution of energy, imparted to the structure by the wind-structure 
| See from one principal axes to another has the - tendency to limit the _ 
response amplitude, thereby, imparting indirect damping to the system. This 
Be nea is more pronounced for a structure with equal natural frequencies — 


is s important for practical purposes. force coefficients a are highly 

_ dependent upon Reynolds number, aspect ratio, surface roughness, turbulence — 
intensity and shear in the approach flow. A wide variation of forces can be 
measured on apparently identical cylinders in different test facilities. One of 

_ the major difficulties is the simulation of high Reynolds number flow around 
wind tunnel” models “of circular crosssection. Augmenting the model surface 
with discrete or homogeneous, or both, roughness elements may contribute 
towards simulation of the features of high Reynolds number flow. In this manner, 

a wake may be developed artificially, at subcritical Reynolds numbers, which - 

is similar to the one naturally produced at the fullscale Reynolds number. A 
- literature survey indicates that when cylinder oscillates, the flow field is tonal 


from that of a stationary cylinder and may have less dependence on the Reynolds 


number. In this manner small scale models of cylinderical structures could be 


employed to reproduce the dynamic behavior; nevertheless, the steady drag : 
_ forces would be dependent on the Reynolds number. The authors have provided 

a good basic procedure, using small scale models, which can be further semen 

toa account nt for an number ‘simulation. The 


Arrenowx. 


Kareem, A., “Wind Excited Motion of Buildings, thesis presented to Colorado State 
= University, at Fort Collins, Colo., in 1978, in partial fulfillment of the requirements 
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DistRIBUTION AND DouBLe Skin WALL” 


Discussion by Theodore Theodore Stathopoulos,’ M. M. ASCE 
_ The author refers to the significance of the wind load distribution for the 
design of double skin wall panel constructions. He emphasizes the lack of relevant — 
research information and indicates the large discrepancy between cosine based 
“February, 1981, by Raymond M. L. Ting (Proc. Paper 16055), 
* Asst. Prof. of Engrg., Centre for Building Studies, Concordia Univ., 1455 de Maison- 
neuve Blvd. W., Montreal, Quebec, Canada H3G IMB. 
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| procedure outlined by the authors. 3) 4 


on commonly made seeaiae with respect to the percentage of wind load 
acting on each component of the double skin wall system. cuuey a 


_ A significant parameter for the wind loads to be reasonably estimated is 


the permeability of both facia and liner components of the wall structure. The | 
_ Statement in the paper that ‘‘the liner construction is reasonably airtight, while - 
the facia construction is not airtight” gives perhaps. a preliminary —— 
but is not enough. Also under question is the assumption of static wind loading; — 
the natural wind is always unsteady, 
_ Based on his own measurements of both external and internal pressures acting 
‘on low- -rise building (3), the writer has that wind pressures on 
“pressures sures acting: on the inside surface of the - wall “depend drastically on the 
wall porosity. This is true of mean and fluctuating components of both external — 
and internal pressures. The implication is that the total wind load [(external 
" pressure — internal pressure) x area] acting on the wall is a function of the | 
permeability of structure and the wall openings. 


1 of the small internal volume of ‘the: spacing between facia and liner but in 
principle the problem remains the same. Pressure or suction on each side of 4 
_ the facia has to be found for the calculation of the design wind load. — erautlis 

Pe Finally it is worth mentioning that the problem of wind load ae 


3. -Stathopoulos, T., Di, Davenport, G., “Internal Pressure Characteristics 
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RESPONSE or RC un UNDER ‘Grounp Motions" 


Medhi Saiidi,* A. M. ASCE end 
‘The authors should be commended for studying a real structure, at = 
Y 2 in the laboratory environment, to evaluate their timely analytical model. The 
writer would like to draw attention to the basic assumptions made in the analysis — 
and the procedure used to evaluate the calculated responses. Newt 
* _ The authors took into account several characteristics of reinforced concrete r 
to model the structure ‘‘realistically.”’ However, two important factors which 
_ generally reduce the “apparent stiffness’’ of reinforced concrete systems were 
not considered. One of these factors is stiffness degradation effects — 7 
in cyclic testing of concrete specimens. Experimental studies have shown that — 
. concrete specimens loaded past the compressive strength exhibit reduction in 
’ stiffness upon unloading and reloading (Ref. 32). Exclusion of this factor in 
the analysis [Fig. l(a)] results in an average : structural stiffness ss larger call 
The other important effect not included in the analysis is slippage of rein- 
- forcement caused by partial loss of bond. Deformations caused by bond alip + 
- usually constitute a significant portion of the total deformation. Ngo and Scordelis - 
developed a linkage element to simulate bond slip in their finite element study — , 
of reinforced concrete beams (Ref. 30). Similar to stiffness degradation of | 
concrete, bond slippage results in softening of the structure. 
_ The fact that the model was “‘stiffer’’ than the actual test structure can be 
well observed in Table 2. The calculated maximum displacements are smaller 
than those measured, while the predicted base shear is larger than the measured 
value. Perhaps more . detailed explanation would be possible if the en 
response histories were reported in the paper. A. 
The writer realizes that the inclusion of stiffness ness degradation and vali dig 
effects further complicates the analysis. However, these effects need to be 
= into account either explicitly or implicitly, if a detailed complex analysis, 
whatever its use, is to be successful. To obtain bry estimates of response, _ 


of the equivalent ‘ rr -degree’’ oscillator are listed in “Table 3. The assumed 
deflected shape is the same with the fundamental mode shape of the wall with 


_ cracked moment of inertia used for the first story. The primary force- displacement | 


curve was calculated based on a static analysis of the wall subjected to lateral 
loads proportional with the height, applied at floor levels. The ‘‘yield’’ ‘oat 


_ “February, 1981, by Awadh B. Agrawal, Leslie G. Jaeger Aftab ‘ Mufti (Proc. 
“Asst. ‘Prof. of Civ. Engrg. , Univ. of Nevada, Reno, Nev. 89557. _ 
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“for inthe hysteresis model. 


- -Tepresents | the load and diaglacement corresponding ‘to the y yielding of the wall 


at the base. The post-yielding branch was defined based on the loads causing 
cf a maximum strain of 0.003 in concrete. This strain was used only to define | 
the slope, and not to limit forces. od Jsbow baa” 


Assumed deflected shape norm normalized with respect to top 


Base moment at the “‘yield’’ point a3 kip-in. 
of ‘post branch of the primary curve | 95 kip-in. fin 


wee 4.—Maximum Measured and Calculated Responses 


Qmodel 
0.23 


shear, in kips 4.000 


moment, in kips-inch Ph 188 


/_Note: 1 in. = 25.4 mm; | kip = 4.45 kN; | kip-in. = OR - 
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10- 10.— —Top Level Response of Shear Wall 

_ The base acceleration and the measured top-level displacement fempeunt were 
obtained from a crude. ; digitization of the waveforms presented in Figs. 6 and © 
7. Although high- frequency accelerations could not be digitized accurately, they 3 


Note: 1 in. = 25.4 mm; | kip = 4.45 KN; | kip-in. = 113.0N-m. 
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were considered insignificant because their contribution to the seaponss ° was 
4 Fig. 10 represents the top-level response of the shear wall based on the Q-Model 
and the authors’ model. It can be seen that the waveforms from both models 
were generally similar. With respect to the maximum responses (Table 4), the 
results from the Q-Model were either closer to the experimental data or 
ee: the same with the authors’ results, even a the oonee Model ¥ 


i calculated responses for 1.5 sec and refer to measured response histories a 
presented in Ref. 14 without reciting them. In addition, the calculated crock 7 
_ pattern for the first simulated earthquake is compared with the observed creck 

_ pattern after the third simulated earthquake. The writer believes that, to evaluate i 
a complex nonlinear model, duration of the analysis should be sufficiently long 
to include both small-amplitude and large-amplitude responses, because some | 
of the assumed characteristics of the model may influence these responses 
differently. With respect to crack patterns, it is inappropriate to compare the 
theoretical result for the first earthquake run with the experimental result for 


the third run, because cracks are generally rtre from one motion ced a 


Neo, D., and A. C. Scordelis, “Finite Element. Analysis. of Reinforced Concrete 
_ Beams,” ACI Journal, Vol. 64, No. 3, Mar., 1967, pp. 152-163. 
31. Saiidi, M., and M. A. . Sozen, “Simple and Complex Models for Nonlinear Seismic 
oe = of Reinforced Concrete Structures,’’ Civil Engineering Studies, SRS 465, 
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SPEED DisTRIBUTIONS AND AND RELIABILITY Estimates" 


1005, Table 2, heading: Should read instead of 


Page 1005, Table 2, heading: Should read = 0. of =1.10° 
a Page 1005, Table 2, column ‘ ‘Distribution,”’ lines | and 5: Should re: ae ud “Rayleigh 
Page 1005, Table 2, column ‘ “Distribution,” ’ lines 3 and 7: Should 7 _— 


‘ 
| 
| 
> 
&g 
) q 
| 
23). 


: 
> 
| > 
: 


